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The  paper  describes  a tentative  approach  the  problem  of  hull- 
equipment  interaction  in  an  under  water  shock  environment,  by 
the  way  of  mechanical  impedance  techniques.  Two  different  methods 
are  used  to  analyse  point  impedance  of  submarine  hulls,  by  sinusoidal 
and  transient  excitation.  Measurements  give  a spring-like  behaviour 
of  the  hull  in  the  frequency  range  of  interest  for  rigidly  fixed  equip- 
ments. Underwater  slack  experiments  on  a shock  test  section  have 
confirmed  this  result.  However  this  description  is  no*  sufficient.  The 
hypothesis  of  a viscous  damper-like  behaviour  of  the  hull  at  higher 
f requencies  is  invoked  to  account  for  the  phenomena  that  are  observed 
during  the  first  while  of  the  shock  response. 


1 - INTRODUCTION 

A research  program  has  been  undertaken 
these  last  years  by  the  French  Navy  to  improve 
its  knowledge  of  the  shock  environment  of  sub- 
marine equipments.  The  purpose  is  to  work  up 
better  methods  for  shock  resistance  design  and 
particularly  to  obtain  better  shock  input  values. 
Such  a process  has  been  made  possible  by  impro- 
vements in  shock  measurements  techniques  and 
in  analysis  of  the  dynamic  behaviour  of  structures. 

The  amount  of  work  devoted  to  these  pro- 
blems has  not  permitted  to  develop  an  extensive 
program  of  full  scale  experiments  and  thus  has 
induced  us  to  concentrate  on  a more  fundamental 
aspect  of  the  dynamic  interaction  phenomena  bet- 
ween the  hull  and  the  supported  equipments. 

Though  we  are  conscious,  from  the  experience  of 
other  Navies,  that  this  procedure  is  quite  less 
sure,  nevertheless  we  hope  that  we  shall  "et  a 
very  useful  understanding  of  the  basic  principles. 
Obviously  the  risk  is  eventually  to  get  unable  to 
process  a to  much  sophisticated  information. 

Our  first  approach  of  the  interaction  pheno- 
mena is  concerned  with  the  case  of  medium  weight 
equipments  (up  to  a few  tons)  fixed  to  the  hull 
through  sufficiently  stiif  supports  (rubber  mounts 
being  excluded). 


The  problem  is  considered  under  follo- 
wing aspects  : 

- Impedance  characterization  of  the  hull  at  the 
fixture  points, 

- Modal  identification  of  f xcd-base  equipment, 

- Mechanical  coupling  between  the  two  connec- 
ted parts. 

- Determination  of  shock  input  signals  in  unloa- 
ded hull  condition. 

The  study  is  conducted  in  an  experimen- 
tal way  by  various  dynamical  analysis  tech- 
niques on  one  hand,  and  by  underwater  explo- 
sions in  a full  scale  submarine  test  section  on 
the  other. 

To  validate  the  hypothesis  of  our  ap- 
proach resulting  from  dynamical  analysis  of 
hull  and  equipment,  we  have  extracted  the  major 
part  of  the  useful  information  directly  from  the 
time  history  of  the  signals  recorded  during  shock 
tests.  So  we  have  been  brought  to  improve  signi- 
ficantly our  various  filtering  methods  in  the  fre- 
quency and  also  in  the  space  domains  and  we  have 
consequently  fully  realized  up  to  what  point  the 
philosophy  of  some  shock  measurements  specia- 
lists of  US  NAVY  was  pertinent. 
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The  work  is  now  in  progress.  First 
results  have  made  it  possible  to  point  out  some 
simple  ideas  that  we  are  going  to  present  ; 
many  questions  remain  yet. 


II  - POINT  IMPEDANCE  OF  THE  HULL 

Point  impedance  measurements  have 
been  made  on  the  hull  of  the  submarine  section. 
This  experimental  facility  consists  in  the  full 
scale  section  of  a conventional  submarine  (dia- 
meter : 5 m,  length  : 16  m).  The  hull  was  im- 
mersed up  to  the  level  of  the  hatch  and  the  flo- 
tation tanks  were  full  of  water. 

Two  types  of  excitation  of  the  hull  have 
been  used  : 

- Sinusoidal  force  signal  provided  by  an  electro- 
dynamic shaker  of  limited  strengh  (200  N) 
(Figure  1). 


Fir.  la  - Sinusoidal  Excitation  of  Hull 


Fig.  lb  - Sinusoidal  Excitation 
of  a Rigidly  fixed  Mass 

- Shock  impulse  force  supplied  by  the  percusion 
of  an  anvil  spring  device  fastened  to  the  hull, 
by  a guided  drop  mass  (Figure  2).  In  this  case, 
the  impedance  is  obtained  by  a FOURIER  ana- 


lysis of  the  transient  signals  : force  and  acce- 
leration. Results  are  displayed  on  following 
diagrams  for  the  case  of  an  exciting  force  nor- 
mal to  rib,  45°  downward,  in  the  middle  of  the 
immersed  section. 


Fig.  2 - Transient  Excitation  of  Hull 


- Figure  3 shows  the  force,  acceleration  and 
displacement  signals  delivered  by  piezoelec- 
tric force  and  acceleration  transducers.  It 
immediately  appears  that  the  displacement  is 
roughly  proportional  to  the  force.  So  in  the 
frequency  range  of  interest  the  hull  has  a 
spring  like  behaviour.  After  the  first  pulse, 
a low  frequency  oscillation  (a  few  Hz)  is  occu- 
ring  in  the  displacement  signal.  It  tends  to 
confuse  with  the  erroneous  integration  of  the 
piezoelectric  transducer  noise.  That  low  fre- 
quency behaviour  has  been  analysed  with  other 
kinds  of  transducers  : piezoresistive  and  servo 
accelerometers. 

The  FOURIER  analysis  of  the  signals 
has  allowed  to  assess  the  accuracy  oi  these  first 
conclusions.  On  Figure  4a  we  can  see  the  com- 
pliance function  computed  with  an  FFT  algorithm 
of  1024  points.  This  compliance  is  rather  com- 
plex in  its  details  but  the  spring-like  behaviour 
is  easily  pointed  out  again  at  medium  frequency. 
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Pig.  3 - Shock  excitation  of  Hull 
Tine  signals 


Fig.  4 - Point  Impedance  of  Hull 
Shock  and  Sinusoidal  Excitation 


The  Figure  4b  gives  the  compliance  func- 
tion measured  at  a similar  point  under  a sinusoi- 
dal excitation  by  the  usual  techniques  of  mecha- 
nical impedance.  The  results  are  in  agreement 
with  the  preceding  ones.  We  notice  a resonance 
tendency  around  the  hoop  frequency  of  the  hull. 
Above  that  frequency,  the  compliance  can  be 
considered  as  resulting  from  the  cumulative 
effects  of  a lot  of  higher  frequencies  modes. 

The  average  behaviour  is  to  be  related  to  the 
point  impedance  of  an  equivalent  infinite  plate 
which  is  known  as  being  real.  So  the  global 
behaviour  of  the  hull  on  that  range  of  frequency 
may  be  roughly  described  by  a viscous  damper. 

An  additional  verification  of  the  actual 
existence  of  the  medium  frequency  stiffness  has 
been  carried  out  by  loading  the  hull  with  a ri- 
gidly secured  mass  of  2. 5 t and  by  measuring 
the  resonance  frequencies  of  the  system. 

Results  have  credited  the  value  of  the  stiffness 
previously  estimated  (Fig.  5). 

From  another  point  of  view,  we  have 
found  out,  by  measuring  a^ain  the  compliance 
of  the  hull  in  the  vicinity  of  the  fixture  points 
of  the  mass,  that  this  function  is  quite  similar 
to  that  ia  unloaded  conditions  (Fig.  5).  So  the 
input  impedance  of  the  hull  may  be  considered 
as  being  quite  independant  of  the  near  supported 
equipments. 


Fig.  5 - Point  Impedance  of  Hull 
Influence  of  Loading  Mass 
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Presently,  similar  tests  are  undertaken 
in  different  types  of  French  submarines.  Measu- 
rements are  extended  to  the  low  frequency  range 
with  percussion  excitation. 


ffl  - MODAL  IDENTIFICATION  OF  THE  FIXED 
BASE  EQUIPMENT 

The  choice  of  the  fixed  base  modal  des- 
cription of  equipments  has  been  made  by  the  US 
Navy  many  years  ago.  We  have  adopted  this 
point  of  view.  The  main  advantage  in  our  sense 
is  that  this  description  fits  well  with  the  calcu- 
lation of  the  response  to  a given  input  shock  mo- 
tion of  translation  of  the  base.  It  is  suited  also 
to  formulate  the  equipment-hull  dynamical  in- 
teraction on  the  basis  of  the  measured  hull  com- 
pliance matrix. 

In  the  case  of  a base  of  limited  dimen- 
sions, the  modal  fixed  base  parameters  may  be 
determined  experimentally.  The  base  is  mate- 
rialized by  the  table  of  an  electrodynamical 
shaker,  the  equipment  being  secured  to  it  through 
force  transducers,  the  signals  oi  which  are  sum- 
med up  vo  get  the  resultant  force.  The  modal  pa- 
rameters : frequency,  modal  mass,  damping 
factor  are  obtained  from  the  mass  apparent  func- 
tion by  an  identification  process  (Fig.  6a). 
Transfer  functions  displacement/force  and 
st  rain/for ct>  are  measured  simultaneously  in 
order  to  give  the  modal  shapes  and  modal  parti- 
cular strains  by  the  same  separation  procedure. 

Such  a technique  succeeded  when  applied 
to  light-weight  models  of  structures  (Fig.  6b) 
and  can  be  extended  easily  to  medium  weight 
equipments  (1  t).  Presently,  we  are  rather  orien- 
ted to  perform  the  same  analysis  by  shock  tech- 
niques by  fixing  the  equipment  on  the  anvil  of  a 
shock  machine  and  making  a FOURIER  analysis 
of  the  transient  wave  to  get  the  related  transfer 
functions. 

Medium  weight  shock  machines  are  so 
able  to  accomodate  equipments  of  up  to  3 t 
weight.  For  heavier  ones  it  seems  to  us  that  the 
measurements  must  be  made  in  situ,  by  point 
impedance  of  the  equipment  fixed  to  a submarine 
hull,  the  fixed  base  hypothesis  being  of  quite  less 
a value  in  such  a case. 

An  important  limit  to  these  general  prin- 
ciples of  description  and  analysis  of  the  dynamics 
of  equipment  is  their  non  linear  behaviour  which 
we  often  met  with  at  high  level  of  excitation. 
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Fig. 6b 


Fig.  6 - Fixed  base  impedance  measurement 
on  Electrodynamical  shaker 


IV  - INTERACTION  WITH  HEAVY  OSCILLATORS 
1)  Coupling  with  a single  oscillator 

As  the  spring-like  behaviour  of  the  hull 
corresponds  to  the  frequency  range  of  fixed  base 
first  resonances  of  equipments  we  are  dealing 
with,  we  consider  in  a first  step  the  simplified 
following  linear  model  of  interaction  (Fig.  7a) 
with  an  SDOF  oscillator. 

8o(t),Vo(t)  are  then  the  shock  motions 
that  can  be  measured  at  the  fixture  points  in 
unloaded  conditions.  With  a loading  oscillator 
of  mass  M and  natural  frequency  W,  , the 
FOURIER  transform  of  the  acceleration  ^ 
is  : 


ff(w) 


lfo(W) 
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Fig.  7a 


LOADED  SHOCK  SPECTRUM 


Fig.  7 - Coupling  with  Heavy  SDOF 
Spring  Effect 


The  right-hand  member  of  the  equation  is 
the  response  of  an  unloading  oscillator  of  fre- 
quency Q . . The  loaded  case  can  be  deduced 
from  it  by  the  frequency  transform  (!)  , With  a 
unit  step  velocity  V0  H(t)  the  acceleration  shock 
spectrum  is  : (Fig.  7b). 

WC9)-  np=f  flLV/z 

\Q,2  + 1 j 

For  a practical  use  with  a definite  hull 
stiffness  K,  we  obtain  a Design  Shock  Spectrum 
to  the  following  example  (Fig.  8). 


DESIGN  SHOCK  SPECTRUM 


(I09) 


Fig.  8 - Desigr  Shock  Spectrum 
for  Underwater  Explosion 


At  low  frequencies,  the  spectrum  design 
value  is  equal  to  V0  . At  high  frequencies,  the 
limit  acceleration  value  is  Vo \TK. 


The  shapes  and  values  obtained  in  that 
manner  are  close  to  the  shock  standards  used 
by  Navies  of  several  nations. 


If  we  admit  that  for  higher  frequencies 
the  impedance  of  the  hull  tends  to  be  equivalent 
to  that  of  a viscous  damper  of  coefficient  C 
(GJ>-js-),  the  model  turns  to  be  : (Fig.  9a). 


The  transfer  function  is  : 

ffM = y„M  m -1- 

1 +J-C*  W 


With  a velocity  step  input  V^iti  the  shock 
spectrum  of  the  loaded  acceleration  has  the 
following  shape  (Fig.  9b). 


This  simplified  description  of  the  two 
principal  effects  of  the  hull  impedance  shows 
how  the  Design  Shock  Spectrum  may  be  modi- 
fied according  to  the  relative  values  of  the  cha- 
racteristic frequencies  \fK  and  • , compa- 
red with  the  transition  frequency  . When 
these  phenomena  are  not  sufficiently  separated, 
it  seems  better  to  compute  directly  the  res- 
ponse of  the  system  with  the  complete  model. 
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LOAOEO  SHOCK  SPECT3UM 
VELOCITY  STEP  INPUT 


Fig.  9 - Coupling  with  Heavy  SDOF 
Damping  Effect 


2)  Coupling  with  a multi-degree  of  freedom 
system 

In  a simplified  approach,  we  consider  a 
system  with  a predominant  mode  loading  the 
hull,  the  other  ones  having  negligible  reaction 
effects  at  the  fixture  lerel.  We  consider  then 
the  shock  spectrum  of  the  acceleration  . 

For  a velocity  step  V0  H(t)  the  velocity 
V,  (t)  has  the  following  shape  : (Fig.  10a). 

V.(e).V0H(e)|uQ02cos(fi0  o)j 

Qo  s — ^1 — and  0 = t \/Tj[ 

1 + iCii 

(0) = V0  H (e)  [i  (9)  (1  - Q *)  + Q03sin  (Q0  0)] 


We  are  induced  to  define  the  stock  spec- 
trum of  fl#(e)  . 


A conservative  evaluation  of  max(Q„) 
consists  in  the  sum  of  the  maximal  absolute 
values  of  the  two  terms  in  the  right-hand  mem- 
ber of  0 . The  related  shock  spectra  are  given 
on  figure  10b. 


Fig.  10  - Shock  Spectrum  at  fixture  point 
of  a loading  SDOF 


We  can  notice  from  it  that  : 

-if  Qn  ^ £2o  the  same  velocity  step  spectrum 
is  preserved, 

- if  Qn  » &o  the  oscillator  is  excited  at  reson 
nance, 

• if  the  shock  spectrum  value  is  that 

defined  for  the  preyaling  mode, 

- if  the  spectrum  tends  toward  a re- 

duced velocity  step  spectrum  : 

v„(i  - n3  -v.(_ L_J 

li*  Q?/ 


For  either  case  : Qn  <fi0orQn^Q,it  is 
possible  under  certain  conditions  to  take  into 
account  the  loading  effects  of  secondary  modes 
on  the  source  impedance  {stiffness  K and  stif- 
fness K in  parallel  with  stiffness  M LJ,2  , by 
a procedure  similar  to  that  described  before. 

In  the  most  general  case,  the  interac- 
tions between  the  modes  through  the  hull  impe- 
dance must  be  taken  into  account,  it  is  then 
necessary  to  set  in  action  analytical  procedures 
in  order  to  obtain  modal  responses. 

For  this  purpose,  we  have  choosen  the 
FOURIER  technique.  The  transfer  functions  of 
the  whole  system  are  calculated  and  give  the 
possibility  to  obtain  the  responses  to  a given 
input  signal  by  FF1  Algorithm. 

The  figure  11  shows  the  velocity  res- 
ponse of  the  three  modes  system  (Fig.  1 la)  to 
the  "vertical  rise  triangular"  velocity  input. 
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Fig. 11b 


Fir.  II  - Velocity  response  of  a three-modes 
System  by  FFT  algorithm 


V - EXPERIMENTS  IN  SHOCK  ENVIRONMENT 

1)  Measurements 

The  problem  of  shock  measurements  is 
known  to  be  a difficult  one.  The  difficulties  are 
due  to  the  severity  of  the  severity  of  the  trans- 
ducers mechanical  environment,  and  to  the  need 
to  discriminate  the  useful  part  of  the  informa- 
tion among  predominant  erroneous  or  non  signi- 
ficant signals.  So  the  measurement  system  must 
be  organized  specifically  for  each  specific  case 
with  a complementary  use  of  the  following  trans 
transducers  and  signal  processing  techniques  : 

- Piezoresistive  accelerometers  rather  than 
piezoelectric  ones  in  order  to  avoid  shift  and 
low  frequency  cut  off  problems, 

- Mechanical  filters  to  protect  the  transducer 
and  to  be  adapted  to  the  frequency  range  of 
interest, 

- Spatial  filtering  by  mixing  of  several  trans- 
ducers to  point  out  some  particular  ensemble 
effects, 

- Use  of  strain  gauges,  and  force  transducers 
at  fixture,  to  identify  the  shock  input  motion 
in  its  effects, 

- Time  history  signal  processing  with  filtering 
techniques  to  separate  modal  responses  and 
with  integration  to  obtain  velocity  information, 

- Shock  spectrum  analysis  of  the  signal  pre- 
viously conditioned  by  the  mentioned  method. 

2)  Unloaded  shock  characteristics 

Unloaded  shock  spectra  have  been  ac- 
quired by  various  experiments  with  underwater 
explosions  in  our  shock  test  section.  They 
assume  a constant  velocity  shape  in  frequency 
range  from  a few  Hz  to  a few  kHz  as  it  is  usual- 
ly accepted. 

The  analysis  of  the  integrated  accelera- 
tion signal  has  confirmed  that  fact  and  led  us  to 
admit  that  the  shock  input  is  an  approximate 
velocity  step  with  a short  rise  time  (<0.5  ms ) 
and  a long  decay  time  (Fig.  8). 

3)  Experimental  investigation  of  coupling  effect 

Experiments  have  been  mad®  with  a 
heavy  loading  mass  (a  few  tonnes),  and  various 
fixture  conditions.  We  describe  one  of  them 
which  is  quite  significant  of  the  hull  behaviour. 

A 2. 5 t hass  is  fastened  to  the  hull  through 
a very  stiff  seating  (Fig.  12a).  According  to  the 
presented  model  of  the  hull,  the  equivalent  dia- 
gram is  (Fig.  12b). 
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Acceleration  and  velocity  measured  at 
the  top  of  the  mass  are  plotted  on  Figure  12c. 
The  velocity  is  oscillating  around  the  step  input 
Vo  with  two  periods  relative  to  two  succes- 
sive phases  of  the  phenomena. 


Fig.  12  -Underwater  Explosion 
Experimental  Results  on  2.5  t Mass  (Rigid  Mount) 


In  the  first  milliseconds,  the  velocity 
step  is  directly  transferred  through  the  equiva- 
lent viscous  damper  C to  the  oscillating  system: 
M,  16K  and  the  resulting  maximal  acceleration 

Then  the  mass  oscillates  on  the  two 
springs  in  series  (#  K)  and  the  acceleration  is 

k~"v°Vpr* 

Thus  we  are  induced  to  talce  into  account 
successively  two  aspects  of  the  impedance  which 
bring  about  two  quite  different  values  of  the  ac- 
celeration in  opposite  directions. 

The  viscous  damper  is  considered  as 
active  during  the  little  while  when  the  hull  acts 
like  an  infinite  plate. 

When  the  boundar  y effects  take  place  due 
to  reflection  and  propagation  around  the  cylin- 
der, the  behaviour  of  the  hull  turns  to  be  that 
of  a simple  spring. 

With  this  hypothesis,  the  duration  of  the 
damper  effect  might  be  estimated  to  be  equal  to 


a characteristic  propagation  delay  of  the  elastic 
waves  of  the  order  of  the  hoop  mode  period. 

Such  an  experiment  gives  another  possi- 
bility to  validate  the  measured  values  of  K et  C. 

This  case  can  be  reasoned  by  using  ener- 
gy transfer  concepts.  In  the  first  phase,  the 
energy^Mvi^  transferred  from  the  stif- 

fness 16K  to  the  mass.  In  the  second  phase  this 
kinetic  energy  is  transferred  back  from  the 
mass  to  the  stiffness  K :£MV|siU££3  (Fig.  12). 

^ l\ 

From  accelerations  and  it  is  then 
possible  to  deduce  the  reaction  forces  at  the  fix- 
ture point. 

This  way  of  reasoning  is  suited  to  non 
linear  stiffness  behaviour  as  in  the  following 
experiment  (Fig.  13a). 

The  same  mass  is  fastened  to  its  seating 
with  Belleville  springs.  The  acceleration  and 
velocity  signals  (Fig.  13b)  show  the  successive 
actions  of  variable  stiffness  in  series.  Each 
time  the  acceleration  Jf(t)  (or  the  transmitted 
force)  is  controlled  by  the  most  flexible  spring 
in  operation. 


AS:  Scaling 

SC,  BaHcvllla  Springs 

CD,  Hull 

OE  , Belts  plasticity 


Fig  13b 

Fig.  13  - Underwater  Explosion 
Experimental  Results 
on  2.5  t Mass  (Non  Linear  Mount) 
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Energy  transfer  calculation  makes  it 
possible  to  access  the  stress  level  which  is 
reached  in  every  part  of  the  fixture  system  in- 
cluding the  region  of  plasticity. 


VI  - CONCLUSION 

The  results  that  have  been  described 
are  concerned  with  a limited  number  of  experi- 
mental configurations.  However  they  already 
allow  us  to  obtain  a first  view  of  the  hull-equi- 
pment behaviour,  in  a certain  range  of  mass 
and  resonant  frequencies. 

Mechanical  impedance  techniques  have 
shown  a spring-like  behaviour  of  the  point  impe- 
dance of  the  hull  in  the  frequency  range  of  inte- 
rest for  rigidly  fixed  equipments.  Underwater 
shock  experiements  on  a shock  test  section  have 
confirmed  this  result. 

This  first  view  incites  us  to  deal  with 
the  time  input  velocity  signals  as  basic  data  for 
analysing  the  interaction  phenomena.  The  resul- 
ting Design  Shock  Spectra  are  close  to  the  values 
and  shapes  commonly  used. 

However,  this  description  is  not  wholly 
convenient.  The  hypothesis  of  a viscous  damper- 
like behaviour  of  the  hull  at  higher  frequencies 
is  invoked  to  account  for  the  phenomena  that  are 
observed  during  the  first  while  of  the  shock  res- 
ponse. 


The  damping  effect  has  been  related  to 
the  caracteristicimpecanceof  an  equivalent  in- 
finite plate.  So , that  tempory  effect  exists  until 
reactive  behaviour  takes  place,  after  propaga- 
tion and  reflections  of  elastic  waves. 

This  interpretation  is  quite  incomplete 
and  does  not  account  for  the  effect  of  water  on 
imped  mce  and  water-cavitation.  These  pheno- 
mena would  have  to  be  investigated.  Simulations 
on  hull  models  under  water  and  in  air  may  bring 
out  the  predominant  influence  of  hull  or  water  on 
the  source  impedance,  and  ascertain  the  origin 
of  its  stiffness-damper  aspect. 
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The  Navy  large  Floating  Shock  Platform 


Pari.  1:  Physical  Description  and  '’apebilities 


C.G.  Schrader 
West  Coast  Shock  Facility 
S.in  Francisco,  California 


Equipment  to  be  installed  on  Navy  ships  must  be  tested  and/or 
designed  to  meet  Navy  shock  specification  standards.  The  previous 
capability  for  direct  testing  extended  to  items  weighing  40,000/ 
60,000  lb.  The  Large  Floating  Shock  Platform  has  now  been  con- 
structed and  calibrated,  and  its  shock  characteristics  are  compat- 
ible with  those  of  the  previously  existing  test  devices.  Direct- 
testing  capability  for  shipboard  shock  can  now  v»  extended  to 
items  weighing  c-rw;  400,000  lb. 


BACKGROUND 

The  Nan/  requires  that  equipments  intended 
for  shipboard  use  be  shock  tested,  if  possible, 
to  the  specifications  of  MIL-S-901.  The  test- 
ing machines  specified  in  the  current  version 
of  MIL-S-901  are  capable  of  accomodating  test 
loads  of  up  to  60,000  lb. 

In  1969  the  West  Coast  Shock  Facility, 

San  Francisco,  California  was  tasked  by  the 
Naval  Research  Laboratory,  Washington,  D.C.  to 
conduct  a conceptual  study  for  a Large  Floating 
Shock  Platform  (LFSP)  which  would  be  capable  of 
shock  testing  loads  of  60,000  to  320,000  lb. 

This  study  indicated  that  if  the  LFSP  were  to 
be  designed  as  a scaled-up  version  of  the  Float- 
ing Shock  Platfo-ra  (FSP),  some  economies  would 
have  to  be  madi  These  were  accomplished  by 
changing  the  ate  .0.  from  HY  80  to  I ITS  and  bottom 
structure  from  an  inner  deck  to  open  framing. 

In  1971,  WCSF  was  tasked  by  Naval  Ship 
Systems  Command  (as  part  of  Operation  DIVE 
UNDER,  Phase  II)  to  prepare  contract  plans  and 
specifications.  The  LFSP  was  built  in  accor- 
dance with  these  and  was  delivered  to  WCSF  in 
February  1973. 

DESCRIPTION 

The  LFSP  is  a rectangular,  flat-bottcraed 
barge  50  ft  long  by  30  ft  ceam  (Fig.  1).  The 
sides  and  ends  are  shock-resistant  to  a height 
of  12  ft  above  the  botten,  and  are  continued 
a further  6 ft  in  height  by  a lighter  bulwark. 
The  shock-resistant  bottom,  sides  and  ends  are 
constructed  of  1-1/2"  IIT3  plate  strongly  re- 
inforced with  gussets  and  stiffeners:  the 
bulwarks  are  3/8"  NTS  plate.  The  bettor;  3 ft 
of  the  LFSP  is  occupied  by  an  egg-crate 
arrangement  of  longitudinal  and  athwartship 


Fig.  1 - The  Large  Floating  Shock  Platform 

stiffeners  with  cell  size  of  approximately  4 
ft  square.  The  upper  edges  of  these  stiffener! 
are  capped  by  flanges  3 in.  wide  x 3 in.  thick, 
forming  a mounting  plane  to  which  the  test  ioad 
is  xttnrhed  by  a suitable  foundation  structure. 
The  LFSP  is  covered  by  a three  section  canopy 
18  ft  high  at  its  center.  The  maximum  overall 
dimensions  of  test  items  which  can  be  accommo- 
dated are  46  ft  x 26  ft  x 32  ft  high,  although 
for  higher  items  a special  cover  could  be  fit- 
ted instead  of  the  standard  canopy.  Empty, 
the  I FSP  draws  62  in;  with  an  added  load  of 
500,000  lb,  the  draft  would  increase  to  about 
124",  leaving  a freeboard  of  20  in.  on  the 
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shock-resistant  sides.  With  allowance  for 
fixturing,  this  implies  that  test  items  weigh- 
ing up  to  400,000  lb  can  be  handled  provided 
the  weight  distribution  of  the  installation 
satisfies  stability  requirements. 

Figure  2 shows  the  Floating  Shock  Platform 
(FSP)  installed  in  the  LFSP  as  a calibration 
test  load.  The  FSP  is  28  ft  x 16  i*.  x 10  ft 
high. 


Fig.  2 - FSP  installed  as  test  load 


Figure  3 shows  a close-up  of  one  of  the 
three  mounting  plates  forming  the  foundation 
for  the  FSP. 


The  shock  environment  on  board  the  LFSP 
is  developed  by  detonating  a 300  lb  charge  at 
a depth  of  20  ft  a selected  distance  free  it. 
The  size  and  depth  of  the  charge,  ana  the  loca- 
tion of  the  test  site,  are  chosen  to  minimize 
disturbance  to  the  surroundings.  The  charge 
parameters  result  in  venting  of  the  bubble  on 
its  first  expansion,  preventing  interaction  of 
its  pulsations  with  the  bottom.  The  test  is 
conducted  with  the  US?  at  a distance  of  ‘"*00 
ft  from  the  nearest  structure,  the  pier  which 
serves  as  the  base  for  the  test  operation  and 
on  which  the  test  control  station  is  located. 
The  geometry  of  the  test  array  is  tightly  con- 
trolled. The  position  of  the  LFSP  is  fixed  by 
an  LCM  tensioning  through  an  outhaul  line 
against  a hold-back  line  free  the  LFS.-  to  the 
pier.  The  charge  location  Is  fixed  relative 
to  the  LFSP  by  a small  boa,  tensioning  through 
a charge-outhaul  line  against  a standoff  bri- 
dle from  the  charge  float  to  the  LFSP.  Follow- 
ing the  test,  the  LFSP  is  winched  back  to  the 
pier  by  a retrieval  line,  with  assistance  from 
the  LCM.  In  addition  to  the  I CM  and  the  small 
boat  which  handles  the  charge-outhaul  line,  a 
sonar  boat  is  employed  to  patrol  the  test  area 
as  a precaution  against  fish-kills,  and  a Yard 
‘lug  stands  by  to  assist  in  case  of  emergency. 

Operational  procedures  are  described  in 
detail  in  Reference  (l). 

CONCLUSION 


A calibration  test  scries  of  fourteen 
shots  was  completed  during  May  1973.  Ko  sub- 
stantial difficulties  were  encountered  in  op- 
erations as  outlined  above,  and  test  results 
indicate  that  the  LFSP  will  not  be  affected 
structurally  by  repeated  tests  at  the  probable 
specification  severities. 

This  work  was  sponsored  by  Kaval  Ship 
Systems  Command  under  Element  6371**,  Project 
SWl5. 

REFERENCE : 

1.  "Operation  Manual  for  Large  Floating  Shock 
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Fig.  3 - Mounting  plate  aetail 


The  Navy  Large  Floating  Shock  Platform 
Part  II:  Shock  Characteristics 


E.W.  Clements 
Naval  Research  Laboratory 
Washington,  D.C. 


Equipment  to  be  installed  on  Navy  ships  must  be  tested  and/or 
designed  to  meet  Navy  shock  specification  standards.  The  previous 
capability  for  direct  testing  extended  to  items  weighing  1+0,000/ 
60,000  lb.  The  Large  Floating  Shock  Platform  has  now  been  con- 
structed and  calibrated,  and  its  shock  characteristics  an  compat- 
ible with  those  of  the  previously  existing  test  devices.  Direct- 
testing  capability  for  shipboard  shock  can  now  be  extended  to 
items  weighing  seme  1+00,000  lb. 


BACKGROUND 

In  order  to  maximize  the  combat  effec- 
tiveness of  its  ships  the  Navy  requires  that 
shipboard  equipment  be  designed  to  resist  the 
damaging  effects  of  shock  caused  by  the  explo- 
sion of  enemy  weapons,  and  be  tested  (if  pos- 
sible) by  exposure  to  shock  similar  to  this 
before  acceptance.  It  is  thus  necessary  to 
have  a means  for  generating  controlled  shock 
environments  similar  to  those  found  on  ship- 
board; this  is  needed  both  for  acceptance 
testing  and  to  provide  a tool  for  design  and 
development.  The  Navy  has  such  a means  in  the 
family  of  machines  specified  in  MIL-S-901  for 
acceptance  testing.  This  family  is  comprised 
of  the  Lightweight  Shock  Machine  (for  loads  up 
to  1+00  lb),  the  Mediumweight  Shock  Machine 
(250-6000  lb)  and  the  Floating  Shock  Platform 
(6000-60,000  lb).  Those  machines  have  little 
in  common  in  design  or  operation,  but  are 
highly  compatible  in  the  shock  environments 
they  provide  to  their  test  loads.  They  have 
now  been  joined  by  the  large  Floating  Shock 
Platform  (LFSP),  which  has  a test  load  capac- 
ity of  up  to  400,000  lb. 

CALIBRATION  TEST  SERIES 

The  LFSP,  described  in  Part  I of  this 
paper,  was  subjected  to  an  eight-shot  cali- 
bration test  series  in  May  1973.  This  was 
preceded  by  a preliminary  series  of  six  shots 
in  March  1973*  shortly  after  the  IFSP  was 
delivered  to  the  West  Coast  Shock  Facility 
(WCSF),  Both  series  were  performed  by  person- 
nel fran  WCSF  and  the  Naval  Research  Labora- 
tory (NRL),  Washington,  D.C. 

The  first  series  was  conducted  against  the 
empty  LFSP;  following  this,  a Floating  Shock 


Platform  (FSP)  was  installed  in  the  LFSP  as  a 
test  load.  Four  shots  of  the  second  series 
nere  fired  against  the  LFSP  with  this  load  con- 
dition, the  inner  bottom  of  the  FSP  was  filled 
with  water  to  increase  the  load,  and  the  final 
four  shots  were  fired. 

The  charge  size  for  all  shot;  was  300  lb, 
and  the  charge  depth  20  feet.  These  para- 
meters were  partially  set  by  geographical  con- 
siderations, e.g. , the  depth  of  ater  available 
and  the  proximity  of  buildings,  etc.  With 
charges  of  this  size  and  depth  the  bubble  vent- 
ed on  the  first  expansion.  The  gross  displace- 
ment of  the  LFSP  is  thus  much  smaller  than  that 
of  the  FSP,  s<nce  it  is  due  to  the  shock  wave 
rather  than  to  riding  the  water  displaced  by 
the  bubble.  The  primary  control  variable  for 
operation  of  the  LFSP  is  the  charge  standoff, 
the  horizontal  distance  between  the  position  of 
the  charge  and  the  closest  point  of  the  LFSP. 
The  secondary  control  variable  is  the  charge 
orientation,  i.e.,  whether  it  is  closest  to  a 
side  or  an  end  of  tne  LFSP.  The  major  effect 
of  standoff  is  to  control  the  severity  of  the 
shock  produced,  that  of  orientation  to  change 
the  relative  severity  of  its  components  in  the 
athwartship  and  fore-and-aft  directions.  The 
standoff  was  varied  from  120  feet  to  45  feet. 
Charges  were  fired  against  both  sides  and  both 
ends  with  no  added  load  in  the  IF3P  and  against 
the  port  side  and  stem  with  added  load.  De- 
tails of  the  shot  geometries  are  given  in 
Table  I. 
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Table  I 


LFSP  Calibration  Test 
Shot  Schedule 


Shot  No. 

Date 

1973 

Near  Side 

Standoff 

ft 

Added  Load 
klb 

1 

oA  x>V, 

Port 

70 

O 

2 

1 Mar 

Port 

55 

0 

3 

2 Mar 

Port 

45 

0 

4 

5 Mar 

Stem 

45 

0 

5 

6 Mar 

Starboard 

45 

0 

6 

7 Mar 

Bow 

45 

0 

7 

7 May 

Port 

120 

112.6 

8 

9 May 

Port 

70 

112.6 

9 

10  May 

Port 

45 

112.6 

10 

16  May 

Stem 

45 

112.6 

11 

17  May 

Port 

70 

176.9 

12 

18  May 

Port 

45 

176.9 

13 

22  May 

Port 

45 

176.9 

14 

23  May 

Stem 

45 

176.9 

TEST  LOAD 

The  FSP  employed  as  a test  load  was  con- 
nected to  the  LFSP  mounting  plane  by  three 
edgewise  plates  running  athwartships,  one  near 
each  end  and  one  amidships.  Shock -protec ted 
racks  within  the  FSP  contained  recorders,  in- 
strm.  ntation  components  and  firing  control 
circuitry.  Also  within  the  FSP  was  a concen- 
trically mounted  single -degree -of- freedom  (SDOF) 
system,  a 5,000-lb  concrete  block  supported  by 
an  elaborate  spring  system.  This  was  designed 
same  years  ago  by  N3HDC  for  another  project, 
and  its  only  important  mode  was  that  of  simple 
vertical  translation  at  a fixed  base  frequency 
of  30  H7.  The  total  weight  added  to  the  LFSP 
was  112.6K  lb,  increased  to  176. 9K  lb  later  by 
filling  the  FSP  inner  bottom  with  water.  The 
mounting  system  was  extremely  stiff;  even  with 
these  large  weights  the  translation  frequencies 
of  the  total  mass  lumped  on  the  total  spring  of 
the  plates  would  be  about  300  Hz  and  240  Hz. 

The  behavior  of  the  sdof  should  thus  be  essen- 
tially the  same  as  if  it  had  been  attached 
directly  to  the  LFSP. 

INSTRUMENTATION 

For  the  first  series  of  tests,  transducers 
were  installed  to  measure  velocity  at  five  lo- 
cations on  the  LFSP  bottom.  These  transducers 
were  acceleraneters  in  a shock-mitigating  hous- 
ing with  integrating  electronics.  For  the  sec- 
ond series  of  tests  more  of  these  transducers 
were  added  on  the  LFSP  mounting  plane  at  the 
base  of  the  FSP  mounting  plates  and  within  the 
FSP.  Strain  gages  were  installed  on  the  LFSP 
shell  plating  at  its  mid-section  and  on  its 
central  athwartship  stiffener.  The  output  sig- 
nals were  recorded  on  magnetic  tape.  After 
each  shot  the  tapes  were  played  back  through  an 
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integrator  on  a string  oscillograph  providing 
a strip  recording  of  velocity  and  displacement 
vs.  time  for  each  velocity  transducer,  and 
strain  versus  time  for  each  strain  gage.  A few 
of  the  velocity  records  were  also  fed  into  an 
analog  shock  spectrum  analyzer. 

LFSP  VELOCITIES 

The  waveforms  of  velocities  measured  in 
the  LFSP  are  affected  principally  by  the  lo- 
cation and  orientation  of  the  measurement.  The 
magnitude  of  the  velocity  (in  a given  direc- 
tion) decreases  as  the  distance  between  the 
point  of  measurement  and  the  charge  increases. 
This  gives  a general  increase  in  velocity  as 
the  charge  standoff  is  reduced,  and  also  a con- 
siderable variation  in  the  velocities  measured 
throughout  the  USP  for  a given  shot.  The 
shapes  of  the  waveforms  are  mainly  determined 
by  the  orientation  of  the  measurement.  Those 
in  the  vertical  and  horizontal-parallel"  direc- 
tions have  sharp  rises  1 ms)  and  gradual 
declines  50  ®s)  accompanied  by  structure- 
related  frequencies  up  to  ~ 1 KHz.  The  hori- 
zontal-transverset  velocities  consist  of  the 
structural  frequencies  amplitude-modulated  by  a 


* The  horizontal -parallel  direction  is  the  hor- 
izontal direction  parallel  to  the  line  from  the 
plar  position  of  the  charge  to  the  nearest 
point  of  the  USP,  i.e. , athwartship  for  a 
side  shot,  fore-and-aft  for  an  end  shot. 

t The  horizontal-transverse  direction  is  the 
horizontal  direction  perpendicular  to  the  line 
fran  the  plan  position  of  the  charge  to  the 
nearest  point  of  the  LFSP. 
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fairly  fact  rice  10  mr.)  nn'l  slow  decline 
200  ns) (Fi  -.  1).  Velocities  measured  on  the 
shell  plating  an<l  the  IFTJ?  noun  tine  plane  on 
top  of  the  stiffeners  are  similar,  although 
those  on  the  plating  have  higher  magnitudes  and 
sharper  rises,  so  sharp  that  it  is  likely  that 
a substantial  fra^pent  of  the  energy  content 
lies  at  frequencies  above  the  2 KHz  pass  limit 
of  the  velocity  transducer.  Thus  the  measure- 
ments made  on  the  plating  must  be  considered 
less  reliable  than  those  on  the  mounting  plane 
or  FSP.  Regardin'’-  the  mounting-plane  veloci- 
ties, the  variation  of  peak  vertical  velocity 
with  measurement  location  is  such  that  the 
greatest  is  1.75  tines  the  smallest.  This 
variation  is  enough  to  obscure  see/-  of  the 
other  points  of  interest,  and  It  is  convenient 
to  speak  of  the  average  value  of  all  measured 
values  in  a given  direction  and  at  a given  type 
of  location — viz.,  the  average  vertical  mount- 
ing-plane velocity  or  the  average  ath'.ra.rt::hip 
FSP  velocity,  etc. 


Fig.  1 - Typical  velocity  waveforms  measured  on 
the  LFSP  mounting  plane. 


The  variation  is  less  for  horizontal-parallel 
velocities,  and  slight  for  horizontal-trans- 
verse. On  the  average,  the  variation  in  peak 
velocity  with  measurement  orientation  is 
1:.7:.3  for  vertical:  horizontal-parallel: 
horizontal-transverse.  The  essential  effect 
of  charge  orientation  (end  shot  vs  side)  is  to 
determine  which  of  the  LFSP  axes  (athwartship 
and  fore-and-aft)  is  aligned  with  which  of  the 
test  axes  (horizontal -parallel  and  horizontal- 
transverse).  It  would  be  expected  that  the 
average  peak  velocities  would  be  slightly  lower 
for  end  shots  than  for  side,  since  the  average 
point  of  the  LFSP  is  more  removed  freo  the 
change  In  this  orientation.  There  is  such  a 
decrease,  but  it  is  not  consistent  and  is  of  an 
amount  which  could  be  expected  frem  normal  vari- 
ation between  nominally  identical  shots.  The 


peak  velocities  on  the  plating  shows  r.o  effect 
from  an  increase  in  load  from  0 to  112. CK  lb, 
and  a substantial  decrease  for  an  increase 
frem  112.6  to  176. y.  lb.  The  mounting  plane 
peak  velocities  show  a similar  decrease  for  an 
increase  in  load  weight  from  112.6  to  176. UK  lb. 
These  effects  are  s aamarized  in  Fig.  2A  fc  2B, 
showing  average  peak  velocities  vs  standoff. 


LOAD  (FSP)  VFI/XITEC 

The  velocity  waveforms  measured  in  the  >GP 
are  similar  to  those  found  at  the  same  loca- 
tions when  the  FSP  is  operated  alone  (l),  im- 
plying that  because  the  mounting  structure  is 
stiffer  than  the  inherent  structure  of  the  FSP 
(fret-free  beam  node  ~ 120  Hz)  the  details  of 
the  local  FSP  structure  are  dominant  in  fixing 
the  details  of  a particular  measured  waveform. 
Most  of  the  transducers  in  the  FSP  were  placed 
around  its  perimeter  close  to  the  mounting 
plates,  the  regions  most  rigidly  coupled  to  the 
LFSP  mounting  plane.  The  waveforms  here  resem- 
ble those  of  the  associated  LFSP  waveforms  with 
the  addition  of  a discernible  sinusoidal  compo- 
nent. Those  measured  towards  the  FSP  center 
are  located  over  cells,  and  have  a more  strong- 
ly sinusoidal  nature.  The  magnitudes  and  the 
variations  of  magnitude  are  greater  in  t«£  FSP 
than  at  the  LFSP  mounting  plane.  The  highest 
peak  velocities  are  those  from  the  central 
locations  above  cells,  but  those  about  the 
perimeter  are  also  higher  and  more  variable 
from  place  to  place  than  those  in  the  LFSP. 

In  the  vertical  direction,  the  highest  peak 
velocity  is  2.5  times  the  smallest,  and  the 
average  peak  is  1.5  times  the  average  peak  at 
the  LFSP  mounting  plane.  In  the  two  horizontal 
directions  the  variation  over  location  is  much 
leas,  and  the  average  peak  values  are  substan- 
tially the  same  in  the  FSP  and  the  LFSP.  The 
ratio  of  average  peak  velocities  in  the  FSP  is 
1:. 1*5: -15  (vertical :horizontal-parallel hori- 
zontal-transverse). The  Influences  of  the  test 
parameters  are  much  the  same  as  on  the  LFSP  ve- 
locities, save  that  the  vertical  FSP  velocity 
is  more  strongly  reduced  by  increasing  load 
weight  (Fig.  3). 


LOAD  (SDOF ) RESPONSES 

The  SDOF  response  was  measured  by  a sin- 
gle accelerometer  sensing  vertical  notion 
centrally  located  on  top  of  its  mass.  The 
signal  from  this  transducer,  its  integral  and 
double  integral  were  played  hack  on  strip 
charts.  The  acceleration  curve  carries  a small 
amount  of  high-frequency  hash  for  the  first 
couple  of  cycles,  but  after  this  dies  down  the 
curve,  like  those  of  velocity  and  displacement, 
is  a clean,  highly-damped  sinusoidal  at  a fre- 
quency of  29.1*  Hz.  The  measured  peak  values 
are  tabulated  in  Table  II. 
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Fig.  2 - The  variation  of  average  peak  velocity 
with  charge  standoff  - A,  ui  the  shell -plating; 
B,  on  the  WSP  mounting  pl'ne. 


Fig.  3 - Variation  of  average  peak  load  (FSP) 
velocities  as  a function  of  charfe  standoff. 
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Fig.  4 - A typical  shock  spectrum  derived  free 
a velocity  recording. 


Fig.  5 - Design  shock  spectra  for  the  LFSP-A, 
the  variation  of  shock  spectrin  value  with 
modal  frequency  and  conpone.  t direction;  B,  the 
variation  with  modal  weight. 


LFSP  SHOCK  SPECTRA 

Shock  spectra  were  derived  from  the  re- 
cordings of  seme  of  the  vertical  velocities  at 
the  USP  mounting  plane.  These  generally  show 
a velocity-shock  nature  below  about  Uo  Hz,  and 
a few  interaction  peaks  above,  as  in  Fig. a . 

The  dips  at  around  30  and  ho  Hz  are  present  In 
sene  degree  in  all  of  these  spectra  and  in  the 
absence  of  any  knowledge  of  the  nodal  descrip- 
tion of  the  test  load  structure  the  average 
level  of  the  naxinax  spec t mi  in  this  regicr  is 
taken  as  the  shock  spectrin  value.  The  dis- 


placement and  acceleration  limits  can  be  8:t 
from  the  integrals  and  slopes  of  the  velocity 
curves,  since  for  a set  of  massless  oscillators, 
such  as  define  the  shock  spectrum,  the  peak 
relative  displacement  for  low  frequencies  is 
the  maximum  displacement  of  the  shock  motion, 
and  the  peak  acceleration  at  high  frequencies 
is  the  maximum  acceleration  involved  in  the 
shock  motion.  The  data  at  hand  do  not  indicate 
a change  in  corner  frequency  with  weight. 
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Table  II 
SDOF  Responses 


Shot  Ho. 

Hoar  Side 

Standoff 

ft 

Added  load 
klb 

Peak  SixjF 
Ac cole r. 

S 

Response 
Velocity 
ft/ sec 

Displ. 

in. 

7 

Port 

120 

112.6 

45 

5.6 

.49 

8 

Port 

70 

112.6 

63 

8.6 

1.35 

9 

Port 

45 

112.6 

38 

13.3 

2.40 

10 

Stern 

U* 

112.6 

106 

11.7 

1.96 

11 

Port 

70 

176.9 

56 

8.9 

1.47 

13 

Port 

45 

176.9 

91 

14.2 

2.07 

14 

Stem 

45 

176.9 

87 

12.0 

2.22 

Oata  for  larger  loads  may  do  so.  Next  the 
variation  of  the  shock  spectra"  value  with 
modal  weight  must  be  determined.  This  has 
been  attempted  on  the  basis  of  the  assumptions 
that  the  shock  spectrum  value  components  are 
related  much  as  peak  velocity  components,  ar.d 
that  the  weight  in  the  lowest  frequency  mode 
is  some  6o£  of  the  total.  The  first  is  rea- 
sonable for  linear  systems,  the  second  a rule 
of  thumb  for  simpl"  systems  in  uniaxial  trans- 
lation. Under  these  assumptions  the  shock 
spectra  curves  of  Fig.  5A  & 5B  are  obtained. 


Fig.  6 - The  effect  of  charge  standoff  and 
load  on  shock  spectrum  .alue 

Figure  5A  shows  the  variation  of  shock  spectrum 
value  with  modal  frequency  for  each  component 
direction.  The  corner  frequencies  are: 
Vertical,  9.5  and  117  Hz;  horizontal-parallel, 
31  and  154  Hz;  horizontal-transverse,  59  and 
171  Hz.  The  shock  spectrum  values  are: 
vertical,  V ; horizontal-parallel,  .7  V&; 


The  variation  of  V with  modal  weight  is  shout 
ir.  Fig.5B  for  45  f?  standoff,  and  is  124  ips 
for  modal  weights  up  to  about  90E  lb,  declining 
to  ICC  ips  at  142K  lb.  This  curve  should  con- 
tinue to  decline,  and  eventually  level  off  at 
sene  value;  its  structure  will  become  clearer 
as  data  are  accumulated  throu^i  use.  A 
straight  line  connecting  the  points  for  the  two 
load  conditions  on  a linear  plot  has  the  equa- 
tion 

VQ(ips)  = 156  - .38  X modal  weight  (K  lb). 

The  variation  of  V with  standoff  is  shown  in 
Fig.  6. 

STRUCTURAL  EFFECTS 

During  the  first  six  shots  (LFSP  unloaded) 
minor  cracking  of  the  welds  associated  with 
some  of  the  side  stiffeners  occurred.  These 
were  repaired,  and  no  cracking  was  observed  on 
the  subsequent  shots.  The  strain  gages  showed 
elastic  behavior  of  the  LFSP  structure  with  the 
exception  of  the  installation  on  the  bottom 
plating  at  the  center  of  the  cell  closest  to 
the  charge.  For  a shot  at  45  ft  standoff,  this 
showed  a set  of  100  to  200  uin/in.  For  a given 
shot  geometry,  the  peak  strains  were  higher 
with  the  greater  load.  However,  the  differences 
between  the  peak  strains  for  the  two  nominally 
identical  shots  were  about  twice  as  great  as 
the  difference  between  the  lesser  of  these  and 
the  shot  with  the  lower  load  and  the  same  standr 
off.  Since  it  is  not  likely  that  a standoff  of 
less  than  45  ft  will  be  required,  the  LFSP 
structure  should  not  be  significantly  affected 
by  normal  use. 

SUMMARY 

The  Large  Floating  Shock  Platform  has  ex- 
tended the  Navy's  shuck  testing  capability  to 
loads  of  about  400,000  lb.  In  its  essential 
shock  cliaracteristics  the  LFSP  behaves  in  much 
the  same  way  as  the  FSP,  although  some  details 
(for  example,  gross  motion)  are  different. 

The  LFSP  will  be  included  as  a prescribed 
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testing  device  in  future  editions  of  M.IL-3-901. 
Under  -what  will  probably  be  its  r.omal  operat- 
ing conditions  the  LFS?  will  be  essentially  an 
elastic  machine,  so  that  its  structure  and 
shock  characteristics  should  change  little 
with  use.  It  would  oe  desirable  to  obtain 
information  on  the  LrSP’s  behavior  with  ex- 
tremely large  loads  in  order  to  extend  the 
curve  of  shock  spectrum  'value  vs  nodal  weight. 
Thin  information  should  becct..  available  as 
large  items  are  tested. 
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DISCUSSION 

Mr.  Mador  (Electric  Boat) : Do  you  Intend  to 

use  a similar  size  charge  in  conducting  MIL- 
S-90IC  teats  for  this  large  barge  as  you  do  on 
the  smaller  barges? 

Mr.  Schrader:  This  is  a 300  lb.  charge. 

Mr.  Mador:  Yes,  I realize  that  you  use  that 

for  the  instrumentation  shock  characteristics. 
But  if  you  were  to  use  this  platform  for 
qualification  of  large  components  what  would 
your  charge  geometry  be? 

Mr.  Schrader:  The  300  lbs.  charge  would  be 

used  with  a 20  ft.  depth.  The  closest  in 
standoff  would  be  50  ft.  and  the  next  would 
be  65,  80,  and  then  1 inch  at  100  ft. 

Mr.  Berezow  (Naval  Ordnance  Laboratory) : You 

mentioned  that  the  buble  was  vented  from  the 
primary  pulse,  was  that  because  of  the  300  lb. 
charge? 

Mr.  Schrader:  Yes,  as  you  probably  know  about 

J5Z  of  the  energy  of  the  explosive  that  we 
used  is  in  the  bubble.  If  we  were  closer  to 
the  bottom  this  bubble  would  then  migrate  to 
the  bottom  and  bounce  on  the  bottom  several 
times  before  it  would  vent,  and  so  we  used 
the  shallow  20  ft.  standoff  to  let  the  bubble 
blow  to  the  surface  to  give  us  some  airborne 
noise  in  lieu  of  exciting  all  of  the  residents 
of  the  area  who  might  think  that  an  earthquake 


was  present. 

Mr.  Roberts  (Raytheon  Company) : Did  you  feel 

that  you  could  have  had  an  oscillating  bubble 
if  you  had  a small  charge?  Did  that  affect  it? 

Mr.  Schrader:  Yes  we  could,  but  that  would 

have  meant  that  we  wouldn't  have  as  good  a 
velocity  distribution  throughout  the  whole 
barge;  and  as  it  v ■>.  the  minimum  to  the 
maximum  velocity  < - ibution  ranged  about  1.75. 

Mr.  Roberts:  In  previous  reports  I noticed 

that  there  were  oscillating  bubbles  with  smaller 
charges.  I was  trying  to  figure  out  why  in 
some  cases  they  were  oscillating  and  in  other 
cases  they  did  not. 

Mr.  Schrader:  With  the  small  floating  shock 

platform,  if  the  water  is  deep  enough,  the 
bubble,  as  you  previously  indicated,  will 
pulsate  about  three  times  before  it  verts. 

This  usually  doesn't  bother  anything  because 
most  ’•rings  are  broken  in  that  first  utili- 
se* ,«d  rise  and  the  only  time  this  bubble  gets 
t > be  annoying  is  if  there  is  a low  frequency 
s> item. 

Mr.  White  (Bettis  Atomic  Power) : Do  you  have 

plans  to  build  an  additional  large  barge?  You 
had  implied  that  this  was  the  only  one  in 
existence.  I was  wondering  if  you  had  plans 
to  build  one  for  the  East  Coast. 

Mr.  Schrader;  No,  not  that  I know  of.  Thesi- 
larger  tests  are  rather  expensive  and  I don’t 
think  that  they  are  planned  as  regular  pro- 
duction tests.  These  tests  that  will  provide 
information  for  prototype  equipment,  such  as 
a large  reduction  gear,  not  every  piece  of 
machinery  will  be  tested  this  way  maybe 
Mr.  Sullivan  has  some  ideas  on  this. 

Mr.  Sullivan  (NAVSEC):  The  requirements  for 

this  barge  have  already  been  written  into  a 
revised  M1L-S-901  that  is  presently  going 
through  the  mill,  so  this  will  be  Included. 

Again  since  it  is  a unique  facility  I don't 
anticipate  that  we  will  automatically  divert 
all  large  pieces  of  equipment  through  this, 
but  it  will  be  alternative  to  a formal  dynamic 
analysis.  We  would  hope  that  more  equipment 
is  diverted  for  tests  rather  than  strictly 
analysis. 

Mr.  Duffy  (General  Electric  Company) : Do  you 

have  any  data  on  relative  costs  between  the 
small  platform  and  the  large  platform? 

Mr . Schrader : The  costs  of  this  test  seemed  to 

run  about  20  or  30Z  mote  than  with  the  small 
platform,  and  this  is  just  conducting  the  tests. 
The  installation  of  the  bigger  equipment  can  be 
pretty  expensive,  it  depends  on  what  kind  of 
equipment.  This  foundation  cost  us  about 
$20,000  to  be  installed. 

Mr.  Welch:  Is  this  large  enough  to  shock  test 

the  largest  item  of  equipment  aboard  a Naval 
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combatent  ship?  I can  answer  for  machinery, 
turbines,  gears,  and  turbine  generator 
sets  and  say  that  it  is.  I have  tone  question 
about  boilers  and  perhaps  reactor  coaponents. 

I think  that  this  large  floating  shock  platform 
rounds  out  the  Naval  shock  testing  and  shock 
analysis  capacity  in  a very  fine  and  complete 
fashion.  It  has  been  the  experience  of  those 
who  have  used  dynamic  analysis  that  it  is 
very  effective  in  the  early  stages  of  design. 

It  is  a design  method  in  that  it  permits  one 
co  make  changes  during  design.  But  because  of 
the  necessarily  mathematical  complexities  it 
tends  to  be  rather  gross  in  the  math  models 
that  one  develops  for  the  dynamic  analysis,  or 
rather  gross  in  their  concept  of  the  dynamical 
system.  Furthermore  it  generally  does  not  take 
into  account  the  cross  effects,  the  cross 
couplings  that  exist  between  the  directions  of 
shock;  on  the  other  hand  the  shock  barge  test 
tends  to  be  very  detailed  and  an  apparatus 
subjected  to  the  shock  barge  test  generally 
gets  an  extensive  shock  in  all  directions,  and 
although  there  may  be  some  criticism  '.bout  the 
arbitrary  nature  of  the  shock  barge  test,  it 
at  least  encompasses  detailed  elements  of  an 
apparatus  that  are  often  neglected  when  one  math 
models  the  apparatus  for  dynamic  analysis. 
Furthermore  it  takes  into  effect  sometimes  very 
valid  cross  effects  which  are  important  to  the 
actual  shock  resistance  of  machinery. 

Mr.  Schrader;  All  of  the  test  loads  that  have 
been  considered  for  tests  on  the  barge  looked  as 
if  they  would  be  accomodated.  There  might  be  a 
stability  problem  with  some  of  the  hardware  t'.at 
has  a very  high  center  of  gravity,  but  off  hand 
I don't  know  of  any  machinery  that  won't  fit 
into  this. 

Mr.  Sullivan;  Originally  we  attempted  to  come 
up  with  a size  that  would  satisfy  our  projected 
requirements.  It  will  take  all  the  boilers 
that  the  Navy  has  made,  including  the  LHA 
boiler  which  is  the  biggest  boiler  ever  made 
for  shipboard  service.  There  may  be  some  other 
components  in  the  reactor  area  that  we  cannot 
handle . 

Mr.  Schrader:  I agree  with  Mr.  Welch  whole- 

heartedly that  these  tests  are  very  ingenious 
in  finding  out  one  weak  point  that  somebody  has 
not  calculated  in  their  dynamic  design  and  that 
6th  mode  which  was  not  considered  to  be  import- 
ant. These  tests  have  been  very  surprising  in 
the  way  that  they  will  point  out  things  to  you. 
As  an  example,  in  the  calibration  series  the 
welding  rod  that  was  used  to  fasten  this  founda- 
tion into  the  FSP  it  started  to  crack  at  the 
joint  and  we  couldn't  understand  it  so  we  took  a 
careful  look  at  the  rods  that  were  used  and  low 
and  behold  it  was  an  Incorrect  rod.  This  had 
started  the  crack  and  you  may  have  noticed  that 
we  put  some  little  gussets  to  spread  out  the 
load  at  the  cracking  point  and  we  had  to  rewe’d 
things;  but  it  is  rather  interesting  how  shock 
in  a test  like  this  can  pick  out  that  one  flaw 
in  a whole  system  just  like  thst. 
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THE  EFFECT  OF  UNIFORM  EXTERNAL  PRESSURIZATION  ON  THE 
DYNAMIC  RESPONSE  OF  ISOTROPIC  CYLINDRICAL  SHELLS 


Frederick  J.  Dzialo 
Civil  Engineering  Department 
University  of  Massachusetts 
Amherst.  Massachusetts 


Following  Flugge's  derivation  for  the  buckling  of  cylindrical  shells, 
the  equations  of  motion  for  dynamic  loading  of  isotropic  cylindrical 
shells  subjected  to  an  external  uniform  pressure  have  been  formulated. 

The  equations  of  motion  neglect  the  effect  of  rotatory  inertia  and 
shear,  and  in  general  are  applicable  to  cylindrical  shells  whose 
structural  wave  length  is  significantly  greater  than  the  thickness 
of  the  shell;  however  as  pointed  out  by  ’.E.  Greenspon  [3],  Flugge's 
theory  which  includes  bending  and  membrane  effects  is  good  at  even 
the  shorter  wave  lengths,  provided  the  thickness  of  the  shell  is  not 
too  great,  that  is  the  ratio  of  inside  radius  to  outside  is  less 
than  0.5.  The  normal  mode  theory  was  utilized  to  provide  dynamic 
solutions  for  the  equations  of  motion. 

Solutions  are  also  provided  for  the  Timoshenko- type  theory,  and 
comparisons  are  made  between  the  two  theories  by  considering  and 
neglecting  in-plane  inertia  forces. 

Comparison  of  results  is  exemplified  by  a numerical  exanple  which 
considers  the  effect  of  an  external  pressure  on  the  dynamic  response 
of  a shell  supported  by  a thin  diaphragm  and  subjected  to  a localized 
unit  impulse. 


INTRODUCTION 

In  a recent  paper  [1],  E.N.K.  Liao  and  P. 
G.  Kessel  investigated  the  dynamic  response  of 
a cylindrical  shell  with  Initial  stress,  and 
subjected  to  general  three-dimensional  surface 
loads.  This  paper  presented  general  solutions 
for  both  Flugge's  and  Donnell's  equations  for 
a simply  supported  cylindrical  shell,  as  well 
as  analytical  solutions  for  a stationary  radial 
point  force  and  a stationary  couple. 

The  purpose  of  this  paper  is  to  investi- 
gate the  effect  of  an  external  uniform 
pressure  on  the  dynamic  response  of  isotropic 
cylindrical  shells.  Following  Flugge's  [ZJ 
derivation  for  the  buckling  of  cylindrical 
shells,  the  equations  of  motion  are  formulated. 
The  equations  of  motion  neglect  the  effect  of 
rotatory  inertia  and  shear,  and  in  general  are 
applicable  to  shells  whose  structural  wave 
length  is  significantly  greater  than  the  thick- 
ness of  the  shell;  however  as  pointed  out  by 
J.E.  Greenspon  [3],  Flugge's  theory  which 
Includes  bending  and  membrane  effects  is  good 
at  even  the  shorter  wave  lengths  provided  the 
shell  is  not  too  thick;  that  is  the  ratio  of 
inside  radius  to  outside  radius  is  less  than 


0.5.  The  equations,  useful  in  calculating 
displacements  and  stresses  for  impulsive 
loads  applied  to  a comparatively  small  region 
of  the  shell,  are  found  by  utilizing  the 
classical  normal  theory  which  provides  dynamic 
solutions  for  the  equations  of  motion  as  well 
as  natural  boundary  conditions  for  a cylindrical 
shell  subjected  to  an  external  uniform  pressure. 

Solutions  are  also  provided  for  the 
Timoshenko-type  theory,  and  comparisons  are 
made  between  the  two  theories  by  considering 
and  neglecting  in-plane  i <ertia  forces. 

Comparison  of  results  is  exemplified  by  a 
numerical  example  which  considers  die  effect  of 
an  external  uniform  pressure  on  the  dynamic 
response  of  a shell  simply  supported  by  a thin 
diaphragm,  and  subjected  to  a localized  unit 
impulse. 
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Equations  of  Motion 

Following  Flugge's  [2]  exact  derivation  for 
the  buckling  of  cylindrical  shells,  the  differ- 
ential equations  of  motion  for  impulsive  load- 
ing of  cylindrical  shells  subject  to  an  external 
uniform  pressure  become: 
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Substitution  of  equation^  (4)  through  (13)  into 
equations  (1)  through  (3)  yields: 
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Orthogonality  and  Modal  Vibrations 

For  free  vibrations,  equations  (14)  through 
(16)  may  be  written  as  follows 
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} = i-p  (1  ♦ 3k)  - q2 
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Equations  (17)  through  (19)  yield  the 
natural  frequencies  and  mode  shapes.  The  ortho- 
gonality condition  for  the  mode  shapes  is 
derived  by  assuming  that  the  displacements  u, 
v,  and  w for  the  nth  mode  have  the  form 
1w  t 1(1)  t 
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Finding  the  orthogonality  condition  Involves 
the  following  steps: 

1.  the  nth  terms  of  expressions  (20)  are  Insert- 
ed Into  equations  (17)  through  (19),  and  the 
resulting  equations  are  multiplied  by  u^x), 
v_(x)  and  wm(x),  respectively.  Integrated  over 
the  domain,  and  added, 

2.  the  irth  terms  of  expressions  (20)  are  In- 
serted Into  equations  (17)  through  (19),  and  the 
resulting  equations  are  multiplied  by  u.(x), 

v.  (x) , and  wn(x),  respectively,  integrated  over 
tne  domain,  and  added; 

3.  the  equation  resulting  from  Step  2 Is  sub- 
tracted from  the  equation  resulting  from  Step 

1;  It  Is  Integrated  by  parts,  and  use  Is  made  of 
equations  (4)  through  (13)  to  obtain  the  final 
orthogonality  relation.  The  orthogonality 
condition  may  be  written  as  follows: 
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and  t are: 

Fixed: 

Hi  nged: 

u = v = w = 0 Mx  * 0 
Simply  Supported: 


N - -5i  - p « 0 N + pw  - P » 0 
xa  a 3x  x ^ 3x 


(25) 

Integral  of  the  Square  of  Eigenfunctions 

The  Integral  of  the  square  of  the  eigen- 
functions Is  evaluated  from  equation  (21)  by 
a limiting  process.  For  any  prescribed  boundar 
condition,  the  evaluation  of  the  Integral  may 
be  determined  as  follows: 
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To  determine  the  natural  frequencies  and  mode 
shapes,  equations  (17)  through  (19)  may  be 


w-n 

solved  by  assuming 

, x-u 

u=Ae''x^acos(m$)ei“nmt 

(22) 

v*Be 1 x^a s 1 n (m* )e * “nm* 
w=Ce>'x^acos(m$)e1“Jnmt 

(27) 

(23) 

Inserting  equation  (27)  Into  equations  (17) 
through  (19)  yields  equation  (28): 
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Free: 
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The  characteristic  equation  is  found  by 
setting  the  determinant  of  the  set  of  equations 
(28)  equal  to  zero.  To  determine  the  eigen- 
values, the  following  method  is  utilized: 

A value  of  guessed  and  Inserted  into  the 
characteristic  equation.  The  characteristic 
equation  will  yield  eight  roots.  For  unequal 
roots,  equations  (27)  may  be  written  as 
follows: 
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where  for  each  x^  there  exists  a relationship 
between  the  amplitudes  A;,  Bj  and  Cj  from  the 
determinant  of  equation  (28). 


Equations  (29)  with  the  necessary  boundary 
conditions  will  lead  to  a determinant  | j | = 0 
versus  wj^,.  The  eigenvalues,  , are  those 
for  which  | a - j | - 0.  At  a point,  wjL,  when 
|ay | = 0,  the  ratio  of  the  amplitudes  Ay  B^ 
andJC^  can  be  calculated  from  the  determinant 
of  eqution  (28). 


Substituting  the  above  equations  into  equations 
(14)  through  (16),  and  utilizing  the  ortho- 
gonality condition  (21)  yields  the  following: 
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where 
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For  an  impulsive  loading  as  shown  in  Figure  1, 
equation  (31)  becomes 


%{t) 


_ [(rc])/afc-c|fotf>si,wn(t-i)dXdxd» 

TJo'i'^tu^v^Jdxdt 


(32) 


For  a concentrated  impulsive  loading  equation 
(31)  becomes 

...  f’+tJ/a  r«. 

qn(t)  = lim  r 1 

ei"°  Ve,)/a 


'2  P 


sinun(t-i)dxdx  d* 

X wnfZ7ip"h(u"*VH,dxd* 

jo  Jo 


(33) 


where 


Px  ' 4e,e 
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Solution  for  Impulse 

Consider  an  impulse  per  unit  area,  ix(x,ip) , 
i (x,  ) and  ir(x,  ) acting  on  the  cylinder  for 
an  infinitely  short  time.  The  cylinder  may  now 
be  considered  to  be  vibrating  freely  with  the 
following  initial  conditions: 

At  t = 0: 


u * v - w = 0 


(34) 


3U 

at 


ix(x,  t) 

pH 


av  _ y«»  ♦> 
at  pTi 


aw 

at  pH 


(35) 


Solutions  for  Forced  Vibrations 


Equations  (17)  through  (19)  may  be  solved  by 
assuming 

U * jo  Un(X’  ° q"(t)  V = Jo  Vn(x*  *]  qn(t) 

oo 

W ' n=0  Wn<X’  ^ q,,(t)  (30) 


The  displacements  for  free  vibrations  are 
given  as: 


/ 
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where 


U * Jo  Un(A"  C0S  ^ + B,b  Sl"  “"t} 

V “ Jo  V"(A*  C0$  “"l  + B|" s1n  “-0 

w * Jo  cos  V + B.«  s1n  V}  <36) 


Substituting  the  initial  conditions  into 
equation  (36)  and  making  use  of  orthogonality 
yields  the  following 

■ju  <z> 

u * l u B sin  rt  v * l v B sin  u„t 


w » l w_B_  sin  j_t 
«K>  B"n  ■ 


(37) 


where 


f2»  ri 

(1 vu_  + 1 v + i w ) dx  d* 
J a J a x hi  9 m r in 


B 

ID  ui  tC  If  (1 


n )V  L ph(u* + v* + ^ dx  dt 


(38) 


For  a distributed  impulse  as  shown  in 
Figure  1,  equation  (38)  becomes 


/(n+ei)/a  k+c2 

[1  *um  + 1 . v ♦ ^ -W_]  dx  d* 
J(n-c,)/a  ic-c,  X m ♦ M r ■ 


“m 


i;  f.  °h(um  + vm  ♦ dx  d* 


(39) 


For  a concentrated  impulse,  equation  (38) 
becomes 


lim  f(n+ei)/a  (•c+e2 

^ [i„Um  + Vm  + Vj  dx  d* 

ci"KjI/  I,  xid  9 m rin 

J(n-ei)/a  JC-f.2 


an  \J  f*  oh<um  + vm  + wm>  dx 


(4D) 


Equations  of  Motion  for  Timoshenko  Theory  f4] 

Equations  (1)  through  (3)  can  be  reduced 
to  those  presented  by  Timoshenko  by  assuming 
the  following  conditions: 


a.  The  circumferential  strain  e^,  e and 
y»*  are  equal  to  zero  in  calculation 
o“x#  and  Xx,. 

b.  Membrane  forces  are  not  affected  by 
heading  stresses,  nor  bending 
moments  by  membrane  stresses. 


Assumptions  (a)  and  (b)  yield  N.x  - Nx.  and 
M$x  = Mx..  Assuming  e*  = (v'  + w)/a;  u', 
u"  * 0,  and  y * * (u'  * v')/a  = 0,  equations 
(1)  through  (3)  become: 


"i  ♦ aN’x  + pa(v"  + w')  + a2px  - pha2  0 

aV  aNi*  - a V Pv' + 8 V pha 20 

-aQ*  - aQ'  + aN  - pa(w"  ♦ w)  - Pw"  - a2p 
v x 9 r 

* pha2  0 (41-43) 

The  membrane  forces  and  moments  from  equations 
(6)  through  (13)  become 


•N  * Vv*  + w + vu') 

9 a 

H.k  * J (^)(u*  ♦ v’) 

9X  a C 

M = (“V*  + w”  + vw") 

9 a 

Mtx  = fr(1  - v)(w”  - »■) 


(44) 


N = ^ (u‘  + W + vw) 

X a 

IL.  * 7 ITU''  + V’> 

X9  a c 

M “ T-  (w"  + VW"  - w') 

X 3 

\ - b (1  ' v)(w"  " y,) 
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Substitution  of  equations  (44)  into 
equations  (41)  through  (43)  yields  the 
following: 

u"  ♦ ‘ ♦ i^u-  ♦ uw*  + q^v'*  ♦ w') 

. fl2px  _ oha2  32u 

+ _ir  * TTIF 

(!_£_«)„••  ♦ v ♦ (-!-£-%"  ♦ w‘  - k[w"‘  ♦ w— ] 

♦ k[(l  - v)V  ♦ V]  - q2v"  ♦ ^ f£ 

vu*  ♦ v*  ♦ w ♦ k[w"*+  2w"''  ♦ tr~]  - k[v“‘ 


+ (2  - v)v"‘]  ♦ q2w"  ♦ q^(w  * w) 


*2pr  _ -oha2  32w 

“H5 T“3F 


(45-47) 


Equations  (45)  through  (47)  may  be  written 


+ a 


2 px 

r 


(i-j-^)u'*  ♦ b3v-  ♦ e4v"  ♦ w-  - k(w-  ♦ W") 


vu‘  ♦ v*  ♦ 8gW  ♦ k(w“  + 2w""  ♦ W •) 

- k[v‘‘*  ♦ (2  - v)vM<]  ♦ q2w“  ♦ qf" 


fl2pr  _ -oha2  s2* 

- nr  “T-  w 


(48-50) 


where 


u “ H U«»i  cos  "♦  cos  i 

ib  n l 

v * T 7 V sin  in*  sin 

■ n m * 

w * I [ cos  m*  sin  ^ 


To  determine  the  natural  frequencies  and 
mode  shapes  the  determinant  for  the  frequency 
equation  becomes: 


.[«3(!5i)-i,U|5V(a;»),i,(aji)3]wm . o 

^(•>0 V»  - tV-il2?’2 

. Pha2^ 

D rw,J  mn 


-[a3(^y-  ^V)  (E*  )*kp£h  Um 
♦[«4«*k(^)(*)(5^)2]»m  = 0 

♦[l-(2k*l-a4)m2-(Ei)2(l-a,.) 

♦ (52) 


6,  - ♦ q1  B2  = V ♦ q, 

B3  ■ k ♦ 1 B4  * (1  - v)(k  + 0.5)  - q2 
65  ^ 1 + ql 

Illustrative  Example  for  Cylinder  Supported  by 
Thin  Diaphragm 

Flugge's  Theory 

For  a cylinder  supported  by  a thin  dia- 
phragm the  following  displacements  satisfy  the 
natural  boundary  conditions  as  derived  from  the 
orthogonality  conditions: 


The  ratios  of  the  mode  shapes  are  defined 


as: 


_ Umn  . -CD  ♦ BE 

"mn  nr ' mmF 

mn 


mn  _ -AE  ♦ BC 
'mn  BT  ‘ AD  - 


(53) 


where 


mn 


ha2  2 
'mn 


* _ „ #n*a»Z  2 , pha^  ; 

A *a5'T”'  ~ m a2  * ^ 

B 


4^m(^) 
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— 


c = a3(nr)  - + M~>3 

n - „ -2  „ /tua>2  . pha2  2 

D V -“lhr > +T“^ 


E = -a4m  - k(^-^)m(£-p.)2 


Solutions  for  Forced  Vibrations 

From  eauatlons  (30)  and  (31)  the  dynamic 
displacements  become 


u ' I I Umn  cos  m^cos 
m n 

v “ Z Z vmn  s1n  n,^s1n  nr^mnW  (54) 

m n 

» * U Wmn  cos  m*(sin  £r)<W(t) 


%n(t)  = 

r2"  rl 

Jo  Jo  Jo 


WA-.phTT£l 

on 


v(t> 


w-V 

■f 

'»  «»  T1’ 

^nmi  ^ 

- x)  dx  dx  d<t> 

8 

u2  ~”snn 

on  on; 

phrr 

C fo  C tpx(x'  *’  A)u",n  cos  ** cos  ^ 

+ P.(x,  *,  xjV^  sin  m$  sin  ~ 


+ Pr(x,  4,  A)Wmn  cos  m*  sin 

sin  “mo(t  • *)  d*  <Jx  d$ 
mrr 

X “mnPh  TT  (U2  + VL  * W2  ) 
mn  Z iim  ron  mn 


Solution  for  Radial  Impulse 

a.  Unit  Step  Load  Distributed  over  Area 
(2c,  x 2c2) 

u ■ j,  i,  “ 

* i,  i j. 

n=l  m=l  1=1 


-J,  J,  As-uw-t*!- 

■■j,  V-1'”'-5? 


4 J,  J , l °-(t)  s,n  ^ tos  - 

(5 

where 

W>  ■ 

’ 


°nm1  nml 


CM  ? «t»  T1 


(59) 


s = 1 with  in-plane  inertia  neglected 
s = 3 with  in-plane  Inertia  included 


b.  Solutions  for  Unit  Impulse 

Solutions  for  a unit  impulse  can  be  found 
by  differentiating  with  respect  to  time  the 
solutions  for  a unit  function.  Q-n-ift)  and 
Q„.,(t)  become: 


Qnoi^ 


* Jo. /I.  -j-J-t  ,in  Sfi- sin 
n*7  a l"  “nol  5 1 1 ' 


fsin  "noi(t) 


'V 


"noi 


} 


(60) 


W1* 


8 

phv 


{^T^Ta-rTT"5?5’"5? 

1,5 in  “nm1(th 


S1„  Sfcj.  s1n  |- 


nmr 

"nmi 


7 


(61) 


c.  Dynamic  Strains  and  Stresses 

Having  determined  the  dynamic  displace- 
ments, the  dynamic  middle  surface  strains  and 
stresses  are  determined  from  the  following 
relationships  [1]: 


cx  ■ » 

= li“  + II 

Yx<()  a sg.  ax 

1 — ^~7  (cv  + wt.J 

X t - v X 9 

% = T^V  (s  + vtx} 


nx*  = 2(1  + v)'  YX* 


(62) 
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Data  for  Illustrative  Example 


1,2  Inches 

a ■ 60  inches 

24  inches 

; * 12  inches 

2 inches 

c2  * 2 inches 

0 radians 

v = 0.33333  n * 1 - I 

0-29  P 

c 

(Flugge)  = 5083-855  ^ 

P ,T.  . . . - 5071.793  v&i 

c(Timoshenko) 

E * 30  x 106  ps.i. 

TABLE  I 

EFFECT  OF  EXTERNAL  UNIFORM  PRESSURE  ON 
FUNDAMENTAL  FREQUENCY  {CYCLES/SEC. ) 

Flugge's  Theory 

h/2a  = 0.01  s /2a  =0.2  P*  5083.855 

c ps.i. 


P/Pc 

n,  m 

f 

0 

1.  6 

461.46 

0.2 

1,  7 

422.42 

0.4 

1.  7 

374.72 

0.6 

1.  8 

313.82 

0.8 

1,  8 

230.72 

0.95 

1.  9 

115.63 

0.98 

1.  9 

73.13 

1.00 

1.  9 

0.00 

TABLE  II 

EFFECT  OF  EXTERNAL  UNIFORM  PRESSURE 
ON  LOWER  FREQUENCIES  (CYCLES/SEC.) 

Flugge's  Theory 

n = 1 


P/Pc  * 0 

m 

fl 

f2 

f. 

j 

0 

572.40 

2579.58 

4471.88 

1 

565.22 

2601.76 

4508.18 

2 

545.75 

2666.71 

4615.37 

3 

519.28 

2770.30 

4788.65 

4 

492.44 

2907.31 

5021.12 

5 

471.53 

3072.65 

5304.90 

6 

461.46* 

3261.81 

5632.13 

7 

465.25 

3470.96 

5995.59 

8 

483.95 

3696.87 

6389.04 

9 

517.04 

3936.83 

6807.24 

10 

563.17 

4188.57 

7245.89 

P/Pc  - 0.5 


fl 

f2 

f3 

537.79 

2572.29 

4467.77 

529.15 

2594.27 

4503.96 

505.12 

2658.64 

4610.85 

470.58 

2761.32 

4783.66 

431 .69 

2897.20 

5015.50 

394.67 

3061.23 

5298.54 

365.24 

3248.95 

5624.92 

348.42 

3456.57 

5987.48 

347.99* 

3680.88 

6379.98 

365.62 

3919.18 

6797.20 

400.58 

4169.22 

7234.85 

* Fundamental  frequency 


TABLE  III 

EFFECT  OF  EXTERNAL  UNIFORM  PRESSURE  ON 
HIGHER  FREQUENCIES  (CYCLES/SEC.) 

Flugge’s  Theory 

n = 3 


m 

fl 

P/Pc  * 0 
f2 

f3 

0 

1893.26 

7738.74 

13404.75 

1 

1902.11 

7745.76 

13416.82 

2 

1910.69 

7766.78 

13452.99 

3 

1925.04 

7801.68 

13513.04 

4 

1945.25 

7850.26 

13596.67 

5 

1971.41 

7912.26 

13703.44 

F/Pc  * 0.5 


1807.77  7716.86  13392.16 

1810.44  7723.81  13404.20 

1818.46  7744.65  13440.26 

1831.89  7779.25  13500.14 

1850.85  7827.42  13583.54 

1075.45  7888.90  13690.02 
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TABLE  III  (continued) 


TABLE  IV 


n = 5 

P/Pc  - 0 


EFFECT  OF  EXTERNAL  UNIFORM  PRESSURE  ON 
FREQUENCIES  (CYCLES/SEC.)  IN  8UCKLING 
MODE  COMPARISON  OF  THEORIES 


n> 


0 

1 

2 

3 

4 

5 


5091.63 

5094.85 

5104.51 

5120.61 

5143.16 

5172.16 


12897.91 

12902.10 

12914.67 

12935.60 

12964.84 

13002.34 


22339.84 

22347.08 

22368.80 

22404.95 

22455.46 

22520.24 


Parameters:  h/2a  = 0.01,  t/ 2a  = 0.2 
Buckling  Mode:  n = 1,  m = 9 
Buckling  Pressure:  Pc  = 5083.855  p.s.1. 

(Flugge's  Theory) 

P = 5071.793  p.s.i . 

(Timoshenko's  Theory) 

Flugge's  Theory 


n 

1 


P/Pc  - 0. 

r 

Neglecting 

5 

Including  Axial 

Axial 

fl 

f2 

f3 

P/P 

fl 

Inertia 

f2  f3 

Inertiu 

fl 

4998.41 

12861.43 

22318.82 

c 

L 0 

5001.56 

12865.59 

22326.04 

0 

517.04 

3936.83  6807.24 

519.78 

5011.04 

12878.05 

22347.70 

0.2 

462.47 

3929.78  6803.23 

464.91 

502b. 84 

12898.80 

22383.75 

0.4 

400.52 

3922.72  6799.21 

402.62 

5048.97 

12927.79 

22434.12 

0.6 

327.03 

3915.64  6795.19 

328.74 

5077.44 

12964.97 

22498.72 

0.8 

231.25 

3908.55  6791.18 

232.46 

1.0 

— 

— — 

... 

m - 

10 

p/pc  - o 

P/Pc 

Timoshenko's 
fl  f2 

Theory 

f3 

fl 

f. 

f„ 

f 

0 

517.05 

3937.01 

6807.29 

519.82 

1 

2 

3 

0.2 

462.47 

3933.98 

6807.90 

464.94 

563.17 

4188.57 

7245.89 

0.4 

400.51 

3930.95 

6808.52 

402.65 

1196.55 

6123.37 

10597.00 

0.6 

327.02 

3927.92 

6809.14 

328.76 

2195.92 

8410.40 

14562.79 

0.8 

1.0 

231.24 

3924.90 

6809.75 

232.47 

3600.27 

5414.29 

10830.29 

13310.62 

18756.15 

23052.99 

P/Pc  - 0.5 


400.58 

1077.53 

2088.67 

3498.39 

5315.33 


4169.22 

6100.93 

8382.92 

10796.81 

1327U.66 


7234.85 

10584.18 

14547.01 

18736.89 

23029.97 
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EFFECT  OF  EXTE^-AL  UNIFORM  PRESSURE  ON  DYNAMIC  RESPONSE  FOR  UNIT  IMPULSE 
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Figure  2.  Radial  Displacement  versus  Time  forfor 

* * 0,  x - t/2  at  Various  Values  of  Pressure 
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External  uniform  prassures  of  large  mag- 
nitude and  small  variations  of  impulsive  area 
widths  greatly  affect  the  dynamic  response  of 
cylindrical  isotropic  shells  subjected  to  a 
localized  impulsive  loading. 

For  free  vibrations  pressures  of  large 
magnitude  reduce  the  lower  frequencies  sub- 
stantially while  the  higher  frequencies  are 
not  appreciably  affected.  Uniform  pressures 
in  the  neighborhood  of  50  percent  of  the 
buckling  pressure  can  reduce  the  fundamental 
frequencies  by  30  percent,  while  the  higher 
frequencies  normally  the  membrane  mode 
frequencies  of  the  n,  m mode  will  have  no 
appreciable  change. 

Comparison  of  frequencies  with  the  Flugge 
and  Timoshenko  theories  show  good  agreement  for 
n = 1,  m = 9,  Table  IV;  however  for  low  order 
bending  modes  (n  = 2,  m = 1)  there  may  Indeed 
be  c moderate  difference  in  the  frequencies. 

For  forced  vibrations  as  Illustrated  by 
a localized  unit  impulse,  the  following  con- 
clusions can  be  made: 

a.  Large  external  un  form  pressures  have 
predominantly  large  effects  on  longitudinal 
displacements  and  strains.  Consequently  the 
longitudinal  stresses,  ox  will  be  more  sensitive 
to  change  while  the  circumferential  strains 

and  stresses  will  Increase  moderately. 

b.  Shearing  stresses  also  decrease  as 
demonstrated  in  Table  VII. 

c.  Radial  displacements  and  response 
times  will  have  considerable  Increases  as  shown 
in  Figure  2. 

d.  Increasing  the  width  of  loading  can 
produce  large  changes  in  the  displacements  and 
stresses  as  shown  in  Table  IX. 

e.  Comparison  of  the  theories,  with  or 
without  in-plane  inertias  indicates  that  good 
agreement  exists  between  the  radial  and 
circumferential  displacements,  while  the 
greatest  discrepancy  occurs  in  longitudinal 
displacements  and  strains. 
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ON  DEFINING  TIHE  DOMAINS  FOR 
RADIATION  DAMPING  AND  ADDED  MASS  EFFECTS 
IN  FUIID-STRUCTORAL  INTERACTION 


Alfred  V.  Clark,  Jr. 
Naval  Research  Laboratory 
Washington,  D.C.  20375 


A study  Is  made  of  the  behavior  of  an  oscillator  composed  of  a rigid  sphere 
Immersed  In  an  acoustic  fluid  of  infinite  extent,  and  coupled  to  a fixed  base  by 
an  elastic  spring.  The  sphere  is  impulsively  excited  and  its  subsequent  motion 
calculated.  For  high-frequency  oscillators,  it  is  found  that  radiation  damping 
by  the  fluid  determines  the  form  of  motion  for  all  times.  For  low-frequency 
oscillators,  added  kass  effects  dominate  the  solution  for  times  on  the  order  of 
a period  of  the  oscillator.  The  assumption  of  radiation  damping  is  shown  to 
approximate  the  exact  solution  for  low-frequency  oscillators  at  early-times,  but 
underestimates  peak  response. 


INTRODUCTION 

In  the  study  of  fluid-structural  inter- 
action problems,  attempts  are  often  made  to 
approximate  fluid  effects  due  to  structural 
motion  by  "radiation  damping"  or  "added  mass." 

In  general,  the  statement  is  made  that  radia- 
tion damping  is  valid  for  "high-frequency, 
early-time"  phenomena,  and  added-mass  (or 
accession  to  inertia)  for  "low-frequency,  late- 
time" phenomena.  For  vibration  problems,  where 
transient  effects  can  be  ignored,  one  can  use 
the  ratio  of  structural  to  acoustic  wavelengths 
to  determine  whether  high-  or  low-frequency 
effects  are  dominant. 

For  underwater  shock  phenomena,  "early 
time"  is  often  described  as  being  much  less 
than  the  time  necessary  for  a shockwave  to 
transit  the  body,  so  that  the  time  scale  is 
based  on  geometric  considerations.  However,  it 
can  be  shown  that  there  are  cases  for  which 
radiation  damping  influences  the  motion  of  the 
structure  for  times  much  later  than  early-times. 

To  obtain  estimates  of  the  relative  impor- 
tance of  radiation  damping  for  various  time  (or 
frequency)  regions,  a study  has  been  made  of 
the  response  of  a simple  single-degree-of-free- 
dom  oscillator  submerged  in  an  infinite  acous- 
tic medium,  and  subjected  to  an  impulsive  exci- 
tation. It  can  be  shown  that  the  loading  due 
to  structural  motion  through  the  fluid,  i.e., 
the  "radiation  loading"  can  be  represented  by 
the  convolution  of  the  oscillator  acceleration 
with  a "hydrodynamic  function."  This  hydro- 
dynamic  function  represents  the  load  on  the 
oscillator  when  given  a step-velocity  change, 
so  that  it  is  a generalized  impulse  response. 

For  an  arbitrary  body,  it  is  difficult  to 


find  the  exact  form  of  this  function,  so  vari- 
ous schemes  for  approximating  it  have  been  pro- 
posed; these  schemes  take  advantage  of  the 
early-  and  late-time  behavior  of  the  function, 
and  attempt  to  approximate  its  behavior  at  in- 
termediate times.  Since  the  presence  of  the 
convolution  integral  in  the  equation  of  motion 
of  the  oscillator  Increases  its  complexity,  it 
is  desirable  to  know  u'nen  It  s permissible  to 
use  radiation  damning  and  added  mass  theories 
in  underwater  shock  problems. 

By  examining  the  equation  of  motion  In 
dimensionless ■ form,  it  can  be  shown  that  a rele- 
vant parameter  Is  the  dimensionless  time  (or 
frequency)  ratio  e,  defined  to  be  the  ratio  of 
the  period  of  the  oscillator  to  the  time  for  an 
acoustic  wave  to  pass  around  the  body  perimeter. 
If  e « 1;  i.e.,  for  a high-  frequency  oscilla- 
tor (by  this  definition),  radiation  damping  will 
be  the  dominant  mechanism  in  determining  struc- 
tural motion  for  all  times  of  interest.  For 
e » 1,  low-frequency  oscillators,  added  mass 
effects  will  be  dominant  for  times  which  are  of 
the  order  of  magnitude  of  oscillator  period. 

This  paper  considers  the  case  of  a rigid, 
oscillating  sphere,  since  it  is  possible  to 
solve  t'  r the  hydrodynamic  function  exactly,  and 
examln  .ts  effect  on  oscillator  motion.  It  is 
showr  chat  there  is  no  "boundary  layer"  in  time 
for  which  radiation  damping  is  dominant,  and  is 
to  be  matched  with  added-mass  effects  occuring 
at  a later  time.  Rather,  in  those  cases  for 
which  added  mass  effects  are  dominant,  these 
effects  are  present  at  all  times;  for  cases  in 
which  radiation  damping  is  dominant,  it  will  be 
dominant  at  all  times.  The  frequency  ratio  e 
is  the  factor  which  determines  the  character  of 
the  solution;  it  is  ehown  that  indlscrlmate  use 
of  early-time  radiation  damping  in  conjunction 
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with  low-frequency  oscillators  or  late-time 
added  mass  with  high-frequency  oscillators  can 
lead  to  erroneous  results. 

EQUATIONS  OF  MOTION  AND  INTERACTION 

Consider  a rigid  sphere  submerged  in  an 
acoustic  medium  and  connected  to  a rigid  base 
by  an  elastic  spring  as  shown  in  Fig.  1.  It  is 
assumed  that  the  acoustic  impedance  seen  by  the 


placement  pattern  and  x'  a generalized  coordin- 
ate, the  fluid  force  can  be  related  to  body 
motion  by  the  following  equation  [ll: 

t' 

Ff  ^ “ ■ j X'  (t'  - O X (t')  dT'  (2) 
t'  - 0 

It  can  be  verified  by  direct  substitution 
that  the  Impulse  response  function,  X',  is  pro- 
portional to  the  fluid  force  acting  on  the  body 
when  the  body  Is  given  a step-velocity  change 
in  displacement  pattern  cp  (o');  l.e.,  when 
y'  (o',t')  - 8'  6'  (t')  to  (o').  Except  for 
simple  geometries,  we  do  not  have  exact  expres- 
sions for  X'  (t');  however.  It  Is  possible  to 
genera te  approximations  which  will  give  early- 
time  tadletlui,  and  late-time  added  mass  behnvion 

For  the  rigid,  oscillating  sphere,  we  can 
calculate  X ' by  first  determining  the  pressure 
acting  over  its  surface  due  to  a step  velocity 
change.  When  the  sphere  is  given  the  velocity 
§'  (o',t')  * 8'  (t')  co  (o'),  the  pressure  over 
its  surface  £an  be  computed  from  [ll 

P'  (o'.t'l  - j z ' (o ' , t ' - t')  x'  (t')  dt'(3) 
t'  - 0 

For  the  sphere,  let  y'  (o',t')  » °o  f>’  (t')cosc 

where  0 is  the  angle  between  the  outward  normal 
and  the  x-dlrection  in  Fig.  2. 


Fig.  1.  Rigid  Sphere  Immersed  In  Acoustic  Medium 
and  Connected  to  Fixed  Base  by  Elastic 
Spring 


sphere  is  that  of  an  infinite  medium.  The 
sphere  is  initially  at  rest  and  is  given  an 
impulse  at  t'  - 0*.  The  resulting  motion  of 
the  sphere  is  determined  from  the  equation 

m'  x'  (t')  + lc'  x'  (t')  = + Fj  (t')  + 

I'  (t')  (1) 

where  m'  ■ mass  of  sphere 

x'  » displacement  of  sphere 
k'  * stiffness  of  spring 
Ff  « fluid  force  on  sphere 

1^  “ Impulse 

6'  « Dirac  delta  function. 

Obviously,  we  need  the  value  of  Fj,  to  specify 
the  problem  completely. 

The  fluid  force,  Ff,  arises  due  to  the 
sphere's  motion  through  the  fluid.  Let  c ' 
denote  the  surface  of  the  body,  and  y'  denote 
displacement  of  the  body  along  its  outward 
normal.  If  y'  » cp  (o')  x'  (t')  with  cp  a dis- 


y* 


SPHERE 


Fig.  2.  Defines  Angle  0 

Then  the  associated  pressure  Pg  (o',t')  is 

Pj  (c',t')t- 


I- 


(a'.t'-r')  6'  (t')  dr' 


(4) 


*We  use  primes  to  denote  dimensional  quantities; 
unprlmed  quantities  are  dimensionless. 


z 

- 0 


which  shows  that 

z'  (c',t')  “ P'  (o', t')  / v'. 

The  surface  pressures  fo?  the  sphere  given 
a velocity  step  have  been  calculated  to  be 

* * 

Pj  (o ' , t ' ) = ?c  v^  c w 1 cos  (w*'  t ') (cos  fl)  5) 

where  Wj*  ■ c/a,  c « sound  speed  in  fluid,  and  a 
is  the  sphere  radius  [2]. 


F = F'/(m'a(w^)a) 

e “ w,'/w' 
o 

x « dx/dt 

* i ' 

m » m /m 
w 

The  system  of  equations  (1)  and  (9)  assume  the 
dimensionless  form 

x (t)  + x (t)  - F (t)  = Iq5  (t)  (10) 


The  force  associated  with  this  pressure  is 
computed  from 

F'  (t')  “ F*  (o',t')  cos  0 do'  (6) 

* 

where  the  integration  is  performed  over  the 
sphere  surface;  from  equations  (2),  (3),  and  (5), 
it  is  found  that 

X'  (t')  « ££  z'  (o',t')  cos  6 do' 
for  the  oscillating  sphere.  Consequently, 

X ' (t')  « -^a-  pc  e 1111  cos  v't ' H (t')  (?) 

with  H(t')  the  Heaviside  step  function. 

Substituting  this  result  into  equation  (2), 
taking  the  Laplace  transform  of  (2)  for  zero 
initial  conditions,  and  rearranging  leads  to: 


o 

F (t)  + 2e  F (t)  + 2c*  F (t)  - 

-m  « [°x°+  t XT  (11) 

Note  the  quantity  ■ t is  a measure  of  the 
coupling  between  the  fluid  and  the  sphere;  thi 
quantity  has  several  physical  irterpretations 
Writing 

me*  (®^)  / (m'w')  * R 

shows  that  R is  the  ratio  of  two  impedance-like 
terms.  The  denominator  represents  the  imped- 
ance of  the  sphere-spring  system  in  vacuo, 
whereas  the  numerator  is  the  early-time  imped- 
ance of  a piston  with  surface  area  equal  to  1/3 
of  the  sphere  surface.  Hence,  R is  the  ratio 
of  the  "fluid  impedance"  to  the  "oscillator 
impedance  in  vacuo."  For  simplicity,  wc  let 
m ~ 0 (1)*  in  this  paper. 


+ 2wj' s ' 


+ 2wi3 


F'  (s')  * 


^ (s')a  (s'  + Wj' ) x'  (s';  (8) 

with  A = the  area  of  the  sphere. 

This  can  be  inverted  to  give  the  differential 
equation 

F'  (t')  + 2w1'F'(t')  + 2w*  F'(t')  - 


provided  that  appropriate  initial  conditions 
vanish.  We  can  guarantee  this  by  letting  our 
system  respond  at  t'  * 0 , a short  time  after 
the  application  of  the  impulse. 

Equation  (9)  is  the  interaction  equation 
which  relates  fluid  force  to  oscillator  motion. 
Combined  with  (1),  it  completely  defines  the 
problem,  giving  rise  to  a 4th  order  set  of 
equations. 

It  is  convenient  to  cast  the  problem  in 
dimensionless  form;  to  this  end,  we  set 

w^  « .fk'/m',  natural  frequency  In  vacuo 

t = w't' 
o 

x = x'/a 

m'  = paA/3 

io  - i ;/(■»'«*;) 


Since  wf  is  the  frequency  (and  decay  con- 
stant) for  X'(t'),  we  have  another  interpreta- 
tion of  R.  The  quantity  e measures  how  fast 
\'(t)  is  oscillating  (and  decaying)  relative  to 
the  sphere-spring  oscillator.  When  e is  small, 
then  « w'  and  we  have  a "high-frequency 
oscillator."0  X'(t')  remains  essentially  con- 
stant (a  Heaviside  step-function),  while  the 
oscillator  executes  many  vibrations.  For  large 
g,  we  have  a low-frequency  oscillator;  X'(t') 
oscillates  (and  decays)  rapidly,  while  the 
oscillator  is  slowly  responding. 

In  the  literature,  it  is  sometimes  stated 
that  the  submerged  body  moving  through  the 
fluid  radiates  plane-waves  fot  "ear’y-times," 
or  that  the  fluid  can  be  represented  as  an 
added-mass  for  "late-times."  It  is  proposed  to 
deduce  from  the  equations  of  motion  and  inter- 
action for  the  sphere  what  the  effect  these 
phenomena  have  upon  oscillator  motion.  Since 
the  parameter  e is  a measure  both  of  frequency 
and  fluid-oscillator  coupling  in  the  interaction 
equation  (11),  it  is  natural  to  consider  what 
happens  to  this  equation  (and  the  equation  of 
motion  (10))  when  e becomes  asymptotically 
small  or  large.  That  is,  we  look  at  the  be- 
havior of  the  equations  for  high-  and  low-fre- 
quency oscillators.  We  will  find  the  exact 
solution  for  oscillator  displacement  x by 
taking  the  Laplace  transform  of  the  coupled  set 
of  equations  (10)  and  (11)  then  rearranging  to 
obtain  an  equation  for  x;_this  will  have  the 
form  x = N(s)/D(s),  with  x the  Laplace  transform 

*The  notation  m ~ 0 (1)  means  that  (m)  is  of  the 
order  of  magnitude  of  unity. 
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of  x.  We  use  the  residue  theorem  to  solve  for 
x;  to  do  so,  we  must  find  the  roots  of  the 
denominator  D(s).  We  will  discover  that  the 
location  of  the  roots  is  easily  estimated  for 
c » 1 and  e « 1;  i.e.,  for  low-  and  high- 
frequencv  oscillators.  One  set  of  roots  will 
be  of  order  of  magnitude  unity;  the  other  set 
will  be  of  order  €.  The  first  set  of  roots 
correspond  to  vibrations  on  the  time  scale  of 
the  oscillator  period  in  vacuo.  The  second  set 
of  roots  correspond  to  motion  on  what  we  call 
the  "intermediate-time"  scale,  t'  ~0(a/c).  We 
will  show  that  the  residues  associated  with  ttis 
latter  set  of  roots  become  vanishingly  small 
for  low-  and  high-frequency  oscillators.  In 
the  asymptotic  analysis  section,  we  show  why  wc 
should  expect  this  behavior.  Furthermore,  we 
show  that  using  radiation  damping  to  predict 
the  response  of  low-frequency  systems  is  only 
valid  for  early-times,  and  that  maximum  oscilla- 
tor excursions  occur  (for  € » 1)  at  late-times. 
The  early-time  damping  theory  will  be  shown  to 
underpredict  peak  oscillator  displacement  by  a 
factor  proportional  to  1/e,  for  € 1. 

EXACT  SOLUTION  FOR  OSCILLATING  SPHERE 

Returning  to  the  system  of  equations  (10) 
and  (11),  taking  the  Laplace  transform  for  zero 
initial  conditions,  and  combining  the  two  trans- 
formed equations  results  in 

I r(s  + e)a  + eE  T 

x (S)  = -2 (12) 

[(s3+l)(sa+2es+2e9)+sa  (m)t(s+e)] 

Here  s = s'/w',  s'  being  the  (dimensional) 
Laplace  trans?orm  variable,  and 

or 

x (s)  = l e”St  s ( t)  dt  (13) 

t = 0 

We  rewrite  our  equation  on  x as 
I r(s+e)a  + ea  1 

x (s)  = (14) 

(s-Sj  ) (s-sr)  (s-sa)  (s-s4) 


This  gives  us  an  undamped  sinusoidal  vib- 
ration for  low-frequencies  at  late-times.  The 
ratio  of  this  submerged  frequency  to  the  in 
vacuo  frequency  is: 


;UD  » [^2  + m I , which  shows  how  the 

o 

"added  mas3"  effect  lowers  the  frequency. 

Now  let  us  consider  the  low-frequency  case 
at  intermediate  times;  we  have  t'  ~ 0(a/c),  so 
that  t — 0(w'a/c),  or  t — 0(l/«).  For  motion 
to  occur  on  ?he  intermediate  time  scale,  we 
require  that  s — 0(e).  The  denominator  becomes, 
approximately, 

s‘  + e s3  (2  + m)  + eas*  (2  + m)  =s  0 (17) 

which  has  roots 

slts  - - | [_(2  + m)  ± {(2  + m)a  -A  (2  + m)}^ 

(18) 

If  the  mass  of  displaced  water,  m'  , is  greater 
than  2m',  then  s,  ,8  will  be  real;  this  gives 
exponentially  decaying  solutions  to  the  equation 
of  motion  at  intermediate-times.  If  m < 2m', 
then  the  sphere  will  experience  damped  vibration 
for  intermediate-times. 

Consider  the  high-frequency  case  e « 1, 
and  let  us  look  at  intermediate-time  solutions; 
s ~ 0(e)  or,  s « 1.  The  denominator  reduces 
to* 

s3  + 2es  + 2ea<w0ors  . » - e ± is  ; (19) 

3 •* 

we  see  that  intermediate-time  solutions  corre- 
spond to  lightly-damped,  slow  oscillations.  For 
early-times,  t'  « a/c;  in  the  case  of  high- 
frequency  oscillators,  w^  » a/c.  Thus,  early- 
times  for  high-frequency  oscillators  are  those 
for  which  t ~ 0(1),  Motion  occurring  on  this 
time-scale  will  be  associated  with  roots  for 
which  s ~ 0(1).  For  e <<  1 and  s ~ 0(1),  we 
have  : 

s3  + s m s + 1 0 (20) 


where  st , ss  etc.,  are  roots  of  the  denominator 
in  (12). 

Now  late-times  are  those  for  which 
t'  » a/c;  recall  that  our  dimensionless  time 
t is  given  by  w't'.  For  lov?- frequency  oscilla- 
tors, late  times  will  correspond  to  t ~ 0(1), 
since  w'  « a/c;  that  is,  late-times  correspond 
to  times  oi.  the  order  of  magnitude  of  the 
oscillator  period.  Thus,  for  low-frequency 
oscillators  at  late  times,  we  set  £ » 1 and 
s — 0(1)  in  (12);  the  denominator  in  (12)  will 
then  be  approximated  by 

2es  (sa  + l)  + sV  (m)  a 0 (15) 

which  gives  us  the  roots 


which  has  roots 

Sj  ,2  = - h [e  m ± {ca  m - 4}  (21) 

Depending  on  the  sign  of  the  discrim' nant,  we 
can  have  underdamped,  critically  damped,  or 
over-damped  solutions.  However,  since  we  have 
stipulated  that  e « 1,  it  will  take  a large 
mass  ratio  m to  give  us  an  overdamped  solution, 
even  for  early-times. 

Roots  for  tnv  specific  case  of  m = 2.5  are 
given  in  Table  A.  The  early  time  roots  for 
e « 1 are  given  as  s5  and  Sj , as  are  intermedi- 
ate-time roots  for  e » 1;  intermediate  and  late- 
time  roots  for  c « l and  £ » 1 are  s and  s4  , 
These  are  in  good  agreement  with  the  asymptotic 


* We  assume  se  m « 1,  £a  m •'<  1 
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roots  listed  above 


early-time  solution 


For  the  high-frequency  cases  (e  * .01,  .10) 
we  have  a solution  in  the  form  of  decaying  sin- 
usoids. 


X (t)  - xe  (t)  + (t) 

I al  t 

- ^ 6 [Qi  sin  bj 
+ eS3  C [q  sin  b 

b 3 3 

3 

where  we  have  set 


t + Qj  cos  bj 
t + Q.  cos  b 

3 


(22) 

t] 

t] 


X ’ To  e'mC't/2  sin  t H(t-0+)  - >.g 

A plot  of  (22)  is  shown  in  Fig.  1 for  ? = ,01 
and  IQ  = 1;  it  is  seen  that  X does  indeed  have 
the  form  above. 


° o 


08- 


o o 


«=001 


Sj  ” aj  + ibj  s3  = a-  - ibj 

s » a + ib  s.  = a - ib 

3 3 3 3 3 

and  the  constants  through  Q4  are  defined  in 
Table  B.  Note  that  aj  and  a are  negative. 
Here  X (t)  is  the  early-time  solution  corre- 
sponding to  roots  Sj , st;  X^  (t)  is  the  inter- 
mediate-time solution. 

For  the  intermediate-  (e  “ 1.0)  and  low- 
frequency  cases,  we  have 
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(23)  Fig.  3.  High-Frequency  Oscillator  Response 


I0[*j  + e3gt 

(a3  -aj  ) (4  + ) 

I e33t  [g  sin  b t + h cos  b tl 


3 


The  coefficients  are  defined  in  Table  C;  X^  is 
the  intermediate-time  solution,  and  Xf  is  the 
late-time  solution. 

For  the  high-frequency  oscillator,  we  can 
show  that: 


bj  = 1 b - s (24) 

3 

Qi  /t>!  = 1 Qs  = % “ Qs  = o 

We  see  from  (21)  that  the  conditions  on  aj  and 
bj  are  fulfilled.  To  show  e.g.  that  Q4 
vanishes,  consider  that  Table  B gives  Q,  as 

Q.-,  “ [2(aj+e)biH1  + Aj  Jj  (es  -b?+(aa  te)s  ) 1/M,  (23 

Since  a = -e  and  b = e (as  shown  by  (19)),  we 
see  thal  Q*  vanishes;  similar  results  hold  for 
Qj  and  0 . It  can  also  be  verifiea  that 
Qj /b^  = 1,  reducing  (22)  to  the  high-frequency, 


In  the  instance  of  the  low-frequency 
oscillator,  we  wish  (23)  to  reduce  to  the  added- 


10' 

X 

08 

- 

o 

o 

% 

06 

o 

o 

o 

l_ 

0 

0 

0 

04 

h 

02' 

A 

1 0 

0 

- 1 - i 1 

0 

1 

1 

02ir  ir  2ir 

3»r 

4 w 

-02 

0 

0 

-04 

0 

-06 

0 

0 

-08 

- 

-to1 

_ 

Fig.  4.  Low-Frequency  Oscillator  Response 
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ataa,  jtte-tiae  solution  giver  by  I b (in  b t* 
H (t-0  Now  the  roots  given  in  (it;  show3 

that  a becomes  vanishingly  small,  and  \?  equals 
(2/(2+«)  for  e » 1.  Using  these  results,  we 
find  that  the  coefficients  multiplying  the  ex- 
ponentials e*1 1 and  e in  (23)  are  0 (Io/e); 
the  intermediate-time  solution  then  becomes 
vanishingly  small  as  e becomes  large.  Further- 
more, it  can  be  demonstrated  that  h - 0,  so 
that  the  late-time  solution  depends  only  upon 
(1  /b  ) g sin  b t.  For  the  specific  case  of 
m 2 2^5,  we  have  b?  « 4/9;  also,  g - 9/20.25  as 
e » 1 or  g * l?  so  that  (I  g)/b  - 1 , and 

we  have  reduced3our  general  solution  ?23|  to  the 
low-frequency,  late-time  form.  The  behavior 
of  (23)  is  shown  in  Fig.  4 for  e * 100  and 
I * 1.  In  Fig.  5,  we  show  X and  for  « » 10 
over  the  time  domain  for  which  they  are  the 
sam<-  order  of  magnitude;  this  clearly  indicates 
that  our  solution  is  dominated  by  late-time, 
added-mass  effects. 
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Fig.  6.  Early-time  and  Exact  Solutions  for 
Low-Frequency  Oscillator 

the  maximum  oscillator  excursion,  from  the 
exact  solution,  is  given  by  I b.  Hence,  if  we 
were  to  try  to  predict  peak  oscillator  displace- 
ments using  radiation  damping,  we  would  under- 
estimate the  peak  by  a factor  of  1/bme  for 
<5  » 1. 
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Fip  3..  Intermediate-  and  Late-Time  Solutions 
for  Low-Frequency  Oscillator 

It  is  interesting  to  compare  the  exact 
solution  £'jf  c » 1 with  approximate  solutions 
generated  by  assuming  early-time  radiation 
damping.  In  this  case,  our  interaction 
equation  (11)  reduces  to  F(t)  “-me  X(t) , 
giving  ar  equation  of  motion  of  the  form 

X(t)  + m e J(t)  + X(t)  - 1Q  6(t)  (26) 

For  large  e,  the  solution  to  this  ig  given 
(approximately)  by  X * IQ  (1  -e  1,15  )/(me). 

In  Fig.  6,  we  show  the  exact  solution  X,  (from 
(23)),  and  X plotted  for  early-times.  Note 
that,  although  gotd  agreement  exists  between 
the  two  for  early-times,  these  times  are  much 
shorter  than  the  time  necessary  for  the  os- 
cillator to  reach  its  peak  excursion.  Now  X 
has  an  upper  bourn  given  by  Io/rae,  whereas 


ASYMPTOTIC  ANALYSIS 

We  previously  noted  how  the  poles  of  our 
expression  for  X(s)  could  be  found  by  assuming 
motion  on  time  scales  t —0(1)  and  t — 0(1  /«)  . 
Note  that,  although  the  denominator  in  (12)  is 
a quartic,  this  assumption  allowed  us  to  sim- 
plify the  determination  of  roots  of  the  denom- 
inator to  the  problem  of  finding  roots  of  a 
quadratic  (see  e.g.  equations  (15),  (17),  (19), 
and  (20)),  We  could  have  further  reduced  our 
labor  in  finding  the  exact  solution  had  we 
known  that  residues  associated  with  poles  for 
which  s — 0(s)  would  become  vanishingly  small. 

We  will  use  some  ideas  from  perturbation  theory 
to  show  why  our  system  of  equations  (10)  and 
(11),  (which  give  rise  to  (12))  exhibits  this 
behavior. 

We  will  consider  the  behavior  of  the  coupl- 
ed system  of  equations  (10)  and  (11)  for  various 
time  scales.  We  make  the  assumption  that,  for 
each  time  scale,  all  the  time  derivatives  of  F 
are  the  same  order  of  magnitude;  a similar 
assumption  is  made  concerning  derivatives  of  X. 

We  then  solve  the  coupled  system  (for  each  time 
scale)  end  examine  the  character  of  the  solution. 
If  it  ii  a physically  meaningful  solution,  our 
assumptions  above  are  valid;  physically  impos- 
sible solutions  imply  that  our  assumptions  are 
Invalid.  In  the  former  case,  we  have  derived  a 
set  of  equations  which  govern  the  behavior  of 
the  system  for  that  particular  time  scale.  In 
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the  latter  case,  our  systeai  of  equations  is 
discarded  as  inapplicable.  We  will  use  a scale 
factor  s,  defined  as  follows:  We  let  t • st, 

with  t ~ 0(1)  and  t * w^t';  thus  s determines 
the  order  of  magnitude  of  (dimensionless)  time 
we  are  considering.  We  will  discover  that  the 
magnitude  of  s for  which  we  derive  meaningful 
solutions  corresponds  to  the  magnitude  of  the 
roots  of  the  denominator  in  equation  (12). 

Low-Prejueocy  Oscillators 
Case  1.  Late- tines 


which  states,  as  expected,  that  F « -m  eX  on 
this  time  scale. 

Substituting  this  into  our  equation  of 
motion  (10)  gives 

?— +m*s  — + X-I  6(t)  (28) 

dP  dt  ° 

and,  assuming  all  derivatives  are  the  ^ame 
order  of  magnitude  and  recalling  that  s » e, 
we  can  rewrite  (28)  in  the  fora 


Considering  the  low-frequency  case  first, 
we  look  at  times  ouch  greater  than  the  time  for 
a sound  wave  to  travel  arcur<d  the  body.  For 
e » 1,  this  corresponds  to  a time  scale 
t -0(1);  i.e.,  we  consider  tlsies  on  the  order 
of  a few  periods  of  vibration  o»  the  oscillator. 
Looking  at  the  interaction  equation  (11),  we 
assume  all  derivatives  of  F wisi  respect  to  t 
are  the  same  order  of  jgagnltude,  as  are  all  de- 
rivatives of  x;  i.e.,  F — 0(F)  and  °X°  ~ 0 ()U. 
Since  e » 1,  this  allows  us  to  simplify  the 
interaction  equation  to  read  F - -o  X.  Recall- 
ing that  F(t)  • -fX(t-T)  X°  (t)  dt  and  that  the 
added  mass  ra'  is  *(l/2)  m'  for  the  sphere,  we 
see  that  X(tf  w (a'Jm')  for  low-frequency 
oscillators  on  this  time  scale.  We  have  re- 
covered our  low-frequency,  late-time  approxima- 
tion on  time  scale  t ~ 0(1). 


dp 


»s  — — I 6(t) 
dt  0 


(J2) 


(i-e  >e)  with 


The  solution  is  given  by  X _ , 

3 * n t ; this  physically  impossible  solution 
corresponds  to  the  solution  of  (26)  and  we  have 
shown  its  behavior  in  Fig.  6.  Since  our  solu- 
tion describes  a physically  unreasonable  dis- 
placement of  the  sphere,  our  coupled  set  of 
equations  will  not  admit  a solution  which 
varies  on  the  early  time  scale.  We  conclude 
that,  although  radiation-damping  is  the  dominait 
fluid  mechanism  at  early-tlmes,  these  times  are 
too  short  to  significantly  affect  the  oscill£ar 
motion.  Furthermore,  we  cannot  attempt  to 
"match''  this  early-time  solution  to  a late- 
time  solution,  since  the  early-time  solution 
approaches  the  value  X£  « I Jot. 


Inserting  this  form  of  F into  the  equation 
of  motion  gives 

(24m)  8°(t)  + 2 X(t)  - 2 I 6(t) 

o 

for  t ~ 0(1);  this  has  a solution  XI  b • 
sin  bt  H(t-0rf)  with  \?  ■ 2/(2  + m) . We  use  the 
notation  Xj  to  denote  the  solution  to  our 
equation  of  motion  on  the  late-time  scale.  Note 
that  Xj  begins  at  t ■ 0+  to  allow  satisfaction 
of  the  initial  conditions. 


For  times  other  than  late-times,  we  define 
a new  time  variable  £ by  t " s t,  where  s is  a 
dimensionless  scale  factor.  Substituting  this 
into  the  Interaction  equation,  we  have 


?■—  + 2cS  ^ + 2eaF 


d? 


d£ 


-a  r-  d3F  a.  4a F * 
iff  [s  + s JC7) 

dt3  dP 


We  will  examine  the  behavior  of  this  equation 
for  various  orders  of  magnitude  of  s. 

Case  2.  Early-times 


Defining  early-times  as  those  for  which 
radiation-damping  is  the  dominant  fluid  mech- 
anism, so  that  -Fj  «*  (pcA/3)  x'(t),  requires 
that  X'  * (pcA/3)  H't')  for  early-times.  We 
see  that  X'(t’)  will  have  this  asymptotic  form 
when  t’  « a/c,  or  when  our  dimensionless  time 
t is  of  order  1/s  with  s » «.  Our  interaction 
equation  (27)  then  has  the  approximate  form 

d*F  « dsX 

— * -ra  €8  

d»t  d£3 


Case  3.  Intermediate  times 

This  motivates  a search  for  a solution 
where  time  derivatives  are  the  same  order  of 
magnitude  for  intermediate  times;  i.e.,  when 
t ~0(l/e),  or,  in  terms  of  dimensional  time, 
t'  ~ 0(a/c) . 

We  take  £ « st,  where  now  s ~ 0(e).  Recall 
that  our  late-time  solution  X.  was  given  by 
Iob  sin  bt,  which  will  be  non-zero  even  for 
t ~ 0(1/ e) ; we  must  account  for  X^  and  its 
derivatives  in  (10)  and  (27).  Accordingly,  we 
set  X ■ Xj  + X where  X^  will  be  the  solution 
for  intermediate  times,  and  will  vanish  for 
late-times,  as  we  will  subsequently  show. 
Combining  (27)  and  (10)  gives  the  fourth-order 
differential  equation 

s4  — - + (24m)  es  — - + (24m)  t3  s* — - 4 

d£4  df  dP 


P — + 2eS  ^ 4 2eaX  = 
dP  d£ 


I (P  - 4 2ti--4  2ca  6 (t) ) ; (29) 

° d?  d£ 

with  the  choice  X = X(  4 Xj,  we  reduce  this  to 
a simpler  form 
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dX. 

d®X 

d*X 

s*  -*■  + 

(2-Hn) 

es3  + (2+m) 

.*?— U 

dp 

dt3 

dP 

(29a) 

I (P 

2 ss  — + 2ca6) 

0 

dP  + 

dt 

This  has 

a solution  X.  = (1  It) 
i 0 

- . -P.  t 

.a,  e * 

+ a2  e*Pst!  H(t-0+) 


where  X^  Is  unknown  as  yet.  For  this  time  scaltj 
It  Is  straightforward  to  verify  that  all  terms 
In  the  Interaction  equation  must  be  retained, 
provided  we  made  the  assumption  that  all  deri- 
vatives of  F land  all  derivatives  oe  are  the 
same  order  of  magnitude.  Combining  the  Interac- 
tion and  motion  equations  leads  to  (29).  This 
c;  n be  simplified  by  recalling  that  X s X + Xj, 
and  that  X satisfied  the  differential  equation 
(32).  Using  this  fact  allows  us  to  reduce  (29) 
to  the  form 


where  Pi  and  P2  are  roots  of  the  quadratic 

P*  + (2  + m)  eP  + (2+m)  e*  - 0;  (30) 

which  Is  the  same  as  equation  (17), 

It  can  be  shown  that  the  solution  corresponding 
to  P3  ■ 0 vanishes  Identically.  Note  that  X^-0 
as  t ~ 0(1).  Since  X^  does  not  decay  exponent- 
ially as  Intermediate  times  are  approached,  we 
cannot  define  a region  in  clme  where  X.  and  X; 
are  to  be  matched.  Furthermore,  X^  and  X.  are 
the  same  order  of  magnitude  for  t ~ 0(l/ej,  so 
that  no  "boundary  layer"  in  time  exists*,  Fig.  5 
compares  X,  and  X over  Intermediate-times  for 
e - 10. 

The  constant  a{  and  a2'  are  0(1)  so  that  X^ 
Is  0(lo/c)  and  becomes  vanishingly  small  for 
low-frequency  oscillators,  t » 1.  Our  solution 
reduces  to  the  modified  Impulse  response 
I b sin  bt;  the  intermediate-time  solution,  which 
arises  due  to  the  complicated  fluid-structural 
Interaction  on  time  scale  t ~ 0(l/e),  can  be 
neglected. 

High-Frequency  Oscillators 
Case  1.  Early-tlmes 


d*X,  dX 

? — + 2cS  — + 2 ea  X,  * 0 (e3)  (33) 

dP  dt  1 

Note  that  the  characteristic  equation  for  (33) 
is  the  same  as  (19).  From  the  above,  we  con- 
clude Xj  — 0(e)  so  that  X^  - 0 for  e « 1. 

Case  3.  Late-tlmes 


There  Is  no  late-tlme,  added  mass  solution 
for  high-frequency  oscillators.  This  can  easily 
be  verified  by  looking  at  time  scale  t ~ 0(l/s) 
with  s « e;  on  this  time  scale,  X and  Its 
derivatives  vanish  due  to  the  exponential  term 
e ’ . Then  the  inter  action  equation  re- 

duces to  the  added-mass  form 


F « - h m d*Xf/d? 
where  X.  Is  the  late-tlme  solution, 

7 


We  assume  X.  varies  on  the  late-tlme  scale 
sc  that  all  Its  derivatives  are  the  same  order 
of  magnitude.  The  equation  of  motion  becomes 


t e 

? (l+m/2)  + X.  + P + Xa  = I 8(t) 

dP  6 ° (J4) 


The  high-frequency  oscillator  corresponds 
to  e <•'  1,  or  wf  <•'  w*.  Setting  t » st  for 
early-tlmes  again  gives  the  Interaction 
equation  (27). 

With  s - 0(1)  and  e « 1,  (27)  reduces  to 

F = -m  es  — (31) 

dt 


and  since  X satisfied  the  early-time  equation 
(32),  we  have 

dX  <fx 

X,  **  (m/2)  c.3  — 1 + 0 (?  -)  (35) 

' dt  dP 

We  conclude  that  X.  ~ 0 (P);  for  s « t,  we  see 
that  X,  a 0> 


This  shows  that  the  fluid  force  on  the  sphere 
Is  due  to  radiation  damping  at  early-tlmes  for 
the  high-frequency  oscillator.  Using  (31)  In 
the  equation  of  motion,  we  have 

^ + me  — + X = I 6(t)  (32) 

d?  dt  0 

with  solution  X = IQe  (rae^)t  spn  t H(t  - 0+) . 
The  characteristic  equation  for  (32)  Is  the  same 
as  equation  (20) . 

Case  2.  Intermediate-times 

For  intermediate-times,  we  take  t =*  s t 
where  now  s - 0(e);  we  also  let  X = Xg  + X^, 


Then  we  are  left  with  only  our  early-tlme 
solution;  we  conclude  that  although  radiation- 
damping Is  the  dominant  fluid  mechanism  for 
early-times  only,  Its  effect  determines  the  form 
of  solution  for  much  longer  times. 

Recapitulating  the  above  results; 

1)  For  late-tlmes  at  low-frequencies  (s  —0(1), 
e » 1),  added-mass  effects  are  dominant  In  the 
Interaction  equation;  t —0(1),  or  late-tlmes, 
correspond  to  times  comparable  to  a period  of 
the  oscillator. 

2)  For  early-tlmes  at  high-frequencies  (s  — 0(1)  > 
e « 1),  radiation  damping  Is  dominant  In  the 
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Interaction  equation,  and,  since  t ~ 0(1), 
early-times  are  comparable  to  a period  of  the 
oscillator. 

3)  For  intermediate-times  at  high-  and  low- 
frequencies,  it  is  necessary  to  set  s ~ 0(e), 
and  t — 0 (1/e).  All  terms  are  retained  in  the 
interaction  equation,  so  that  radiation  damping 
is  not  dominant  for  low-frequency  oscillators 
at  intermediate-times,  and  added-mass  effects 
are  not  dominant  for  high-frequency  oscillators 
at  intermediate-times. 

4)  Although  radiation  damping  is  the  dominant 
fluid  mechanism  for  early-times,  we  found  that 
coupling  the  early-time  interaction  equation 
with  the  early-time  form  of  the  equation  of 
motion  led  to  a physically  impossible  solution 
for  low-frequency  oscillators.  The  late-time 
added-mass  dominated  form  of  the  coupled 
equations  was  found  to  be  0 (?)  for  high  fre- 
quency oscillators;  we  had  s « e « 1 in  this 
case. 

5)  For  the  low-frequency  oscillator,  the 
dominant  solution  occurs  on  the  late-time  scale 
when  added-mass  theory  is  valid.  It  is  errone- 
ous to  attempt  to  use  radiation-damping  to 
account  for  fluid  effects  for  this  oscillator. 
There  is  no  boundary  layer  in  time. 

6)  For  high-frequency  oscillators,  the  motion 
is  accurately  treated  by  assuming  radiation- 
damping as  the  dominant  fluid  mechanism;  even 
though  this  is  strictly  true  only  for  e.rly- 
times,  its  effect  on  the  oscillator  is  felt  for 
much  longer  times. 

CONCLUSION 

We  have  shown  that  the  effect  of  the  fluid 
can  be  calculated  by  convolving  the  hydrodynamic 
function  with  the  acceleration  of  the  submerged 
oscillator.  This  convolution  can  be  approximXed 
by  radiation-damping  for  early-times,  defined 
by  t'  « a/c,  or  times  much  less  than  a transit 
time;  for  late-times,  if  the  oscillator  is 
low-frequency,  then  added-mass  effects  prevail. 

For  a given  mass  ratio,  we  use  the  value 
of  dimensionless  parameter  t “ c/(aw)  to 
estimate  whether  our  oscillator  is  h?gh-  or  low- 
frequency.  We  found  that,  by  appropriately 
scaling  the  dimensionless  time  and  keeping  in 
mind  the  magnitude  of  e,  we  were  able  to  couple 
our  interaction  equation  and  motion  equation  to 
give  meaningful  results.  In  this  way,  we  dis- 
covered that  we  cannot  use  the  early-time  ra- 
diation-damped form  of  the  interaction  equation 
in  conjunction  with  the  motion  equation  for 
e » 1,  low-frequency  oscillators.  The  correct 
form  was  shown  to  be  the  intermediate-time  form 
of  the  interaction  equation,  in  which  no  simpli- 
fication was  possible.  The  same  result  was 
shown  to  be  true  for  late-times  when  we  con- 
sidered high-frequency  oscillators.  We  defined 
intermediate-times  as  those  for  which  t'  ~ 0(a/(), 
the  time  for  a wave  to  transit  the  body. 


It  vus  demonstrated  that  the  intermediate- 
time  solution,  derived  from  Che  coupling  of  the 
full,  unsimplified  form  of  the  interaction 
equation  with  "he  equation  of  motion,  became 
vanishingly  small  for  both  high-  and  lev.— fre- 
quency oscillators.  This  left  us  with  the 
early-time,  radiation  damped  solution  for 
t « 1,  showing  that  radiation  damping  made  its 
effect  felt  long  after  early-times  for  high- 
frequencies.  For  e >>  1,  the  late-time  added- 
mass  solution  is  dominant  over  virtually  all 
times;  if  we  wish  to  find  maximum  oscillator 
excursions  (for  maximum  stresses,  etc.),  it  is 
sufficient  to  use  added-mass.  The  urt.  of 
radiation-damping  in  conjunction  with  low-fre- 
quency oscillators,  though  accurate  for  early- 
times,  led  to  underestimates  of  peak  oscillator 
response  by  a factor  of  1/e.  We  also  showed 
that  the  exact  solution  for  our  sphere  model 
reduced  appropriately  to  high-  and  low-frequen<y 
solutions  predicted  by  our  asymptotic  analysis. 
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a + ibo 

3 


* a - ib 

3 3 


-15.0371,  sa  - -29.9630 

- -.0000  ± .6671 
-150.0037,  sa  - -299.963 

- -.000  ± .6671 


*■“**.-  b33 
*ia  “ + Ao 

g - (FaHla  - G3b3Ala)/(A»s  + b3a)(Ae3  + b3a) 
h - (G3H,a  + b3AlrF3)/(Ala  + baa)(Afa  + b^) 


(For  definitions  of  A^  , Ag,  etc.,  see  Table  II) 


TABLE  II 

Constants  for  High-frequency  Solution 

»i  “ *i  + ibi  - «i  * i*h 

s - • + ib  s4  * a - ib 

3 3 3 3 3 


a +t 

3 


*1  - «1  - 8, 


a - a, 

3 1 


®i  “ + b3 

B •>  b.  + b 


b,  -b 

aty3  + Ca  - bj3 

or  a + c3  - b3 

3 3 

2*1  th 


H,  - Al3  - BjBj 
H - A 3 - B 0 

3 3 3 3 

-T1  ■ Bl  + Bj 

J3  ” B3  + S3 

Ml  “ (Al3  + 0i3)^a  + Bj3 ) 
M ■ (Aa+Ba)(Aa  + Ba) 

3 3 3 3 >1 

Oi  ■ (Fi«3  - a3g,.j3)/m3 

Qs  “ (Gi»  ♦ FXA  J )/M 
% ' (Fg11!  * A1G3J1)/Ml 
Q,  - (G3Hl  + F3AaJl)/Ml 
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ENVIRONMENTS  AND  MEASUREMENTS 


HARPOON  MISSILE  FLIGHT  ENVIRONMENTAL 
MEASUREMENT  PROGRAM 

V.  S.  Noonan,  J.  L.  Gubser,  R.  D.  Harmening 
McDonnell  Douglas  Astronautics  Company— East 
St.  Louis,  Mo. 


This  paper  presents  a summary  of  acoustic  and  vibration  measurements  made 
during  the  Design  Phase  flights  of  the  Harpoon  missile.  The  data  is  presented  in 
the  form  of  power  spectral  density  and  shock  spectrum  plots. 


INTRODUCTION 

The  Harpoon  is  an  anti-ship  missile  capable  of 
being  launched  from  various  aircrafts  and  surface 
ships.  The  main  features  of  the  missile  are  shown  in 
Figure  1.  With  the  award  of  the  Harpoon  Missile 
Program  to  McDonnell  Douglas  Astronautics 


Company-East,  an  inflight  dynamics  data  measure- 
ment program  was  implemented  to  evaluate  the 
acoustic  and  vibration  predictions  made  for  the 
missile  and  to  serve  as  a data  base  for  future  tactical 
missile  designs.  This  paper  presents  a summary  of 
the  flight  data  obtained  during  the  Design  Phase  of 
the  program. 


FIGURE  I HARPOON  MISSILE 
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The  objective  of  the  Design  Phase  of  the  Harpoon 
program  was  to  prove  the  adequacy  of  the  overall 
missile  system  concept  and  to  verify  the  compati- 
bility of  the  missile  with  aircraft  and  surface  ship 
launch  platforms.  To  accomplish  these  objectives  re- 
quired several  types  of  vehicles.  The  first  vehicles 
flown  in  the  program  were  Aerodynamic  Test 
Vehicles  ( ATV ) and  Blast  Test  Vehicles  ( BTV ) which 
were  of  boilerplate  construction.  The  ATV’s  were  air 
launched  from  a P-3  aircraft  while  the  BTV's  were 
launched  from  an  ASROC  launcher  located  at  a shore 
site  (Naval  Missile  Center,  Point  Mugu,  California  lor 
on  a DE-1052  ship.  After  successful  boilerplate 
flights  Control  Test  Vehicles  (CTV)  and  Guidance 
Test  Vehicles  (GTV)  were  flown.  These  vehicles  were 
used  to  verify  the  performance  of  individual  and  over- 
all missile  systems.  The  CTV's  and  GTV's  had  no 
warhead  and  were  essentially  identical  with  one 
exception.  The  GTV’s  had  a functional  seeker  while 
the  CTV's  utilized  a dummy  seeker  which  was 
dynamically  similar  to  the  real  seeker.  The  air 
launched  and  surface  launched  CTV’s  and  GTV's 
utilized  a turbo-jet  engine  during  cruise  flight  and  the 
surface  launched  vehicles  used  a solid  propellant 
booster.  CTV's  and  GTV's  were  launched  from  P-3 
aircraft  and  from  an  ASROC  launcher  located  at  a 
shore  site  or  on  a DE-1052  ship. 

Internal  and  external  acoustic  and  vibration 
measurements  were  made  on  ATV’s,  BTV's,  CTV's 
and  GTV's.  Data  presented  in  this  paper  is  in  the 
form  of  power  spectral  density  and  shock  spectrum 
plots.  The  data  provides  a good  cross  section  of 
dynamic  environments  resulting  from  P-3  aircraft 
captive  carry,  solid  propellant  boost  flight,  turbojet 
engine  flight,  high  and  low  altitude  drops,  maximum 
aerodynamic  pressures  and  transient  conditions 
resulting  from  aircraft  launch,  booster  ignition  and 
separation  and  turobjet  engine  ignition. 

DATA  MEASUREMENT  SYSTEMS 

Acoustic  and  vibration  measurements  were  made 
onboard  and  offboard  of  the  Harpoon  missile  using 
piezoelectric  type  sensors.  The  offboard  acoustic  and 
vibration  sensor  systems  had  a 10  to  8000  Hz  and  10 
to  2000  Hz  bandwidth  respectively.  The  onboard 
acoustic  and  vibration  sensor  systems  had  a 10  to 
2000  Hz  bandwidth  which  was  dictated  by  the  missile 
telemetry  system.  Onboard  acoustic  and  vibration 
data  were  measured  continuously  using  IRIG 
channels  15,  17  and  19.  During  the  ATV  captive 
flight  tests,  data  was  “hard  wired"  to  the  P-3.  For 
BTV,  CTV  and  GTV  flights,  the  data  was 
telemetered  to  a ground  station  for  recording. 

SPECTRAL  ANALYSIS  TECHNIQUE 

Data  were  obtained  for  each  flight  and  recorded  on 


magnetic  tape.  The  resulting  time  histories  were  re- 
viewed for  peak  amplitude  values  and  used  to  select 
data  intervals  for  subsequent  PSD  and  shock 
spectrum  analysis.  Upon  review  of  the  time  histories 
the  data  were  classified  as  either  stationary  or  non- 
stationary for  spectral  analyses.  Sationary  data  were 
analyzed  by  the  use  of  power  spectral  density 
analyses  and  transient  data  were  analyzed  by  shock 
spectrum  techniques.  Whenever  possible,  sample 
lengths  and  bandwidths  were  chosen  to  yield 
maximum  frequency  resolution  and  sufficient  degrees 
of  freedom  for  statistical  confidence  in  amplitude 
levels. 

Some  difficulties  were  experienced  in  interpreta- 
tion of  the  power  spectral  density  data  due  to  low 
signal  to  noise  (S/N ) ratios.  The  lowest  possible  max- 
imum amplitude  setting  for  the  vibration  sensor  sys- 
tems was  5.0  g's  peak.  Most  of  the  vibration  levels 
measured  in  the  forward  portion  of  the  missile  were 
less  than  1.0  g peak  for  steady  flight  conditions. 
Therefore  most  of  the  steady  state  data  was  below 
20%  of  full  scale.  The  vibration  measurement 
systems  showed  a nominal  noise  floor  of  approxi- 
mately 7%  of  full  scale.  The  combination  of  the  low 
signal  level  and  the  nominal  noise  level  results  in  a 
S/N  of  3 or  less. 

A S/N  ratio  of  approximately  1 .5  results  in  only 
about  10  dB  of  spectral  resolution.  If  the  portion  of 
the  spectrum  above  300  Hz  is  of  lower  amplitude  than 
the  low  frequency  portion  of  the  spectrum,  this  low 
S/N  will  result  in  masking  the  spectrum  by  the  noise 
floor  leaving  telemetry  cross  talk  hamonics  at  400, 
800,  1200,  1600  and  2000  Hz  as  the  predominate 
peaks.  This  effect  can  be  seen  in  some  of  the  PSD 
plots  presented  later.  When  the  S/N  increased  to 
approximately  2.5  the  spectral  resolution  improved 
to  25  or  30  dB,  but  the  telemetry  noise  at  400  Hz  still 
is  visible.  For  good  resolution  in  power  spectral 
density  analyses,  a S/N  of  at  least  4 is  desirable. 

PSD  analyses  were  performed  on  a Time/Data  100 
Time-Series  analyzer  by  the  Pacific  Missile  Range. 
The  system  uses  a Fast  Fourier  Transform  technique. 
The  analysis  bandwidth  used  on  the  Time/Data  100 
is  a direct  function  of  the  digitizing  rate  used  to  pro- 
cess the  flight  analog  data.  For  the  onboard  data, 
with  a 2KHz  frequency  response,  the  narrowest 
analysis  bandwidth  possible  was  10  Hz.  For  offboard 
data,  with  an  8 KHz  response,  the  narrowest  analysis 
bandwidth  possible  was  20  Hz. 

The  plot  format  utilized  by  the  Pacific  Missile 
Range  for  PSD  analysis  presented  spectral  values  in 
reference  to  a psuedo  peak  unit  squared  per  Hz  rather 
than  a RMS  unit  squared  per  Hz.  The  amplitude 
values  used  were  based  on  a sinusoidal  peak  unit  (i.e., 
1.4  times  the  rms);  thus,  all  PSD  amplitudes  were 
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plotted  a factor  of  two  greater  than  the  usual  format. 
To  prevent  confusion,  the  plots  presented  in  this 
paper  are  copies  of  the  originals  with  the  scales 
changed  to  reflect  the  spectral  values  based  on  rms 
rather  than  peak  values. 

OFFBOARD  MEASUREMENTS 

Measurements  offboard  of  the  missile  were  made 
during  the  BTV,  CTV  and  GTV  flight  tests. 

During  air  launches  of  the  CTV’s  and  GTV's  the 
vibration  environment  inputted  to  the  Harpoon  Data 
Processor  (HDP),  mounted  in  the  P-3,  was  measured. 
The  measurement  was  made  at  the  interface  of  the 
HDP  and  aircraft  structure.  The  HDP  was  mounted 
inside  the  fuselage  structure  away  from  the  missile. 
Thus,  all  vibration  inputs  to  the  HDP  were  P-3 
generated.  Figures  2 and  3 present  PSD  plots  of  the 
longitudinal  and  vertical  (referenced  to  aircraft ) HDP 
vibration  inputs.  A dynamic  analysis  identified  the 
PSD  peaks  at  15,  25,  65,  135,  205  and  275  Hz  as  local 


resonances  of  the  HDP  mounting  structure. 

During  launches  from  the  shore  based  ASROC 
launcher,  acoustic  measurements  were  made  in  the 
live  and  empty  adjacent  cells  and  at  various  ground 
level  locations  in  front  of  the  launcher.  During  ship 
launches  from  the  ASROC  launcher  acoustic 
measurements  were  made  in  the  live  and  empty 
adjacent  cells  and  at  several  deck  locations. 
Locations  of  the  ground  and  deck  sensors  are  shown 
in  Figure  4 and  5. 

A summary  of  the  maximum  overall  sound 
pressure  levels  measured  offboard  during  shore  and 
ship  launches  is  shown  in  Figured  6.  Data  from  5 
shore  launches  and  5 ship  launches  is  represented  in 
the  table.  The  data  measured  in  the  live  cell  during 
shore  launches  was  similar  to  that  obtained  during 
ship  launches.  This  was  also  the  case  for  the  adjacent 
cell  measurements.  The  live  cel)  measurements  are 
somewhat  questionable  due  to  the  heat  generated 
during  missile  booster  ignition  which,  in  effect, 
saturated  the  acoustic  sensor  system. 
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FIGURE  2 POWER  SPECTRAL  DENSITY 
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ON  BOARD  MEASUREMENTS  - ATV‘S 

The  objective  of  the  ATV's  was  to  evaluate  the 
inflight  proximity  effectgs  of  the  Harpoon  and  P-3 
aircraft  and  to  demonstrate  safe  captive  carriage  and 
release  from  the  aircraft.  The  ATV  was  a boilerplate 
structure  which  simulated  the  external  geometry, 
weight,  center-of-gravity,  and  inertias  (pitch  and 
yaw)  of  the  Harpoon  missile. 

One  ATV  was  utilized  for  captive  flight  acoustic 


measurements.  Six  flush  mounted  microphones  were 
located  at  various  points  on  the  missile  shown  in 
Figure  7.  As  the  ATV's  were  boilerplate  no  vibration 
measurements  were  made.  All  data  were  collected 
and  “hard-wired"  to  a frequency  modulated  tape 
recorder  located  in  the  aircraft.  A total  of  six 
A TV/P-3  captive  carriage  flights  were  flown  with  the 
aircraft  performing  various  specified  maneuvers.  In 
general,  the  tests  were  conducted  for  various  speeds 
and  altitudes,  turns,  near  stall,  side  slip  and  dive 
conditions. 


PARAMETER  ACOUSTIC  SENSOR 

DESIGNATION  LOCATION 


VA02 

MS25.5.  0° 

VA03 

MS25.5.  180° 

VA04 

MS25.5,  210° 
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MSI  11.75.  270' 

mi 

MSI  55,  0° 

FIGURE  7 - EXTERNAL  ACOUSTIC  SENSOR  LOCATIONS  - ATV 


A table  of  the  maximum  acoustic  overall  sound 
pressure  levels  for  the  various  maneuvers  and 
locations  is  presented  in  Figure  8.  Review  of  the 
sound  pressure  levels  shows  the  highest  level  to  be 
155  dB  which  was  measured  during  a low  altitude 
dive.  Acoustic  levels  tend  to  be  1 to  5 dB  higher  on 
the  top  (0°)  of  the  missile  as  compared  to  the  same 
station  but  on  the  side  or  bottom  of  the  missile.  This 
effect  was  due  to  turbulence  from  the  missile/rack 
interface.  Also,  the  levels  toward  the  rear  of  the 
vehicle  tend  to  be  higher  due  to  turbulence. 

Review  of  the  CTV  and  GTV  acoustic  data  also 
showed  the  aft  measurements  to  have  a higher  fre- 
quency content  than  the  forward  measurements. 

ON-BOARD  MEASUREMENTS  - BTV'S 

The  objective  of  the  BTV's  was  to  evaluate  booster 
ignition,  performance  and  separation  during  surface 
launched  Harpoon  flights.  Three  BTV’s  were  shore 


launched  from  an  ASROC  launcher  and  one  from  a 
ship  based  ASROC  launcher.  The  BTV’s  were  boiler- 
plate construction  having  the  same  geometry,  mass 
center-of-gravity,  and  inertias  as  the  Harpoon  missile 
but  having  alive  booster. 

The  BTV’s  were  instrumented  to  investigate  the 
external  acoustic  environment  and  the  longitudinal 
vibration  levels  imparted  to  the  missile  by  the 
booster  at  the  booster/missile  interface  during 
booster  operation  and  separation.  The  external 
acoustic  environment  at  stations  M.S.  25.5 
(Parameter  VA10)  and  M.S.  155.4  (Parameter  VA11 ) 
on  the  missile  were  monitored  during  the  flights  of 
two  vehicles.  Both  measurements  were  at  top  vertical 
center.  The  booster/missile  interface  vibration  levels 
were  monitored  by  two  accelerometer  systems;  one 
with  a high  peak  acceleration  level  to  investigate  the 
booster  ignition  and  separation  shocks  and  one  with  a 
low  peak  acceleration  to  investigate  the  booster  burn 
vibration.  Vibration  measurements  were  made  on  two 
BTVs. 
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FIGURE  8 - SUMMARY  OF  MAXIMUM  MISSILE  EXTERNAL  ACOUSTIC  LEVELS  - ATV 


The  maximum  overall  sound  pressure  levels 
measured  during  boost  flight  were  143  dB  at  the  front 
of  the  missile  and  144  dB  at  the  rear.  Figures  9 and  10 
i present  PSD's  of  the  BTV  boost  acoustic  measure- 

ments. Due  to  the  short  duration  of  boost  flight  the 
DEGREES-OF-FREEDOM  are  small  thus  statisti- 
cally compromising  the  amplitude  accuracy  of  the 
data. 

Review  of  booster/missile  interface  vibration  data 
indicated  lower  than  expected  shock  levels  due  to 
booster  ignition  and  sepaiation.  The  high  level 
accelerometer  (full  scale  350g  0-pk)  showed  no  data 
above  the  nominal  noise  floor  (approximately  -g 
0-pk).  The  low  level  accelerometer  (full  scale  20g  0-pk) 
appeared  to  just  saturate  at  both  booster  ignition  and 
separation.  The  levels  during  boost  flight  were  near 
the  nominal  noise  floor  (approximately  2g  0-pk)  and 
were  too  low  to  yield  any  spectral  data.  The  level  of 
the  vibration  data  while  the  missile  was  still  within 
the  launcher  cell  was  approx.  8g  0-pk. 


ON-BOARD  MEASUREMENTS— CTV'S  & GTV'S 

The  Control  Test  Vehicle  (CTV)  and  the  Guidance 
Test  Vehicle  (GTV)  flights  were  conducted  to 
evaluate  the  overall  performance  of  the  prototype 
Harpoon  Weapon  System.  Both  vibration  and 
acoustic  environments  were  measured  throughout 
the  missile  flights  to  investigate  aircraft  captive 
carry,  boost  and  cruise  levels.  A total  of  15  CTV  and 
GTV  flights  were  made,  of  which  6 were  surface 
launches  and  9 were  air  releases.  The  surface  launches 
were  pereformed  from  ASROC  launchers  which  were 
either  ship  or  shore  based.  The  air  released  flights 
were  conducted  from  a P-3  aircraft.  The  CTV  and 
GTV  were  structurally  representative  of  the  produc-  < 

tion  Harpoon  but  with  the  warhead  replaced  by  a test 
and  evaluation  section  which  contained  telemetry 
components.  The  CTV's  and  GTV's  were  equivalent 
except  for  the  absence  of  an  operational  seeker  on  the 
CTV.  The  CTV  was  fitted  with  a dynamic  model  of  | 

the  seeker  to  provide  complete  dynamic  similarity  d 
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between  the  CTV  and  GTV.  The  only  configuration 
difference  between  flights  was  due  to  the  type  of 
launch.  Surface  launch  missiles  were  fitted  with  a 
booster  and  booster  separation  system  which  allowed 
the  missile  to  achieve  sufficient  velocity  for  turbo-jet 
engine  start.  After  booster  separation  there  was  no 
significant  configuration  difference  between  surface 
launch  and  air  release  missiles. 

Measurements  were  made  to  determine  the 
vibration  levels  inputted  to  the  Missile  Guidance 
Unit  (MGU),  seeker,  turbo-jet  engine,  missile/boost- 
er interface  control  actuators  and  forward  aircraft 


hook.  All  of  these  items  are  mounted  directly  fo 
missile  primary  structure.  The  vibration  sensors  were 
located  at  the  interface  of  the  item  and  the  missile 
primary  structure.  Acoustic  levels  were  measured 
both  internal  and  external  to  the  missile.  Only  three 
continuous  channels  were  available  during  each 
flight,  thus  making  it  necessary  to  change  the 
measurements  for  different  flights  to  obtain  data  at 
various  locations  in  3 axes.  Figure  11  shows  sensor 
locations.  Dynamic  data  were  recorded  during  all 
phases  of  the  missile  flight.  A total  32  vibration 
measurements,  4 internal  acoustic  and  3 external 
acoustic  measurements  were  made. 


W - VIBRATION  SENSOR 
VA  - ACOUSTIC  SENSOR 


VAW 


FIGURE  11  - ACOUSTIC  AND  VIBRATION  SENSOR  LOCATION  - CTV  AND  GTV 


PSD  analyses  were  performed  for  selected  time 
intervals.  Except  for  those  intervals  occurring  during 
boost  flight  it  was  possible  to  use  durations  of  suffi- 
cient magnitude  to  yield  PSD's  with  an  adequate 
number  of  degrees-of-freedom.  The  duration  of  boost 
flight  was  too  short  to  accomplish  this  but  PSD 
analyses  were  still  performed  to  yield  frequency 
information. 

Figures  12  and  13  illustrate  the  external  acoustic 
environment  measured  during  P-3  captive  flight.  The 
measured  data  indicates  that  the  external  acoustic 
level  was  138  dB  at  the  forward  missile  station  and 
145.5  dB  at  the  aft  missile  station.  Figure  14  is  the 
interna!  acoustic  level  measured  at  the  aft  missile 
station  during  captive  carry.  A reduction  in  noise 
level  of  approximately  8 dB  is  realized  through  the 
missile  structure. 

The  vibration  spectra  for  the  P-3  captive  carry  are 
shown  in  Figures  15  through  18.  These  plots  are 
typical  for  all  of  the  captive  carry  flights  for  which 


data  was  obtained.  In  general,  the  vibration  is  low 
level  with  the  amplitude  of  the  high  frequencies 
'ncreasing  for  the  sensors  located  in  the  aft  sections 
of  the  missile.  The  data  is  consistent  since  the 
turbulent  boundary  layer  does  injrease  along  the 
missile  and  this  is  the  primary  source  of  the  high  fre- 
quency vibrations.  P-3  propeller  frequencies  of  15  and 
65  Hz  are  also  evident  in  much  of  the  vibration  data. 
The  high  amplitude  spikes  occurring  at  400,  800, 
1200,  1600  and  2000  Hz  are  the  telemetry  harmonics 
previously  discussed. 

Operating  of  the  missile  turbo  jet  engine  had  a 
pronounced  effect  on  the  vibration  levels  throughout 
the  missile.  The  engine  was  started  while  the  missile 
was  still  captive  to  the  aircraft.  Figures  19  through 
21  show  the  captive  vibration  environment  at  the 
seeker,  MGU  and  engine  after  the  engine  has 
achieved  100%  rpm.  The  first  axial  bending  mode 
frequency  at  approximately  209  Hz  is  evident.  The 
engine  operating  frequency  at  approximately  670  Hz 
is  very  evident. 
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Vibration  resulting  from  booster  burning  was 
much  lower  than  expected.  Figures  22  through  24  are 
spectra  of  the  longitudinal  vibration  and  Figures  25 
and  26  are  spectra  of  the  lateral  vibration  during  the 
boost  phase.  All  vibration  for  this  phase  show  less 
than  1 grans  with  the  maximum  level  shown  in  Figure 
24  which  was  a x-axis  sensor  located  on  primary 
missile  structure  at  the  missile/booster  interface. 

Missile  cruise-flight  starts  with  release  from  the 
aircraft  or  separation  from  the  booster.  Once  this 
occurs  the  acoustic  and  vibration  data  becomes 
essentially  stationary.  Any  changes  in  vibration 
levels  during  cruise  can  normally  be  associated  with 
some  change  in  turob-jet  engine  performance. 
External  acoustic  levels  are  very  low  as  shown  in 
Figure  27.  However,  the  internal  acoustic  level  in  the 
aft  section  of  the  missile  is  primarily  influenced  by 
engine  noise  and  is  higher.  This  internal  acoustic 
spectrum  shown  in  Figure  28  indicates  that  an  overall 
level  of  147  dB  can  be  expected  as  a result  of  the 
engine  noise.  Vibration  levels  during  cruise  flight  are 
shown  in  Figures  29  through  40.  The  measured  levels 
show  approximately  1 grans  in  the  forward  section  of 
the  missile  and  3 grans  in  the  aft  section.  The  engine 


frequency  of  approximately  670  Hz  and  its  1st 
harmonic  are  evident  throughout  the  missile  and  can 
usually  be  seen  in  all  the  frequency  spectra  irregard- 
less  of  the  direction  being  measured.  A first  axial  (220 
Hz)  missile  mode  that  developed  during  the  cruise 
flight  can  be  seen  in  Figure  29  measured  at  the 
forward  end  of  the  missile  and  in  Figure  40  measured 
at  the  aft  end  of  the  missile. 

Shock  spectrum  analyses  were  performed  on  some 
of  the  transient  events.  Figures  41  and  42  show  the 
shock  spectra  at  the  seeker  and  turbojet  engine  for 
the  booster  ignition  transient.  Comparing  Figures  41 
and  42  gives  an  indication  of  the  shock  attenuation 
that  can  be  attained  along  the  longitudinal  axis  of  the 
missile.  Figure  43  through  45  ai  > the  shock  spectrum 
for  booster  separation  and  Figures  46  through  48  are 
for  turbojet  engine  start.  These  figures  again  show 
the  shock  response  in  the  forward  and  aft  sections  of 
the  missile.  Figures  43  and  46  show  the  cross-axis 
response  at  the  seeker  to  these  events.  In  general, 
shock  levels  for  the  major  transient  events  are  much 
lower  than  expected  and  may  be  due  to  the  filter 
action  measurement  systems  which  rolled  off  at  24 
dB  above  2 Khz. 
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FIGURE  22  POWER  SPECTRAL  DENSITY 
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FIGURE  29  POWER  SPECTRAL  DENSITY 
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FIGURE  30  POWER  SPECTRAL  DENSITY 
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FIGURE  31  POWER  SPECTRAL  DENSITY 


PARAMETER:  W12  DEGREES  OF  FREEDOM:  122 

LOCATION:  NGU  - * AXIS  EVERT:  CRUSE  FUGHT 

OVERALL  L3S  GRRS 


FIGURE  32  POWER  SPECTRAL  DENSITY 
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FIGURE  37  POWER  SPECTRAL  DENSITY 
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FIGURE  38  POWER  SPECTRAL  DENSITY 
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CONCLUSIONS 


In  general,  good  acoustic  and  vibration  data  was 
obtained  for  all  flights  during  the  Design  Phase  of  the 
Harpoon  Missile  Program.  This  data  will  be  used  to 
update  weapon  system  development  phase  environ- 
mental requirements  and,  if  necessary,  to  support 
design  changes.  The  extensive  dynamic  data 
provides  a good  basis  for  establishing  dynamic- 


environments  for  future  tactical  missile  programs. 
The  Harpoon  missile  vibration  and  acoustic  data 
provides  information  for  turob-ji-t  engine  and  solid 
booster  motor  operation,  high  and  low  altitude  flight, 
maximum  dynamic  pressure  flight  and  P-3  captive 
flight.  The  data  provides  good  coverage  of  all  missile 
sections  for  all  of  the  various  flight  conditions. 
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NARROW  BAND  TIME  HISTORY  ANALYSIS  OF 


TRANSPORT  AIRCRAFT  VIBRAT ION  DATA 


Roger  E.  Thaller  and  Jerome  Pearson 
Air  Force  Flight  Dynamics  Laboratory 
Wright-Patterson  AKB,  Ohio 


The  requirement  to  transport  fragile  aerospace  payloads  in 
large  aircraft  has  created  the  need  to  accurately  define 
the  cargo  deck  vibration  environments  of  such  aircraft  as 
the  C-135,  the  C-141,  and  the  C-5A.  To  determine  the 
probability  of  occurrence  of  damaging  accelerations  during 
cargo  aircraft  operations,  a special  method  was  developed 
for  analyzing  and  displaying  the  data,  which  may  be  non- 
Gaussian.  The  method  was  applied  to  low-frequency  cargo 
deck  vibrations  in  order  to  develop  vibratio..  test 
qualification  curves  for  a large,  flexible  spacecraft. 

The  results  allow  the  probabilities  of  given  levels  of 
acceleration  in  separate  frequency  bands  to  be  computed 
and  provide  a measure  of  how  these  probabilities  may 
differ  from  the  Gaussian  case. 


1NTTR0DUCT  ION 

In  recent  years  the  availability  of 
large  transport  aircraft  has  allowed  the 
transporting  by  air  of  many  large,  fragile 
aerospace  payloads.  Some  of  these  pay- 
loads,  such  as  spacecraft  and  launch 
vehicles  carried  horizontally,  are  sensi- 
tive to  low  levels  ot  vibration  at 
frequencies  lower  than  10  Hertz.  The 
measurement  of  cargo  deck  vibration 
environments  has  a long  history  (REF.  1-4), 
but  most  of  these  studies  have  dealt  with 
the  frequency  range  above  5 or  10  Hz. 
Vibration  data  are  ample  in  the  higher 
frequency  ranges  up  to  2000  Hz  but  are 
nearly  non-existent  below  5 Hz. 

Faced  with  this  lack  of  data  and  the 
need  to  transport  payloads  with  known  high 
responses  to  0.5  to  2 Hz  excitation,  the 
Air  Force  Space  and  Missiles  System 
Organization  (SAMSO)  requested  that  the 
Air  Force  Flight  Dynamics  Laboratory  per- 
form the  necessary  flight  tests  to  measuie 
low-f requency  responses  of  a large  aircraft 
cargo  deck  and  to  provide  a measure  of  the 
risk  of  given  levels  of  acceleration 
occurring  during  operations.  The  magnitude 
of  the  accelerations  of  concern  is  such 
that  if  they  occur  only  once,  the  payload 
will  be  damaged.  It  was  known  that  large 
cargo  aircraft  have  elastic  modes  with 
frequencies  as  low  as  0.73  Hz  (5),  and  it 
was  desired  to  measure  responses  at  least 


this  low  in  frequency. 

FLIGHT  TESTS 

lhe  riight  vibration  tests  weie 
designed  to  provide  a maximum  amount  of 
information  on  the  probability  of 
occurrence  of  different  levels  of 
acceleration  with  a minimum  of  instru- 
mentation, with  no  modification  to  the 
aircraft,  and  with  few  test  flights.  In 
order  to  accomplish  this,  12  low-frequency 
accelerometers  were  mounted  in  pairs  at 
six  locations  on  the  cargo  deck  as  shown 
in  Figure  1.  The  accelerometers  were 
attached  to  aluminum  blocks  mounted  on  the 
permanent  cargo  tie-down  rails  at  the 
s ides  of  the  cargo  deck.  Since  vertical 
accelerations  were  considered  likely  to  be 
the  highest,  all  six  locations  included  a 
vertical  accelerometer.  The  forward  and 
aft  locations  included  lateral  acceler- 
ometers and  the  center  locations  included 
longitudinal  accelei ome ter s . By  combining 
the  signals  from  pairs  of  these  acceler- 
ometers in  real  time,  roll,  pitch,  and  yaw 
angulai  accelerations  were  also  derived. 

It  was  felt  that  this  amount  of  informa- 
tion would  be  sufficient  to  define  the 
vihratory  environment  of  the  cargo  deck. 

The  data  of  interest  were  chiefly  the 
accelerations  below  10  Hz.  Past  studies 
of  cargo  deck  accelerations  were  limited 
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mainly  to  frequencies  above  5 Hz,  and  the 
results  at  the  lower  frequencies  had 
large  errors.  The  data  were  low-pass 
filtered  at  30  Hz  in  order  to  eliminate 
the  higher  frequency  and  higher  g signals 
which  commonly  occur  in  aircraft  flight. 

The  accelerometers  chosen  had  capability 
to  DC;  however,  the  signals  were  high-pass 
filtered  at  0.2  Hz  in  order  to  prevent 
amplification  of  any  DC  off  cot  in  the 
amplifiers.  The  frequency  range  of  the 
data  was  then  0.2  to  30  H2. 

The  12  channels  of  data  were  amplified 
and  recorded  in  flight  by  a compact  instru- 
ment package  shown  in  Figure  2,  which 
contained  the  amplifiers,  power  supply, 
and  a tape  recorder.  The  amplifiers  were 
eqiipped  with  automatic  gain  controls 
which  adjusted  the  amplifier  gain  in  10  dB 
steps  over  a 70  dB  range  to  ensure  that 
the  signals  were  recorded  with  proper 
amplitude.  The  automatic  gain  change  was 
triggered  when  the  signal  level  was  out  of 
the  amplifier  range  for  more  than  0.4 
seconds.  Most  flight  conditions  of 
interest,  including  taxi,  takeoff,  cruise, 
flight  in  turbulence,  and  descent,  could 
thus  be  recorded  without  advance  knowledge 
of  the  vibration  levels.  For  the  landing 
data,  however,  the  amplifiers  were  locked 
into  a precalculated  gain  step  in  order  to 
avoid  clipping  the  response  at  the  instant 
of  runway  impact. 

Tne  test  flights  were  intended  to 
cover  a range  of  variables  that  could 
affect  the  accelerations  experienced  in 
various  flight  conditions.  Although  only 
one  aircraft  was  used,  the  gross  weights 
during  the  flight  tests  ranged  from 
455,000  pounds  to  673,000  pounds,  several 
different  pilots  flew  the  aircraft,  27 
takeoffs  and  landings  were  made  on 
different  runways,  and  four  hours  of 
cruising  time  were  recorded  in  conditions 
ranging  from  smooth  air  to  terrain 
following  at  low  altitude,  which  simulates 
heavy  turbulence. 


DATA  ANALYSIS 

The  data  tapes  were  analyzed  by  the 
procedure  shown  in  Figure  3.  The  initial 
survey  of  the  results  indicated  that  the 
conditions  of  climb,  cruise  in  smooth  air, 
and  oescent  were  not  important  in  terms  of 
the  cargo  deck  accelerations  produced. 

The  conditions  of  landing,  high-speed  taxi, 
and  low  altitude  terrain  avoidance,  which 
simulated  flight  in  heavy  turbulence, 
produced  the  highest  accelerations  and  were 
chosen  for  detailed  analysis. 

Because  it  was  desired  to  know  the 
vibration  environment  as  a function  of 
frequency,  the  accelerometer  data  were 


immediately  divided  into  twenty  one-third 
octave  bands  with  center  frequencies  of 
0.4  Hz  to  31.5  Hz.  In  order  to  use  the 
standard  one-third  octave  band  filters 
which  are  available  with  center  frequen- 
cies of  3.15  Hz  and  higher,  the  recorded 
data  were  played  back  at  eight  times  the 
recording  speed.  The  0.4  Hz  data  were 
thus  filtered  by  the  3.15  Hz  filter,  and 
similarly  for  the  other  bands.  Faeh  of 
these  20  filtered  time  histories  for  each 
accelerometer  was  then  analyzed  statisti- 
cally. Each  time  history  was  sampled  at 
a high  rate  compared  to  the  signal  fre- 
quency and  an  acceleration  histogram  was 
constructed  from  the  absolute  values  of 
acceleration.  To  completely  describe  each 
filtered  time  history,  the  histogram  and 
its  statistical  parameters  such  as  mean, 
variance,  and  peak  value  should  be  plotted. 
Then  the  character  of  the  data  could  be 
assessed  completely.  However,  this  would 
require  20  plots  for  each  of  12  sensors 
for  each  flight  condition  — a prohibi- 
tively large  number. 

The  method  chosen  to  display  these 
data  more  compactly  was  to  plot  just 
three  values  for  each  frequency  band  — 
the  highest  acceleration  observed,  the 
acceleration  level  below  which  45.45'i  of 
the  data  fell,  and  the  level  below  which 
68. 27%  of  the  data  fell.  A sample  of 
data  plotted  in  this  fashion  is  shown  in 
Figure  4 for  a Gaussian  random  signal. 

The  reason  for  choosing  these  particular 
percentage  levels  is  that  they  correspond 
to  2 e and  ; of  a Gaussian  distribution. 

As  a result,  these  two  curves  should  be 
in  the  ratio  of  2:1  for  Gaussian  data. 

The  peak  value  curve  will  be  a function 
of  the  testing  time  in  terms  of  the 
frequency  - that  is,  the  number  of  cycles 
of  data  measured.  Theoretically,  an 
infinitely  long  testing  time  will  give  an 
infinite  peak  value  and  an  exact  2:1  ratio 
between  the  45.45%  and  68.27%  curves.  For 
practical  testing  times  and  for  physically 
realizable  accelerations,  this  ratio  is 
sL.il  veiy  close  to  2:1,  but  the  peak 
value  curve  usually  falls  between  3 o and 
6 a (6) . 

When  this  method  of  analyzing  and 
d; splaying  the  data  is  applied  to  a time 
signal  whose  character  is  unknown,  it  is 
possible  to  see  very  quickly  how  closely 
the  data  meet  this  test  of  being  Gaussian, 
and  to  make  ome  estimate  of  the  probabil- 
ity of  occjrri  ice  of  any  given  level  of 
acceleration.  rhese  acceleration  proba- 
bilities can  be  used  to  determine  the 
probability  of  damage  to  a particular  pay- 
load.  If  the  payload  is  unable  to  with- 
stand fie  environment,  these  results  can 
then  be  used  to  design  isolation  mounts 
tor  the  required  attenuation. 


XI 


TEST  RESULTS 

Some  representative  data  from  the 
flight  tests  during  landing  and  during 
low-level  terrain  following  have  been 
assembled  for  illustration.  Figure  S 
shows  vertical  accelerations  at  the  right 
rear  cargo  deck  for  27  landings.  This 
composite  plot  was  made  of  all  the 
landings  in  the  test  p'ogran,  and  includes 
three  runways,  several  pilots,  and  a range 
of  gross  weights.  The  result  is  a com- 
posite of  the  tvpical  landing  acceleration. 
Figure  6 shows  the  lateral  acceleration  at 
the  left  rear  cargo  deck  and  Figure  7 
the  longitudinal  acceleration  at  the  left 
center  cargo  deck  for  all  the  landings. 
These  results  show  a ratio  of  about  2:1 
between  the  95.45%  and  68.27%  levels,  and 
indicate  that  the  data  are  not  very 
different  from  Gaussian  excitation  for 
most  frequency  bands.  There  appears  to  be 
sinusoidal  excitation  superimposed  at 
discrete  frequencies,  apparently  corre- 
sponding to  airframe  resonances. 
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Figure  o Lateral  Landing  Accelerations 
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Figure  5 Vertical  Landing  Accelerations 
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Figure  7 Longitudinal  Landing 
Accelerations 
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A typical  time  history  of  landing 
accelerations  in  the  three  orthogonal 
directions  is  shown  in  Figure  8.  The 
highest  value  of  acceleration  recorded 
occurred  in  the  10  Hz  frequency  band. 
Angular  acceleration  measurements  showed 
it  to  be  a roll  mode  which  causes  large 
vertical  accelerations  at  the  sides  of 
the  cargo  deck  in  the  gear-fixed  modes  of 
landing  and  taxiing. 

It  was  hoped  that  some  data  would  be 
obtained  from  flight  in  clear  air  turbu- 
lence during  the  test  program,  but  only 
very  light  turbulence  was  recorded  in 
several  hours  of  flight.  In  order  to 
obtain  a measure  of  the  response  accelera- 
tions to  be  expected  during  turbulence, 
data  were  taken  during  low-level  terrain 
avoidance.  In  this  flight  condition,  the 
aircraft  flies  at  a fixed  height  of  1000 
feet  above  the  ground.  A terrain-avoidance 
radar  corrects  the  flight  path  in  response 
to  changes  in  ground  elevation  ahead  of 
the  aircraft.  The  result  is  a high  level 
of  vibration  which  is  similar  to  flight 
in  very  rough  air.  The  responses  recorded 
in  this  flight  condition  were  used  as  an 
indication  of  those  expected  during  flight 
in  severe  turbulence. 

A summary  of  vertical  accelerations 
at  the  right  rear  cargo  deck  over  14 
minutes  of  terrain  following  is  shown  in 
Figure  9;  Figure  10  and  Figure  11  show 
corresponding  values  for  lateral  accelera- 
tions at  the  left  rear  cargo  deck  and 
longitudinal  accelerations  at  the  left 
center  cargo  deck,  respectively.  Since 
all  these  data  were  taken  during  the 
equivalent  of  heavy  turbulence,  these 
curves  should  be  used  in  conjunction  with 
an  estimate  of  the  risk  of  encountering 
heavy  turbulence  during  expected  flight 
profiles  in  deriving  vibration  test 
criteria.  A short  time  history  of  the 
three  orthogonal  accelerations  during 
terrain  following  is  shown  in  Figure  12. 
Measurements  were  made  during  high-speed 
taxiing,  in  addition  to  landing  and 
terrain  following.  These  results  and  the 
results  of  the  derived  angular  accelera- 
tions are  described  completely  in  a limited 
distribution  report  (7). 


CONCLUSIONS 

The  results  of  this  flight  test 
program  and  analysis  are  summarized  by 
the  plots  of  maximum  vibration  amplitude 
versus  frequency  in  Figure  13.  In  this 
figure,  a vibration  test  curve  from 
previous  data  (8,  9)  is  shown  as  a dashed 
line.  The  maximum  accelerations  observed 
in  this  program  have  been  converted  to 
equivalent  double  amplitudes  and  plotted 
as  open  circles.  By  applying  an  appro- 


priate safety  factor  to  account  for  the 
short  testing  times,  a vibration  tes* 
curve  can  be  drawn  using  these  plots. 

Tne  data  were  found  to  be  approxi- 
mately Gaussian  with  sinusoidal  excitation 
superimposed  at  discrete  frequencies  which 
apparently  correspond  to  airframe  reso- 
nances. Assuming  the  data  are  Gaussian, 
the  68.27^?  level  corresponds  to  o , and  the 
probability  of  occuirence  of  higher  values 
of  acceleration  can  be  computed  from  this 
value  and  from  the  testing  times  given  in 
this  paper.  A different  test  curve  could 
then  be  developed  using  a particular 
criterion  for  the  probability  of  failure. 
This  method  could  be  refined  to  apply  to 
specific  locations  on  the  cargo  deck  by 
using  the  complete  data  available  in 
Reference  7. 
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DISCfSSlON 


Mr.  Senn  (Anns’  Test  & Evaluation  C remind): 

You  mentioned  and  your  picture*  showed  that 
your  data  were  taken  on  the  floor  of  the  air- 
craft with  no  load. 

Mr.  Thaller'.  So,  in  fact  the  plane  was  loaded 
and  had  a variety  of  loads.  Our  test  occurred 
during  actual  use  of  the  aircraft.  We  couldn’t 
put  any  accelerometers  along  the  center 
longitudinal  axis  of  the  C5A  cargo  deck  because 
it  was  actually  in  use  at  the  time. 

Voice:  In  line  with  the  last  question  did  you 

find  any  difference  in  the  response  of  the  pay- 
load  deck,  or  the  cargo  deck,  as  a function  of 
the  payload  or  cargo  weight?  Was  it  significant 
enough  to  warrant  setting  up  different  test 
levels  for  different  weight  payloads? 

Mr.  Thaller:  Tills  was  a quick  reaction  for 

SAMSO  to  obtain  this  data  and  to  give  them  what 
they  want.  1 have  just  written  up  a plan  for 
exactly  that  project.  Ue  are  going  to  do  an 
analysis  of  variance  and  other  statistical 
techniques  to  find  out  any  dependence  upon  cargo 
loading. 


Mr.  Van  Ert  (Aerospace  Corporation):  Were 

any  operations  performed  on  this  1/3  octave 
Hand  data  to  put  it  in  the  terms  ot  the  last 
chart  you  showed  which  inferred  a sinusoidal 
type  of  test  node?  Your  data  is  1/3  octave 
band  data  which  infers  there  is  a bandwidth 
and  then  vour  final  slide  shows  double 
amplitude  as  a function  of  frequency  which 
infers  to  many  people  that  a sine  wave  test 
would  ho  appropriate.  Have  you  somehow 
arrived  at  some  process  to  equate  sinusoidal 
and  random  data? 

Mr.  Hinller:  So,  the  actual  data  was  always 

in  "g's"  and  this  was  just  a convenient  way 
of  clotting  it. 

Mr.  Van  Ert:  Then  it  wouldn't  be  convenient 

for  me  to  sinusoidally  sweep  to  those  amplitudes. 
What  I should  do  is  have  a 1/3  octave  band 
data  raised  to  that  double  amplitude  to  that 
V level. 

Mr,  Thaller:  We  have  sent  SAMSO  time  histories 

of  the  most  severe  accelerations  and  they  are 
going  to  use  them  in  a simulation  program. 

Mr.  Van  Ert:  1 would  expect  it  might  possibly 

be  used  in  the  design  of  an  isolation  system 
rather  than  sinusoidal  testing? 

Mr.  Thallert  I’m  not  really  aware  of  what  they 
are  going  to  do  since  I didn't  reallv  get  into 
that  phase  of  it. 
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INTRODUCTION 

It  is  often  difficult  to  separate  product 
damage  due  to  handling  from  damage  due  to  vib- 
ration. Handling  mechanical  damage  usually 
produces  immediate  failure  when  the  product  is 
dropped,  whereas  the  fatigue  failure  caused  by 
vibration  may  be  inmediate  or  delayed.  Success- 
ful design  and  testing  of  a package  requires 
that  both  environments  be  considered.  Packages 
are  routinely  subjected  to  drop  tests  in  which 
the  drop  heights  are  supposedly  representative 
of  the  expected  handling  environment.  The 
ASTM  vibration  test  for  shipping  containers 
(D999)  requires  bouncing  of  the  package  for 
one  hour  on  a constant  displacement  shaker  at 
1.1  G (frequency  and  amplitude  unspecified). 
This  test  is  admittedly  not  representative  of 
the  dynamic  environment;  as  such,  it  is  not 
the  optimal  test  specification. 

An  example  of  vibration  damage  is  a 
laser  that  after  shipping  had  a very  short 
life.  A sine-sweep  test  uncovered  a resonant 
frequency  of  the  long  thin  laser  tube  of  15  Hz, 
which  is  in  the  range  of  truck  vibrations. 
Stiffening  of  the  package  to  reduce  the  vibra- 
tion input  to  the  tube  cured  the  problem.  It 
is  very  unlikely  that  the  laser  manufacturer 
would  subject  his  product  to  a bounce  test  for 
one  hour.  But  he  would  be  willing  to  subject 
it  to  a vibration  test  representative  of  the 
expected  environment  and  thus  cure  the  problem 
before  his  produce  reached  the  customer. 

Naturally,  realistic  testing  cannot  be 
performed  if  the  environment  is  not  known;  it 
is  the  purpose  of  this  work  to  measure  the 
dynamic  environment  of  common  motor  carriers 
for  use  in  package  testing. 

BACKGROUND 

Ostrem  (1)  published  a survey  recently 
(1972)  of  the  transportation  shock  and  vibra- 
input  to  cargo.  This  survey  was  compiled  from 
measurements  iy  other  researchers  and  contains 
data  on  trucks,  railcars,  ships  and  aircraft. 
The  vibrarion  data  is  presented  as  probability 
of  occurrence  of  acceleration  plotted  versus 


frequency.  The  data  on  trucks  comes  from  four 
sources: 

a)  Foley  (2)  reported  on  the  vibration 
environment  of  a 15  ton  radioactive 
fuel  cask  shipped  cross-country  on  a 
flatbed  tractor-trailer. 

b)  Data  (3)  on  a rebuilt  tractor- trai ler 
combination  was  collected  and  showed 
that  rebuilding  substantially  reduces 
the  vibration  levels. 

c)  Foley  (4)  also  reported  on  a flatbed 
truck  subjected  to  rough  roads. 

d)  Schlue  and  Phelps  (5)  reported  the 
vibration  levels  experienced  in  an  air- 
ride  van  carrying  spacecrafts  across 
the  country. 

Foley  (6)  wrote  a description  and  summary 
of  vibration  data  that  is  particularly  oriented 
toward  packaging. 

A main  point  here  is  that  none  of  these 
measurements  have  been  made  in  commercial 
shipper’s  trucks  under  ordinary  conditions. 

MEASUREMENTS  AND  ANALYSIS 

Accelerations  were  measured  on  four  runs 
of  tractor-trailers  operated  by  a local  truck- 
ing firm.  Two  runs  were  used  for  equipment 
checkout.  Run  3 was  a 40  mile  run  primarily 
over  expressway  with  a 45  foot  trailer  carrying 
a mixed  cargo  of  28,000  pounds.  Run  4,  over 
60  miles  of  mostly  two-lane  road  was  with  a 
40  foot  trailer  carrying  6000  pounds  of  foam 
cups.  Run  5 was  a "peddle  run”  around  the  city 
with  a 30  foot  trailer  and  an  initial  load  of 
6200  pounds.  Run  6 was  over  expressway  with 
a 40  foot  trailer  loaded  with  9200  pounds  of 
auto  parts. 

Accelerations  were  measured  in  tho  ver- 
tical direction  on  the  floor  at  the  front, 
middle  and  rear  of  the  trailer.  Piezoresiscive 
accelerometers  with  a frequency  response  of 
from  DC  to  200  Hz  were  mounted  at  the  rear  of 
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the  trailer.  Piezoelectric  uc^eVrometers 
with  a frequency  response  of  2 to  5000  Hz  were 
used  ut  the  front  and  middle  locations,  lateral 
accelerations  were  measured  at  the  floor  and 
ceiling  in  the  rear  location  using  , iezore- 
sistive  accelerometers  identical  to  those  used 
for  the  rear  vertical  acceleration.  It  should 
be  noted  that  these  accelerometers  recorded 
all  inputs  to  the  vehicle  floor  and  sidewalls 
including  those  transients  caused  by  bouncing 
packages.  Such  inputs  were,  however,  modified 
by  the  structural  elements  of  the  floor  and 
sidewall  located  between  the  source  of  the 
disturbance  and  the  accelerometer. 

Acceleration  data  were  recorded  on  two 
INTERMEO  four-track  cassette  tape  recorders 
witn  one  track  on  eacn  recorder  used  fo>"  voice 
comnentary  and  synchroni zation.  Piezoresistive 
and  piezoelectric  accelerometers  were  used. 

The  frequency  response  of  ihe  recording  syster 
was  DC  to  3db-down  at  400  Hz.  A separate 
reproducer  is  required  with  this  recording 
system. 

Standard  kinds  of  analyses  were  made. 

An  IBM  1800  computer  was  used  to  compute 
probabilities  and  a SAICOR  real  time  analyzer 
was  used  for  tne  buK  of  the  analysis. 

RESULTS  AND  DISCUSSIONS 

The  results  presented  and  the  discussions 
of  them  are  pointed  toward  the  purpose  of  this 
research;  namely,  the  measurement  of  the  dy- 
namic environment  for  i se  in  package  testing. 

Figure  1 is  a plot  of  the  rear  vertical 
acceleration  over  the  entire  Run  3.  It  shows 
that  the  distribution  is  sharper  than  would  be 
recorded  for  Gaussian  randomness;  the  moment 
of  kurtosis  is  4.3.  The  standard  deviation  is 
0.46  G,  and  there  is  some  departure  from  sym- 
metry for  the  higher  accelerations.  Although 
not  Gaussian,  the  data  are  similar  .ough  to 
justify  using  the  simpler  Gaussian  approach  to 
planning  and  analysis  of  tests. 

A VI BRAN  plot  of  probability  of  acceler- 
ation versus  frequency  was  obtained  for  a 
portion  of  Run  3 at  55  mph  using  a filter  and 
the  probability  program--see  Figure  2.  The  same 
bandwidth  and  center  frequencies  as  in  Ostrem's 
report  (1)  were  used.  Comparison  with  the  data 
in  that  report  shows  generally  lower  accelera- 
tions were  measured  in  commercial  trucks. 

The  basic  character  of  the  truck  vibra- 
tion is  shown  in  the  PSD  plot  of  Figure  3. 

There  are  several  fairly  sharp  peaks  indicating 
harmonics  in  the  0-50  Hz  range,  a sharp  peak 
at  approximately  70  Hz,  and  narrow-band  random 
vibration  in  the  120-200  Hz  range.  This 
general  picture  appears  for  all  the  runs  with 
variation  of  the  levels  and  frequencies. 

The  effect  of  speed  is  seen  in  Figure  3 
to  be  slight  elevation  of  the  PSD  levels  with 


little  change  in  shape  of  the  plot.  RMS 
values  ( 2-400  Hz)  were  measured  to  be  0.52  G 
at  45  mph,  0.54  G at  50  mph,  0.55  G at  55  mph, 
and  0.63  G »f  f>0  mph.  RMS  and  PSD  measurements 
were  made  on  data  for  steady  speeds  over  con- 
crete expressway.  There  is  no  significant 
change  with  speed,  i.e.  there  are  no  new 
frequencies  introduced  as  speed  is  increased. 

It  seems  reasonable  then  to  take  a PS0  plot  of 
a high-speed  run  as  representative  of  the  most 
severe  dynamic  environment. 

Figure  4 shows  that  the  vertical  acceler- 
ation is  largest  at  the  rear  position  in  the 
trailer  for  Run  3.  Furthermore,  the  lower 
frequency  harmonics  are  present  only  on  the 
rear  signal  which  implies  that  they  are 
generated  by  the  tires  and  suspension  system. 
This  was  also  true  for  Run  6,  but  the  RMS 
values  of  the  vertical  accelerations  of  Run  4 
were  nearly  equal  for  all  three  positions  in 
the  trailer.  Other  results  (3)  have  shown 
that  the  forward  position  ha'  higher  vertical 
acceleration.  For  realistic  vibration  data, 
one  should  record  vertical  accelerations  in  at 
least  the  fore  and  aft  positions  of  the  truck. 

Lateral  accelerations  were  found  to  be 
much  smaller  than  the  vertical  as  evidenced  by 
Figure  5.  This  is  plotted  on  a smaller  fre- 
quency scale  because  the  lateral  data  contained 
little  information  over  100  Hz.  Note  that  the 
shapes  of  t>.e  PSD  plots  are  quite  similar 
demonstrating  the  strong  influence  of  the  vib- 
ration input  through  the  rear  suspension.  One 
need  only  record  the  vertical  accelerations  to 
document  the  most  severe  vibrations. 

Figure  6 compares  the  rear  vertical  vib- 
ration for  the  three  over-the-road  runs.  The 
general  shape  and  levels  are  roughly  similar. 
This  is  encouraging  because  it  leads  one  to 
think  that,  in  spite  of  the  tremendous  number 
of  variables  in  trucks,  cargo,  roads,  etc., 
there  may  be  sufficient  similarity  in  the 
vertical  acceleration  PSDs  to  make  a composite 
PSD  useful . 

It  may  be  premature  with  the  limited 
amount  of  data  we  have,  but  a composite  PSD 
for  high-speed  trailer  vibrations  based  on  our 
results  would  appear  as  in  Figure  7.  Possibly 
the  inclusion  of  more  data  would  fill  in  the 
region  between  30  and  65  Hz.  Foley  (6)  has 
compiled  a composite  derived  from  data  on  7 
truck  types  with  5 load  ranges.  It  is  very 
interesting  to  note  that  his  low  frequency  PSD 

level  is  exactly  the  same  (1.5  x 10’^  G^/Hz) 
as  ours.  This  fact  lends  support  to  the 
argument  that  a composite  can  be  developed  that 
is  representative  of  many  truck  types  anu 
conditions. 

A truck,  in  the  course  of  its  run,  will 
experience  several  transients  from  inter- 
sections, railroad  crossings,  pavement  cracks, 
etc.  These  show  up  on  the  floor  of  the  truck 
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as  an  Increase  in  the  level  of  low  frequency  REFERENCES 

components,  not  as  a sharp  pulse.  In  other 


words,  the  acceleration  signal  as  seen  on  a 1. 

strip-chart  recorder  builds  up  as  a nearly 

pure  sinusiod  at  low  frequency  which  then  dies 

out.  Figure  8 is  a plot  of  wave  analyses  taken 

in  the  "capture"  mode  with  the  real  time 

analyzer.  Each  analysis  is  for  one  second  of  2. 

data.  The  truck  drove  onto  a bridge  at 

60  mph.  The  wave  analysis  taken  while  the 

trailer  is  “ringing"  in  response  to  the 

transient  input  shows  the  greatest  increase  in 

an  acceleration  level  at  very  low  frequencies.  3. 

It  does  not  appear  necessary  to  test  a package 

separately  for  vibration  transients;  the  drop 

test  which  is  much  more  severe  is  adequate. 

In  fact,  the  largest  peak  G for  a transient  for  4. 
any  run  was  5 G;  most  peaks  were  around  2 G. 

An  interesting  characteristic  of  trailer 
vibration  can  be  seen  in  Figure  9.  It  presents 
wave  analyses  of  the  “peddle  run"  truck  as  It  5. 

started  and  accelerated  to  40  mph.  The 
basic  shape  of  the  frequency  distribution  is 
established  as  soon  as  the  truck  begins  to 
move;  increasing  the  speed  merely  raises  the 
level  of  vibrations.  6. 

CLOSING 

These  preliminary  results  show  that  the 
data  can  be  acquired  and  managed  Into  a useful 
form  without  too  much  difficulty.  These 
results  plus  those  of  Foley  (6)  disclose  that 
the  vertical  accelerations  are  most  severe 
and  therefore  are  the  only  ones  needed.  This 
Is  Important  If  one  Is  contemplating  recording 
enough  sample  vibrations  generate  a 
statistically  reliable  picture. 

The  package  test  one  would  aim  for  would 
be  the  kind  where  an  electrohydraulic  shaker 
was  programmed  to  reproduce  truck  floor  vib- 
rations at  a large  table.  If  the  goods  were 
shipped  In  stacks  of  five,  then  a stack  of  five 
would  have  to  be  tested  to  examine  the  effect 
of  stacking.  If  the  programmed  vibrations 
were  representative,  this  should  be  a good 
test. 

Even  if  one  has  adequate  data  about 
truck  floor  vibration  in  hand,  the  question  of 
how  to  c mpile  this  into  a simple  usable  form 
is  a difficult  one.  Simply  enveloping  the  PSD 
data  is  selecting  the  worst  possible  case  and 
may  be  overly  severe.  It  is  likely  that  some 
type  of  testing  at  various  PSD  levels  for 
specified  periods  of  time  would  be  better. 
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Figure  1.  Probability  density 
plot  of  rear  vertical 
acceleration  of  Run  3. 
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Figure  3.  Effect  of  speed  on 
the  rear  vertical 
acceleration  of  Run  3. 
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Figure  4.  Effect  of  accelerometer 
location  on  the  vertical 
acceleration  of  Run  3. 
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Sr.  Storav  (Army  Ml»«lle  Command):  Do  you  heve 

eny  plena  In  the  future  to  see  la  up  the  levala 
end  scala  down  tha  tine  for  leboretory  tasting! 

Hr.  Sharpe;  No,  we  eren't  that  far  along  yet. 

In  terns  of  the  problem  of  packaging  It  Is 
necessery  to  heve  e suitable  fellure  criterion, 
depending  on  the  product,  think  thet  type 
of  tine  scaling  would  have  to  go  elong  with  the 
particular  type  of  product. 

Mr.  Certel  (Kinetic  Svstens):  I was  curious  If 

you  hed  considered  trying  e s toiler  aeries 
of  meesuranents  on  a peckage  tester  to  sea  how 
this  conpares  with  some  of  your  actual  meesured 
observations? 

Mr.  Sheroe:  Whet  package  tester? 

Mr.  Gertel:  The  bounce  tester,  so  celled  or 

thet  family  of  testing  devices. 

Mr.  Sharpe:  Are  you  referring  to  the  ASTM  type 

test? 

Mr.  Gertel:  Yes. 

Mr.  Sharpe:  We  did  a quick  and  preliminary  look 

at  that  and  that  generates  a spectrum  In  the 
package  that  consists  of  a series  of  sherp 
lines,  dependent  upon  particular  typa  of 
package.  We  haven't  gone  into  a thorough 
Investigation  of  that  but  it  is  a completely 
different  type  of  spectra  that  you  obtain. 


Mr,  Condourls  (Army  Electronics  Command): 

Did  you  examine  any  of  tha  packages  that  wera 
transports!  to  aee  If  i-hay  wars  damaged.  If 
they  had  been  damaged  wera  thay  proparly 
packagad? 

Mr.  Sharpe:  No. 

Mr.  Condourls:  We  ran  some  tests  on  hallcoptars 

with  Army  electronic  gear  packagad  in  their 
proper  container,  either  woodan  cases  or  card- 
board boxes,  and  we  found  that  very  llttla 
vibration  which  came  through  from  the  helicopter 
deck  get  Into  the  actual  electronics;  If  the 
package  were  properly  protected  we  felt  that 
there  would  be  no  need  for  a vibration  test, 
so  I was  curious  If  you  had  data  of  thla  kind? 

Mr.  Sharpe:  No,  the  thrust  of  this  work  is  to 

obtain  anvlronmental  data  so  we  could  tast 
packaged  cargo  In  the  laboratory. 
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THE  DYNAMIC  ENVIRONMENT  OF  LANDING  CRAFT* 


Mark  B.  Gens 
Sandla  Laboratories 
Albuquerque,  New  Mexico  87115 


Measurements  were  made  of  the  dynamic  environment  found  on 
the  cargo  decks  of  several  amphibious  vehicles  and  assault 
boats  employed  by  the  Navy  and  Marine  Corps  for  assault 
landings.  The  work  was  done  In  cooperation  with  the  Naval 
Weapons  Evaluation  Facility,  Klrtland  Air  Force  Base,  New 
Mexico,  the  Naval  Amphibious  Base,  Coronado,  San  Diego, 
California  and  the  Marine  Corps  Landing  Force  Development 
Center,  Camp  Pendleton,  California.  The  vessels  measured 
were  LVTP7 , LCA2,  LCM8  and  LCU. 

The  reduced  data  show  relatively  low  magnitudes,  but  high 
Incidence  of  discrete,  Intermittent  and  narrow-band  random 
excitation.  The  LVTP7  exhibited  the  highest  amplitudes 
both  during  land  and  water  operation. 


INTRODUCTION 

Supplies  for  ground  troops  during 
an  Invasion  from  the  sea  may  be  brought 
In  by  helicopter  or  by  small  craft  from 
the  fleet  standing  off  the  coast.  The 
dynamic  environment  encountered  during 
such  transfer  Is  a part  of  the  overall 
environment  to  which  tactical  weapons 
may  be  exposed.  The  helicopter  envi- 
ronment has  been  reported  previously 
[1,2,3].  This  document  is  a report  on 
the  landing  craft  environment. 

A study  of  the  landing  craft  envi- 
ronment was  carried  out  with  the  coop- 
eration of  the  Naval  Ordnance  System 
Command  (NAVORDSYSCOM) , Waehlngton, 

D.  C.,  and  the  Naval  Weapons  Evaluation 
Facility  (NWEF) , Klrtland  Air  Foret 
Base,  New  Mexico. 

MEASUREMENT  OF  THE  DYNAMIC  ENVIRONMENT 
General 

Measurements  were  made  on  two 
types  of  vessels.  The  LVTP  and  LCA  are 
amphibious  craft.  Both  the  LCM  and  LCU, 
although  very  different  In  configura- 
tion and  in  size,  are  capable  of 
beaching. 

♦This  work  was  supported  by  the  U.  S. 

Atomic  Energy  Commission. 


Instruments 

Measurements  were  made  by  means  of 
the  ELI  31  recorder  [4]  with  piezo- 
electric and  plezoresi stive  accelerom- 
eters. The  instrumentation  (Table  I) 
was  configured  in  advance  becaur-e  of 
the  necessity  of  rapid  setup  and  take- 
down to  avoid  interference  with  the 
assault  exercise.  Even  though  re- 
stricted by  the  length  of  the  available 
cables,  accelerometer  placement  at  the 
forward  and  aft  areas  of  the  cargo  deck 


TABLE  I 

Instrumentation  for  Landing  Craft 


Channel 

Instrument 

Type 

Axis 

Calibration 

(8) 

1 

PE 

L 

2.0 

2 

PE 

T 

2.0 

3 

PE 

V 

2.0 

4 

PR 

L 

10.0 

5 

PR 

T 

10.0 

6 

PR 

V 

10. 0 

7 

PE 

L 

20.0 

8 

PE 

T 

20.0 

9 

PE 

V 

20.0 

10 

PR 

L 

2.0 

11 

PR 

T 

2.0 

12 

PR 

V 

2.0 
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was  possible  in  all  cases  except  on  the 
LCD. 

One  limitation  became  apparent  when 
data  were  reduced.  The  recorder  had 
developed  a malfunction  which  appeared 
in  the  data  as  a high  number  of  spurious 
peaks  in  the  frequency  band  between  180 
and  240  Hz.  Wherever  this  spurious 
data  completely  masked  real  data,  the 
values  in  that  frequency  band  were 
deleted . 

Data  Reduction 

The  data  were  reduced  in  the  fol- 
lowing steps: 

1.  A real-time  oscillograph  record 

2.  Amplitude  spectral  density  plots 

3.  Vibran  records 

4.  VAIL  summaries  [5] 

Examination  of  the  oscillograph 
recorus  revealed  that  the  amphibious 
vehicles  were  dominant  in  amplitude  over 
a'l  of  the  landing  craft;  therefore,  a 
decision  was  made  to  reduce  a complete 
cross  section  of  data  for  these  vehicles 
and  only  typical  events  for  the  others. 

DYNAMIC  ENVIRONMENT  OF  LANDING  CRAFT 

The  dynamic  environment  of  trans- 
port has  been  described  as  complex  [6]. 
It  consists  of  broadband  random  excita- 
tion distributed  in  a Gaussian  manner, 
as  well  as  some  discrete  elements  which 
may  be  either  superimposed  upon  the 
broadband  base  or  intermixed  with  it. 

All  of  these  elements  are  found  in  the 
dynamic  environment  of  landing  craft. 

The  dynamic  environment  of  landing 
craft  may  be  characterized  as  basically 
broadband  Gaussian  in  nature,  with,  in 
many  cases,  discrete  components  super- 
imposed on  it . When  the  discrete  com- 
ponents are  independent  of  the  basic 
continuous  excitation,  they  are  found 
in  the  top  1 percent  of  the  peaks. 

These  independent  peaks  may  occupy  as 
much  as  80  percent  or  as  little  as  1 
percent  of  the  amplitude  range.  At 
times  discrete  components  are  intermixed 
with  the  continuous  components  to  the 
extent  that  they  may  not  be  identified 
separately  even  though  their  presence 
is  detectable.  Although  this  discussion 
is  devoted  largely  to  the  basic  broad- 
band excitation,  the  presence  of  dis- 
crete elements  (as  shown  on  arrowheads 
on  the  appropriate  frequency  band)  are 
pointed  out  from  time  to  time. 

Examination  of  the  VIBRAN  records 
for  each  event  and  many  of  the  VAIL 
summaries  revealed  major  differences 


between  the  forward  and  aft  positions. 

On  the  amphibious  vehicles  (LVTP7  and 
LCA2)  the  records  for  the  forward  ac- 
celerometers were  generally  Gaussian 
random  in  distribution  while  those  from 
the  aft  position  were  dominated  by  dis- 
crete, often  intermittent,  excitation. 
The  records  for  the  assault  boats  (LCM8 
and  LCU)  were  reversed.  The  aft  records 
were  largely  Gaussian  random  while  the 
forward  records  showed  a preponderance 
of  discrete  excitation.  In  all  cases, 
the  Gaussian  random  excitation  was 
detected  near  the  engines.  This  sug- 
gests that  the  vehicle  structure  tends 
to  damp  out  the  steady  state  excitation 
and  that  the  discrete  measurements  were 
the  result  of  external  excitation  from 
the  tracks,  bow  slap,  etc. 

DYNAMIC  ENVIRONMENT  OF  AMPHIBIOUS 
VEHICLES 

Two  types  of  amphibious  vehicles 
were  encountered.  One,  the  LVTP7,  is 
representative  of  a group  of  vehicles 
which  may  be  characterized  as  floating 
armored  personnel  carriers . These  ve- 
hicles are  propelled  on  land  by  tracks 
and  in  water  by  hydro jet.  The  other 
amphibious  vehicle,  the  LCA2,  is  pri- 
marily a boat  which  has  the  capability 
of  moving  on  land  by  use  of  tracks  and 
in  water  propulsion  by  hydro jet.  Little 
data  are  available  from  the  LCA2  since 
it  broke  a track  during  the  approach  to 
the  ocean.  Because  the  data  obtained 
indicated  values  lower  than  those  of 
the  other  amphibious  vehicle  and  because 
the  LCA2  is  not  in  general  use,  it  was 
decided  not  to  pursue  the  study  of  this 
vehicle  further. 

LVTP7  (Figure  1) 

Transducer  placement  in  the  LVTP7 
included  triaxial  piezoelectric  accel- 
erometers near  the  forward  bulkhead  and 
aft,  near  the  ramp  hinge,  as  well  as  a 
triaxial  cluster  mounted  approximately 
amidships . 

The  LVTP7  exhibits  two  levels  of 
vibration.  The  amplitudes  measured 
during  progression  on  land  fall  in  a 
higher  order  of  magnitude  than  those 
measured  in  the  water. 

Measurement  began  as  the  LVTP7 
moved  across  a concrete  area.  The  ve- 
hicle continued  down  a concrete  ramp 
and  entered  a calm,  land-locked  bay. 
After  crossing  the  bay,  the  vehicle 
ascended  a sandy  beach,  crossed  an  area 
of  loose  sand,  descended  to  the  shore- 
line, and  entered  the  surf  of  the 
Pacific  Ocean.  It  continued  out  through 
the  surf,  turned,  and  returned  to  its 
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Fig.  1 - LVTP7  Amphibious  vehicle 


base  by  following  the  same  route  in  re- 
verse. The  final  record  was  made  as 
the  vehicle  turned  on  the  concrete 
pavement  before  shutting  down. 

The  data  for  the  LVTP7  were  reduced 
in  amplitude  spectral  density  format  by 
use  of  program  MAC RAN  and  in  a peak 
amplitude  distribution  form  by  use  of 
program  VIBRAN.  The  VIBRAN  records  were 
combined  by  program  VAIL  to  give  sum- 
maries of  the  land  events  and  the  water 
events  separately,  as  well  as  a com- 
posite of  all  events. 

The  broadband  base  vibration  re- 
sults for  the  land  events  are  summarized 
in  Figure  2.  Levels  shown  are  higher 
than  any  associated  with  other  types  of 
landing  craft  or  with  the  LVTP7  during 
water  events.  Axis  domination  in  this 
vehicle  is  distributed  typically.  The 
lower  frequencies,  up  to  180  Hz,  are 
dominated  by  the  vertical  axis.  Except 
in  the  300-  to  1000-Hz  banus  ano  above 
1400  Hz.  the  frequencies  above  240  Hz 
are  generally  governed  by  the  transverse 
axis.  The  highest  amplitudes  are  found 
in  the  500-  to  700-Hz  frequency  band 
and  in  the  longitudinal  axis. 

Comparison  of  the  land  event  data 
obtained  from  the  forward  and  the  aft 
accelerometer  groups  reveals  differences 
in  the  amplitudes  and  axis  dominance. 
Figures  3a  and  3b  are  derived  from  the 
VAIL  summaries  for  each  position.  In 


the  forward  position  (Figure  3a)  the 
vertical  axis  dominates  up  to  1000  Hz; 
at  higher  frequencies  it  drops  below 
both  the  longitudinal  and  transverse 
axis.  At  the  aft  end  of  the  deck 
(Figure  3b),  however,  the  vertical  axis 
is  dominant  only  to  350  Hz  and  then 
again  above  1000  Hz.  In  the  forward 
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Fig.  2 - LVTP7  Land  Events,  Continuous 
Excitation  Summary 
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end  Che  highest  amplitudes  are  In  the 
vertical  between  240  and  1000  Hz;  In 
the  aft  end  the  longitudinal  axis  shows 
the  highest  amplitudes  In  a relatively 
narrow  band  between  500  and  700  Hz.  In 
the  same  band,  the  transverse  axis  shows 
higher  amplitude  than  Is  shown  on  any 
of  the  vertical  records. 

Comparisons  between  the  forward 
and  aft  records  by  axis  are  shown  In 
Figures  4a,  4b,  and  4c.  Generally,  the 
highest  amplitudes  are  found  at  the  aft. 
In  the  longitudinal  and  transverse  axes, 
except  for  the  higher  frequencies  (>700 
Hz),  the  higher  amplitudes  are  shown  In 
the  aft  records.  In  the  vertical  axis, 
although  the  high  amplitudes  occur  at 
the  aft  In  the  middle  frequencies,  they 
are  almost  always  shown  In  the  forward 
records  at  the  high  and  low  frequencies. 
The  exception  Is  between  1400  and  1900 
Hz,  where  the  aft  record  shows  high 
amplitudes . 

The  general  characteristics  of  the 
dynamic  environment  of  the  LVTP7  during 
water  events  are  much  the  same  as  those 
during  the  land  events.  Figure  5 is  a 
summary  of  the  water  events.  The  out- 
standing difference  between  the  two 
types  or  events  Is  the  amplitude  of 
overall  dynamic  excitation.  In  all 
instances,  the  data  recorded  In  the 
water  are  lower  in  amplitude  than  those 
recorded  on  land.  Again,  the  excita- 
tion on  the  forward  deck  generally  is 


lower  In  amplitude  than  that  of  the 
after  deck  (see  Figures  6a  and  6b) . 

Figures  7a,  7b,  and  7c  show  the 
relationship  between  tne  forward  and 
aft  records  by  axis.  Once  again  there 
is  a striking  similarity  In  the  shapes 
and  In  the  portion  of  the  vehicle  which 
gives  the  highest  amplitudes.  There 
are  some  minor  differences;  however, 

In  those  frequency  bands  where  the  two 
records  reverse  dominance,  there  is 
little  difference  In  amplitude. 

Inspection  of  the  VAIL  summaries 
Indicates  the  presence  of  a discrete 
component  mingled  with  the  broadband 
Gaussian  random  excitation  In  the  500- 
to  700-Hz  frequency  band.  This  phenom- 
enon Is  particularly  apparent  In  the 
aft  position  In  records  of  the  land 
events . The  plots  of  the  power  spectral 
density  reveal  the  presence  of  a "spike" 
In  the  500-  to  700-Hz  area  In  almost 
every  Instance,  both  forward  and  aft 
and  on  land  or  In  water.  This  confirms 
the  presence  of  a discrete  element 
which  seems  to  be  a narrow-band,  ran- 
domly distributed  excitation.  Since 
it  is  not  possible  to  separate  the 
discrete  component  from  the  broadband 
random  background,  its  presence  has 
been  shown  by  an  arrowhead  on  the  graph. 

In  summary,  the  dynamic  environment 
measured  on  the  LVTP7  showed  higher  am- 
plitude than  those  of  any  other  vessel. 
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Th>_  land  events  gave  higher  readings 
than  the  water  events,  and  the  aft  posi- 
tion showed  higher  amplitudes  than  the 
forward  position.  The  highest  reading, 
2.7g,  was  recorded  in  the  500-  to  700 -Hz 
band  in  the  aft  position  during  land 
events.  The  complex  character  of  the 
excitation  f«  more  evident  in  records 
from  the  water  events  than  from  the 
land  events. 

LCA2 


Measurement  of  the  dynamic  environ- 
ment of  the  LCA2  was  accomplished  by  use 
of  the  same  Instrumentation  as  for  the 
LVTP7.  The  triaxial  clusters  of  piezo- 
electric accelerometers  were  located 
aft  in  the  vicinity  of  the  ramp  hinge 
and  forward  near  the  bulkhead . The 
piezoresistive  instruments  were  placed 
amidship  on  the  plate  directly  above 
the  engines. 

At  the  time  of  boarding,  the  ves- 
sel was  located  on  the  beach,  parallel 
to  the  surf  line  and  above  the  coastal 
ridge.  The  ridge  was  about  5 feet  above 
the  high  tide  level.  The  plan  was  to 
drive  along  the  ridge  until  reaching  a 
convenient  place  to  descend  to  the 
beach,  cross  the  beach,  enter  the  surf, 
and  run  out  several  hundred  yards  into 
the  ocean.  The  route  would  then  be 
reversed,  with  the  vessel  returning  to 
the  same  position.  The  only  portion  of 
the  program  which  was  accomplished  was 
raising  the  ramp,  backing  to  position, 
and  starting  toward  the  point  of  de- 
scent. The  track  broke  at  this  point; 
therefore,  measurement  of  the  entire 
trip  uad  to  be  abandoned. 

The  data  shown  in  Figure  8 were 
derived  from  the  three  records  avail- 
able. The  vertical  axis  is  generally 
dominant  below  350  Hz,  except  in  the 
20-  to  40-Hz  band,  where  the  transverse 
axis  is  slightly  higher.  In  the  fre- 
quencies higher  than  350  Hz  the  trans- 
verse and  longitudinal  axes  essentially 
alternate  in  dominance.  The  presence 
of  discrete  components  is  apparent  only 
in  the  transverse  axis  in  the  700-  to 
1000-Hz  band . This  appears  to  be  a 
contribution  from  the  forward  set  of 
accclcrcmctcrs.  Figures  9a  and  9b  show 
the  forward  and  aft  measurements.  Com- 
parison of  the  two  plots  shows  that,  in 
general,  the  higher  amplitudes  are  found 
in  the  records  taken  from  the  aft  deck. 
This  is  especially  true  below  700  Hz. 
Comparisons  by  axis  (Figures  10a,  10b, 
10c)  show  that  only  in  the  transverse 
axis  is  the  forward  record  frequently 
higher  than  the  aft,  and  then  in  only 
one-third  of  the  frequency  bands. 


The  presence  of  discrete  components 
may  be  detected  by  inspection  of  the 
VAIL  summaries.  Even  though  it  seems, 
in  the  overall  summaries,  that  narrow- 
band  random  excitation  is  present  in 
many  frequency  bands  in  all  axes,  it 


K.TU: 

1.  IlMGITUMMI.  AM)  VERTK.  At  AXES  HAVE 


LANDING  CRAFT  - UA2 
UN)  EVENT*  IK  U Vll 
All  POSmCMM 


I*  100 


APPROXIMATE!  T EQUAI  AMPLITUDES. 
DSC  RET  E COMPONENT  IN  TM  VENTRAL 
AXIS  ONI  T. 

TRANSVERSE  AND  VEATK  Al  AXES  RAVE 
APPRUKIMAtl  I Y EQUAL  AMPLITUDES. 
(JSC RITE  COMPONENT  IN  THE  VERTICAL 
AXIS  OKI  Y. 


1000 


» HEQt'l  N<  Y (Hr) 


Fig.  8 - LCA2  Land  Events,  Continuous 
Excitation  Summary 
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Fig.  9 - LCA2  Land  Events,  Excitation 
Forward  and  Aft 
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may  be  seen  that  this  effect  is  found 
largely  in  the  record  of  the  forward 
accelerometers.  Little  evidence  of 
intermixed  discrete  excitation  is  shown 
in  the  records  from  the  aft  accelerom- 
eters. 

DYNAMIC  ENVIRONMENT  OF  ASSAULT  BOATS 

Assault  boats  are  small  vessels 
which  may  be  used  to  carry  troops, 
vehicles,  and  supplies  from  the  fleet 
to  the  bea<~h.  They  are  characterized 
by  a shallow  draft,  a bow  ramp,  and  an 
ability  both  to  run  up  on  the  beach  and 
to  back  off  from  it.  The  types  in  cur- 
rent use  are  landing  craft  mechanized 
(LCM6  and  LCM8)  and  landing  craft 
utility  (LCU). 

The  LCM6  and  LCM8,  called  "Mike” 
bo  its,  differ  in  minor  details,  but  are 
us,.d  in  the  same  manner.  The  engines 
and  command  facilities  are  located  in 
the  aft  section.  Figures  11a  and  lib 
are  views  looking  forward  to  the  ramp 


and  aft  to  the  quarter  deck.  The  LCl', 
which  is  larger  than  either  of  the  Mike 
boats,  is  configured  with  the  engines 
below  aft  and  the  command  facilities 
contained  in  a superstructure  located 
at  the  aft  portion  of  the  starboard 
side . 

LCM8 

The  measurements  were  made  on  the 
Mike  8 during  a run  for  the  beach  that 
was  part  of  an  assault  exercise  being 
conducted  by  COMPHIBFAC.  The  boat  was 
aboard  the  LPD  (USS  Dubuque)  at  the 
time  of  instrumentation.  Because  the 
cargo  (trucks)  bad  not  been  loaded  tfien 
the  boat  was  available  to  Install  the 
Instruments,  the  accelerometers  were 
placed  >.n  the  gunwale  rather  than  on 
the  rargo  floor.  Figure  12  is  a dia- 
gram f the  final  configuration. 

The  measurements  were  made  during 
the  trip  from  the  dock  ship  to  the 
beach.  The  actions  Included  a period 


Fig.  11  - LCM8  Assault  Boat 


Fig.  12  - Test  Configuration  of  the  LCM8  (no  scale) 
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Fig.  13  - LCM8,  Continuous 

Excitation  Summary 


of  circling  with  the  other  boats  until 
the  signal  was  given  to  run  straight 
for  the  beach.  Measurements  were  made 
during  the  pro-run  maneuvers,  during 
the  run  for  the  beach,  and  when  the 
beach  was  reached. 

When  the  data  were  reduced,  it  be- 
came apparent  that  the  records  taken 
from  the  forward  accelerometers  were 
quite  different  from  those  from  ihe  aft 
cluster.  Although  some  events  and  seme 
frequency  bands  seem  to  have  a Gaussian  - 
like  distribution,  the  number  of  peaks 
counted  during  a 5 -second  interval  are 
too  few  to  represent  a steady-state 
excitation.  It  would  appear  that  the 
environment  measured  in  the  bow  of  the 
LCM8  is  an  intermittent  series  of  dis- 
crete excitations. 

A representation  of  the  overall 
vibration  of  the  LCM8,  shown  in  Figure 
13,  was  developed  from  the  VAIL  sum- 
maries. Comparison  of  Figure  13  with 
Figure  14  shows  that  the  peak  values 
of  the  former  are  much  lower  than  those 
in  Figure  14a.  This  result  is  caused 
by  the  operation  of  program  VAIL,  which 
tends  to  average  the  values  it  combines. 
Since  there  were  so  many  more  peaks 
recorded  in  the  aft  position  (Figure 
14b) , the  few  high  peaks  recorded  in 
the  forward  position  were  averaged  out. 


(a)  (b) 

Fig.  14  - LCM8  Excitation,  Forward  and  Aft 
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Generally,  the  amplitudes  given  for  the 
aft  position  (Figure  14b)  should  be  con- 
sidered to  represent  the  steady-state 
accelerations,  whereas  those  in  the 
forward  record  (Figure  14a)  are  the 
discrete,  intermittent  excitations. 
Figure  14b  is  a representation  of  the 
data  taken  from  the  accelerometers 
located  in  the  aft  portion  of  gunwale 
in  the  cargo  area.  The  longitudinal 
axis,  in  general,  shows  the  lowest  am- 
plitudes and  the  vertical  axit  the 
highest.  The  amplitudes  are  low,  never 
exceeding  0.4g.  Analysis  of  the  nature 
of  the  excitation  shows  that  in  the 
longitudinal  and  vertical  axes  there 
is  evidence  of  the  influence  of  discrete 
or  narrow-band  random  excitation  in 
every  frequency  band  above  40  Hz.  This 
excitation  is  mixed  with  broadband 
Gaussian  random  excitation.  The  excita- 
tion in  the  transverse  axis  appears  to 
be  only  broadband  random  with  a near- 
Gaussian  distribution. 

One  occurrence  as  the  LCM8  was 
approaching  the  beach  showed  an  inter- 
esting comparison  between  the  excitation 
in  the  forward  and  aft  positions.  There 
was  a sand  ridg*  far  enough  out  that 


the  craft  had  not  yet  reduced  speed  for 
beaching.  As  the  boat  struck  the  ridge, 
the  impact  was  great  enough  to  cause 
one  member  ot  the  group  to  lose  his 
footing  and  fall  to  the  deck.  The 
acceleration  time  histories  and  response 
spectra  derived  from  this  event  are 
showi  in  Figure  15  for  the  aft  position. 

In  the  longitudinal  and  vertical 
axes,  the  response  is  greatest  at  45  Hz 
in  the  forward  position,  but  the  aft 
shows  greatest  response  at  170  Hz.  In 
the  transverse  axis  , however,  the  for- 
ward response  is  different  in  shape, 
peaking  at  200  Hz,  and  the  aft  response 
is  essentially  f’.at  across  the  entire 
spectrum  at  0.8?.  The  highest  response 
is  1.25g  at  170  I!z  in  the  vertical  axis 
in  the  aft  position. 

Comparison  between  the  records  in 
the  forward  and  aft  positions  (Figure 
16)  shows  great  dominance  of  the  for- 
ward excitation  in  the  transverse  axis. 
Since  this  does  not  represent  a large 
number  of  peaks,  and  is  intermittent 
discrete  rather  than  steady-state 
excitation,  comparison  is,  perh.-ps , not 
logical.  The  high  amplitudes  of  the 
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Fig.  15  - LCM8  Hitting  the  Beach 
Response  Spectra 
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excitation  in  the  bow,  however,  should 
be  noted.  The  steady-state  excitation 
with  Gaussian  random  distribution  mea- 
sured in  the  aft  position  dominates  the 
other  two  aces. 

The  characteristic  excitation  found 
at  the  bow  of  the  LCM8  is  largely  inter- 
mittent. Farther  aft,  near  the  engines, 
there  is  a broadband  random  base  which 
is  Gaussian  in  distribution,  but  with  a 
generous  mixture  of  narrow-band  random 
excitation.  The  amplitudes  are  low. 

LCU 

The  measurements  were  made  on  the 
LCU  during  a trip  on  San  Diego  Bay.  The 
boat  ran  up  the  bay  for  a period,  then 
beached  at  normal  speed,  backed  off,  and 
returned  to  its  berth  at  the  Coronado 
Amphibious  Base.  Because  of  the  size 
of  the  LCU,  it  was  not  possible  to  in- 
strument both  forward  and  aft . One 
triaxial  cluster  was  located  about  10 
yards  forward  of  midships  on  a tiedown 
cleat  at  centerline.  The  others  were 
placed  aft  between  the  engines . The 
measurements  were  taken  during  a run 
at  flank  speed,  beaching,  backing  off, 
end  full  back. 

Figure  17  is  a composite  derived 
from  all  of  the  records  taken  on  the 
LCU.  While  the  amplitudes  in  the  ver- 
tical axis  generally  dominate,  the 
highest  values,  up  to  0.35g,  are  found 
in  the  transverse  axis  between  20  and 
80  Hz.  High  amplitudes  in  the  longi- 
tudinal axis  appear  above  1000  Hz. 
Otherwise,  the  longitudinal  axis  fre- 
quently shows  the  lowest  values. 

The  summaries  for  the  forward  and 
midships  positions  are  shewn  in  Figures 
18a  and  18b.  A comparison  reveals  that 
the  amplitudes  are  approximately  equal 
below  80  Hz;  however,  above  that  point 
the  record  for  the  midships  position, 
near  the  engines,  generally  dominates. 

It  also  shows  the  prominence  of  the 
transverse  axes  in  the  lower  frequen- 
cies. The  vertical  axis  generally 
dominates  in  the  middle  frequencies  up 
to  1000  Hz,  and  the  longitudinal  axis 
emerges  above  that  point . 

Figures  19a,  19b,  and  19c  provide 
a comparison  between  the  forward  and 
midships  records,  with  peak  amplitude 
in  each  axis  considered.  The  record 
for  midships  dominates  in  all  axes  and 
frequencies  except  in  the  700-  to  1000- 
Hz  frequency  band  in  the  vertical  axis. 

In  general,  the  characteristics 
discussed  under  the  LCM8  appear  in  the 
data  from  the  LCU,  although  not  so 


vividly.  Records  of  individual  events 
show  that  in  the  forward  position  the 
unusual  peak  distribution  across  the 
spectrum  appears,  but  only  when  normal 
or  slow  speeds  are  being  used.  During 
the  high-speed  events,  flank  speed  and 
nearly  full  back,  the  peak  distribution 
is  more  like  th'  normal  situation,  which 
is  exhibited  consistently  by  the  mea- 
surements made  near  the  engines.  The 
large  number  of  peaks  produced  during 
the  high-speed  events  dominates  the 
record  to  the  extent  that  the  VAIL  sum- 
maries do  not  exhibit  this  character- 
istic. It  seems  reasonable  to  conclude 
that  under  normal  speed  conditions  the 
excitation  on  the  forward  cargo  deck  of 
the  LCU  is  largely  discrete  and  inter- 
mittent in  nature,  but  at  high  speeds 
becomes  more  continuous  and  largely 
Gaussian  random. 

The  VAIL  composite  records  for  the 
forward  accelerometers  show  the  presence 
of  narrow-band  random  excitation  only 
in  the  vertical  axis  between  500  and 
1400  Hz;  at  the  midship  station,  how- 
ever, the  characteristics  of  the  narrow- 
band excitation  are  present  in  all  axes. 
In  the  longitudinal  direction,  there  are 
few  indications,  and  only  between  350 
and  1000  Hz.  In  the  transverse  axis, 
again,  there  are  few,  and  only  between 
240  and  700  Hz  in  the  40-  to  80-Hz  band. 
In  the  vertical  axis,  such  phenomena  are 
seen  in  all  frequency  bands  except  the 
lowest.  Comparisons  of  the  number  of 
peaks  reported  in  all  positions  show  the 
most  continuous  condition  to  be  in  the 
midship  position.  For  instance,  in  the 
longitudinal  record  from  the  forward 
position  16,057  peaks  are  shown  in  the 
highest  frequency  band,  whereas  57,235 
appear  in  the  same  band  of  the  longi- 
tudinal axis  in  the  midship  position. 

The  excitation  measured  in  the  LCU 
was  very  similar  to  that  measured  in 
the  ICMo.  Forward  there  was  a great 
deal  of  intermittent  discrete  excitation, 
whereas  the  midship  position  exhibited 
a dominative  broadband,  near-Gaussian 
excitation  with  some  narrow-band  excita- 
tion intermixed.  The  amplitudes  were 
generally  lower  in  the  LCU. 

SUMMARY 

A summary  of  the  land  events  mea- 
sured (Figure  20)  is  based  on  the  VAIL 
summaries . A comparison  of  Figure  2 
with  Figure  2 shows  that,  although  ti.?ri: 
is  little  difference  between  the  LVTi’7 
and  the  LCA2,  the  LVTP7  is  the  dominai  t 
vehicle.  The  vertical  axis  dominates 
up  to  350  Hz.  The  other  two  axes  have 
amplitudes  which  are  slightly  higher 
above  that  frequency;  however,  in  the 
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700-  to  1000 -Hz  band,  Che  longlCudinal 
and  verClcal  axes  are  approximately  of 
equal  amplicude.  A comparison  between 
Che  records  of  Che  LVTP7  and  Che  LCA2 
(Figure  21)  shorn  ChaC  only  In  Cwo  fre- 
quency bands  In  Che  records  of  Che 
longlCudinal  axis  does  Che  LCA2  shew 
higher  ampllCudes  Chan  Che  LVTP7. 

The  summary  of  wacer  events  (Figure 
22)  exhibits  acceleration  ampllCudes 
Chat  are  lower  Chan  chose  measured 
during  Che  land  events  The  transverse 
axis  predominates  In  six  of  Che  fre- 
quency bands,  whereas  Che  verClcal 
predominates  In  only  four.  Examination 
of  Che  records  of  Cue  various  vehicles 
shows  that  the  LVTP7  exhibited  generally 
higher  amplitudes  during  Che  water 
events  than  did  Che  LCU  or  LCM8 . How- 
ever, this  dominance  Is  not  complete. 
Only  In  Che  verClcal  axis  above  20  Hz 
Is  Che  LVTP7  completely  dominant. 

These  data  are  shown  In  Figure  23. 

The  virtually  complete  ascendency 
of  the  record  of  land  events  over  the 
water  evencs  Is  depleted  In  Figure  24. 
Only  In  the  low-frequency  band  of  Che 
transverse  axis  does  Che  amplitude  of 
Che  accelerations  measured  during  Che 
waCer  events  rise  above  those  measured 
on  land. 

CONCLUSIONS 

The  following  general  conclusions 
may  be  made: 
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When  a comparison  calibration  is  made  with  a test  and  a reference 
accelerometer  assembled  on  the  table  of  an  electromagnetic 
vibrator,  significant  error  can  result  if  the  comparison  is  made 
at  frequencies  where  appreciable  transverse  motion  of  the  vibrator 
occurs.  A new  method  is  described  in  which  lateral  motion  is 
restrained  by  using  annular  spring  diaphragms  to  support  the  two 
accelerometers  and  constrain  their  motion  to  a rectilinear  path. 
The  amplitude  of  the  calibrating  motion  is  amplified  by  driving 
the  system  at  resonance.  By  varying  the  thickness,  material  and 
outer  clamping  diameter  of  the  diaphragms  a large  number  of 
calibrating  frequencies  is  provided.  In  the  prototype  calibrator 
described,  these  frequencies  are  in  the  range  98-520  Hz,  and  the 
calibrating  motion  has  less  than  one  per  cent  total  harmonic 
distortion  and  less  than  one  per  cent  transverse  acceleration  for 
displacements  up  to  0.25  am  peak-to-peak.  The  corresponding  values 
of  peak  acceleration  are  47  ms"*  (98  Hz)  and  1330  ms"*  (520  Hz). 


INTRODUCTION 

In  comparison  methods,  the  test  accelero- 
meter is  calibrated  by  comparison  with  a 
reference  accelerometer  that  has  been  calibrated 
in  terms  of  the  relevant  national  standards  of 
measurement.  Laboratories  concerned  with 
vibration  measurement  mostly  use  comparison 
methods,  while  absolute  calibrations  are 
normally  performed  only  :,y  laboratories  main- 
taining primary  measurement  capability  [ 1 ]. 

The  two  accelerometers  to  be  compared  may 
be  assembled  back-to-back  on  a fixture  as  in 
Fig.  1(a),  or  the  comparison  rnay  be  made  with  the 
two  accelerometers  sidc-by-side  on  the  vibrator 
tabic  if  suitable  precautions  are  taken  ( 2 ). 

In  recent  years  increasing  use  is  being  made  of 
the  type  of  reference  accelerometer  designed  so 
that  the  test  accelerometer  can  be  attached 
dinectly  to  the  reference  accelerometer  as  in 
Fig. 1(b).  Descriptions  of  reference  accelero- 
meters of  this  type  and  discussion  of  their 
absolute  calibration  and  the  technique  of  their 
use  can  be  found  in  Ref.[  3 ] Parts  1 and  2. 
Procedures  and  techniques  of  comparison 
calibration,  and  precautions  necessary  to  ensure 
valid  comparison  of  accelerometers,  are  dis- 
cussed in  Ref. [ 3 ] . 

Ideally,  a calibrator  should  generate 
rectilinear  simple  harmonic  motion  and  the  two 


accelerometers  should  be  subjected  to  the  same 
amplitude  of  this  notion.  Departures  from  this 
ideal  that  give  trouble  in  practice  are  dis- 
tortion and  non-rectilinearity  of  the  motion. 
Such  defects  in  the  calibrating  notion  nay 
cause  errors  because  they  can  produce  different 
contributions  to  the  outputs  of  the  two 
accelerometers. 


FIGURE  1 Back-to-Back  Assemblies 


In  practice  the  calibration  may  be  con- 
fined to  frequency  ranges  in  which  the 
vibrator  behaves  satisfactorily,  if  such  ranges 
have  been  identified  in  a prior  investigation 
of  the  performance  of  the  calibrator  under  the 
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sane  conditions  of  loading  and  use  as  in  the 
calibration.  Alternatively  the  quality  of  the 
vibratory  notion  nay  be  nonitored  during  the 
actual  calibration.  For  example,  in  the  com- 
parison calibrator  described  in  Ref. [ 3]  Part  3, 
the  input  notion  is  determined  by  a set  of  three 
reference  accelerometers  which  are  arranged  and 
operated  to  indicate  non-recti linearity  of  the 
notion. 

Some  laboratories  have  constructed 
vibrators  designed  specifically  as  vibration 
input  generators;  for  example,  Dinoff  ( 4 ] used 
an  air  bearing  to  constrain  the  notion  of  an 
electrodynanic  vibration  standard.  More 
recently  there  has  been  coanercial  development 
of  vibrators  designed  for  use  as  calibrators. 

In  the  present  paper,  a different  approach  is 


described  in  which  no  unusual  quality  is 
demanded  of  the  vibrator.  Instead,  the  notion 
derived  from  an  ordinary  vibrator  is  con- 
strained laterally  and  amplified  axially  by  a 
resonance  fixture  interposed  between  the 
vibrator  and  the  back-to-back  accelerometers. 

In  Australia  as  elsewhere,  establishments 
concerned  with  vibration  measurement  maintain 
reference  accelerometers  i.liich  they  use  to 
check  the  calibration  of  their  'working' 
accelerometers.  The  new  calibrator  will  be  used 
to  compare  these  reference  accelerometers  in  the 
low  frequency  range  (100-S00  Hz)  with  a 
reference  accelerometer  that  has  been  calibrated 
at  the  National  Standards  Laboratory  by  the 
interferometric  methods  developed  by  Goldberg 
[S.6.7J. 


FIGURE  2 Cross  Sectional  View  of  Calibrator 
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BASIC  DESIGN  OF  CALIBRATOR 

Referring  to  the  cross  sectional  diagraa 
in  Fig. 2,  the  test  accelerometer  A and  reference 
accelerometer  B are  fixed  at  opposite  ends  of 
the  capsule  C with  their  sensitive  axes  coaxial 
with  the  capsule.  The  capsule  is  supported  on 
three  annular  spring  diaphragms  D,  E and  F that 
are  clamped  around  their  outer  edges  in  a 
housing  G fixed  on  a concrete  block  H in  which 
the  vibrator  J is  fixed.  The  inner  edges  of  the 
diaphragms  are  clamped  to  the  capsule,  which  is 
connected  through  the  driving  assembly  K to  the 
moving  system  of  the  vibrator. 

In  principle  only  two  diaphragms  are 
necessary,  one  near  each  end  of  the  capsule: 
the  third  is  incorporated  so  that  when  the  upper 
diaphragms  and  clamping  rings  are  removed  and 
others  substituted  to  change  the  resonance 
frequency,  as  explained  later,  the  lowest 
diaphragm  holds  the  capsule  in  position. 


The  frequency  of  the  input  signal  to  the 
vibrator  is  adjusted  to  the  fundamental  axial 
resonance  frequency  of  the  system  and,  with  the 
vibration  amplitude  set  at  the  required  level, 
the  electrical  outputs  of  the  accelerometers 
(or  accelerometer,  amplifier,  readout  systems) 
are  compared.  This  is  repeated  at  other 
frequencies  as  required. 

DESCRIPTION  OF  THE  CALIBRATOR 

Fig. 3 is  a photograph  of  the  calibrator 
and  Fig. 4 shows  a typical  assembly  of  the  moving 
system. 

Capsule 

The  capsule  is  of  aluminius  100  mm  long, 

60  am  diameter  and  mass  1.1  kg  The  accelero- 
meters to  be  compared  are  screwed  into  adapt- 
ors L threaded  to  suit  different  makes  of 
accelerometer.  The  test  accelerometer  is  at 
the  immediately  accessible  upper  end  of  the 


*A 


% 


FIGURE  3 The  Back-to-Back  calibrator 


the  operating  frequency  range  of  the  calibrator. 


Driving  Assembly 

The  lower  end  of  the  capsule  is  connected 
to  the  Moving  system  of  an  elect  rod ynamic 
vibrator  of  peak  force  rating  140  N by  the 
driving  asseably  K,  which  is  designed  to  clear 
the  reference  accelerometer  and  its  cable,  and 
to  permit  adjustments  for  alignment  and 
centring. 

Spring  Diaphragms  and  Clasping  Rings 

Two  sets  of  diaphragms  havo  been  provided, 
one  set  made  of  stainless  steel  and  the  other  of 
phosphor  bronze.  Each  set  comprises  twelve 
diaphragms:  four  values  of  thickness  (t«  • 

0.7  mm,  tj  * 0.9  mm,  ti  » 1.2  wm  and  U * 1.6  mm) 
for  each  of  the  three  diaphragms. 

The  effective  outer  clamping  diameter  of 
the  upper  pair  of  diaphragms,  D and  E,  is 
changed  by  using  different  sets  of  the  clamping 
rings  P.  Three  sets  of  outer  rings  provide 
values  of  di  * 107  mm,  dj  = 117  nm,  dj  > 127  am. 
For  the  lowest  diaphragm  F,  the  outer  clamping 
diameter  is  always  di  * 127  nm.  The  effective 
inner  clamping  diameter  (d« ) of  all  diaphragms 
is  50  mm. 


FIGURE  4 Typical  Assembly  of  the 
Moving  Systeii 


capsule.  The  reference  accelerometer,  which 
normally  remains  undisturbed  while  various  test 
accelerometers  are  cumpar<-d  with  it,  is  at  the 
lower  end.  Additional  positions  are  provided 
in  the  sides  of  the  capsule  for  accelerometers  M 
to  monitor  lateral  motion,  or  targets  N for  use 
with  optical  proximity  probes  for  checking 
lateral  motion. 

For  the  prototype  calibrator  two  capsules 
were  designed,  one  of  aluminium  and  one  ceramic. 
The  advantages  of  using  a ceramic  material  such 
as  alumirium  oxide  for  vibration  tables  have 
been  stated  in  Ref.  [ 3 ] Part  3.  Because  of  the 
high  stiffness,  low  density  and  damping  proper- 
ties of  this  material,  the  frequency  at  which 
wave  motion  becomes  significant  is  above  the 
usual  range  of  calibration  frequencies.  This 
material  has  the  further  advantage  that  it  is  an 
electrical  insulator,  hence  the  two  accelero- 
meters are  electrica.il y insulated  from  one 
another  and  from  ground  except  where  connections 
are  deliberately  introduced.  The  development  of 
the  prototype  and  assessment  of  its  mechanical 
performance  described  in  the  present  paper  have 
been  carried  out  with  the  aluminium  capsule. 

The  calculated  axial  resonance  frequency  of  the 
aluminium  capsule,  about  25  kHz,  is  well  above 


When  tightening  the  hexagonal  nut  R at  one 
or  other  end  of  the  capsule  to  load  the  inner 
clamping  rings,  a special  fixture  is  used  so 
that,  when  the  external  torque  (55  Nm)  is 
applied,  the  reaction  is  transferred  directly  to 
the  housing,  and  the  capsule  is  not  strained  in 
torsion. 

To  ensure  uniformity  of  clamping  on  the 
outer  clamping  rings  the  screws  holding  the 
upper  and  lower  flanges  are  uniformly  tightened 
to  a specified  torque  (30  Nm) . 

Housing 

The  shape  of  the  cast  iron  housing  G is 
best  seen  in  Fig. 3.  The  three  large  openings  in 
the  sides  are  for  access  to  the  driving  assembly. 
Four  smaller  holes,  two  of  which  can  be  seen  in 
the  photograph,  align  with  the  cross  holes  in  the 
capsule.  These  provide  access  for  the  proximity 
probes  mentioned  earlier,  and  are  used  as  exits 
for  the  cables  of  transverse  monitoring  accelero- 
meters . 

The  housing  is  located  on  an  adaptor  plate 
S designed  so  that,  when  the  drive  is  disconnect- 
ed and  the  housing  removed  for  access  to  the 
reference  accelerometer  or  the  lowest  diaphragm, 
the  system  can  be  reassembled  without  loss  of 
alignment. 

The  calibrator  is  supported  on  the  concrete 
block  H of  mass  160  kg  resting  on  vibration 
isolators. 

RANGE  AND  PERFORMANCE  OF  CALIBRATOR 

A series  of  tests  has  been  made  to 
evaluate  the  range  of  the  calibrator  and  the 
quality  of  the  calibrating  motion.  Ideally  the 
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calibrating  Notion  should  be  undistorted  recti- 
linear simple  harmonic  notion.  For  practical 
purposes  the  aia  has  ’tee:'  to  determine  the 
range  in  which  distortion  and  transverse  Notion 
are  within  acceptable  limits. 

Frequency  San ye 

The  192  basic  combinations  available  are 
shewn  in  Table  1.  Additional  combinations 
could  be  naoe  by  not  choosing  the  sue  naterial 
and  thickness  for  the  upper  two  diaphragns. 


acceleroaeter  mounted  01  the  calibrator,  using 
an  XY  recorder  calibrates  for  frequency  against 
a frequency  counter.  The  resonance  frequency 
was  read  from  these  response  i tcords. 

A representative  sequence  of  12  resonance 
responses  recorded  in  a succession  of  settings 
is  illustrated  in  Fig. 5.  The  sharpness  of 
resonance  is  indicated  by  Q * 40,  calculated 
directly  from  the  higher  frequency  records. 

The  value  of  Q would  appear  less  if  calculated 


TABLE  1 


Basic  Combinations  of  Diaphragm  Material, 
Thickness  and  Clamping  Diameter. 


Lowest  Diaphragm 

Upper  Pair  of  Diaphragms 

Material 

Thickness 

Material 

S2S  «**-» 

SS* 

ti 

SS 

di  ti  ti  ti  t< 

di  ti  tj  tj  t« 

dj  ti  ta  ti  t< 

SS 

t» 

SS 

12  combinations  as  above 

SS 

t, 

SS 

12  combinations 

SS 

t« 

SS 

12  combinations 

PBt 

t.  ti  ti  ti 

PB 

48  combinations  as  above 

SS 

t.  ti  tj  ti 

PB 

48  combinations 

PB 

ti  t5  t3  t« 

SS 

48  combinations 

*SS  stainless  steel 

tPB  phosphor  bronze 


For  the  experimental  evaluation  of  the 
calibration  frequency  range,  the  calibrator  was 
set  up  and  operated  at  a selection  of  about 
half  of  the  first  96  combinations  listed  in  the 
table;  that  is,  using  the  same  naterial  for  the 
upper  pair  of  diaphragms  a«  for  the  lowest 
diaphraga. 

At  each  setting  the  input  frequency  to  the 
driving  vibrator  was  swept  through  a suitable 
range  and  a record  made  of  the  output  of  an 


from  the  lower  frequency  records,  due  to  the 
influence  of  the  5 Hz  bandwidth  of  the  detector. 

The  numerical  values  of  the  SO  resonance 
frequencies  observed  in  these  tests  are  plotted 
in  Fi^.6.  It  can  be  seen  from  the  spread  of 
the  test  results,  which  represent  about  a 
quarter  of  the  basic  combinations  shown  in 
Table  1,  that  a resonance  frequency  for 
calibration  purposes  could  be  provided  in  the 
vicinity  of  any  required  frequency  in  the  range 
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98-520  Hz. 

If  a particular  precise  value  of  the 
frequency  is  required,  a combination  can  be 
chosen  having  a resonance  frequency  slightly 
higher,  and  mass  added  to  the  capsule  to  lower 
the  frequency  to  the  desired  value.  For 
example,  in  an  observation  at  a resonance 
frequency  of  260  Hz,  addition  of  a 200  gram  mass 
reduced  the  resonance  frequency  to  250  Hz; 
addition  of  500  gram  reduced  the  resonance 
frequency  from  260  to  230  Hz. 

Distortion 

Distortion  in  the  Form  of  harmonics  of  the 
calibrating  frequency  may  arise  from  distortion 


of  the  electric  signal  applied  to  the  vibrator 
or  from  non-linear  stiffness  of  the  mechanical 
system.  The  input  signal  distortion  is  readily 
measured  and  has  been  found  to  be  negligible. 

Distortion  attributable  to  the  mechanical 
system  has  been  investigated  by  driving  the 
calibrator  at  a fixed  amplitude  and  analysing 
the  signal  from  an  accelerometer  mounted  on  the 
top  face  of  the  capsule.  The  analysis  was  made 
with  a heterodyne  analyser,  sweeping  over  a wide 
range  using  a 5 Hz  b.  •'•width.  With  the  set-up 
used,  distortion  components  could  be  measured 
down  to  0,1  per  cent  with  a resolution  of  0.05 
per  cent. 
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FIGURE  6 Observed  Resonance  Frequencies 
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These  tests  showed  that  the  total  harmonic 
distortion  was  less  than  one  per  cent  if  the  dis- 
placement was  liaited  to  a maximum  value  of  0.2S 
an  peak-to-peak. 

One  aanifestation  of  the  mechanical  non- 
linearity has  been  a decrease  in  resonance 
frequency  with  increase  in  displacement 
amplitude.  In  an  experiment  made  with  an  im- 
proved method  of  clamping  the  diaphragms,  the 
resonance  frequency  increased  with  increased 
amplitude,  which  is  in  accord  with  the  theory 
for  centre-loaded  diaphragms  with  perfect 
clamping  ( 8 ]. 

Tests  with  matched  accelerometers,  one  at 
each  end  of  the  capsule, indicated  no  phase  or 
amplitude  difference  within  the  limits  of  the 
frequency  range  of  the  calibrator. 

Another  form  of  distortion  of  the  accelera- 
tion signal  results  from  external  vibration. 

With  the  calibrator  installed  in  the  manner 
described  earlier  there  has  been  no  distortion 
attributable  to  site  vibration. 

Acceleration  Range 

The  acceleration  range  is  determined  by 
the  acceptable  displacement  amplitude,  because 
the  distortion  resulting  from  the  non-linear 
stiffness  of  the  system  becomes  more  serious 
with  increasing  displacement.  With  the  calibra- 
tor in  its  present  form  and  the  displacement 
limited  to  0.25  mm  peak-tc-peak,  the  maximum 
acceleration  (peak)  is  47  ms  *’  at  98  Hz  and 
1330  ms-1  at  520  Hz. 

Transverse  Motion 

The  transverse  motion  was  measured  using  a 
fibre-optics  proximity  probe  to  determine  the 
displacement  of  the  capsule  relative  to  the 
probe.  The  reference  faces  on  the  capsule  were 
polished  brass  targets  N fitted  as  shown  in 
Figs.  2 and  4.  Observations  were  made  with  two 
targets  having  their  faces  vertical  and  mutually 
perpendicular.  The  probe  was  supported  indepen- 
dently of  the  calibrator  in  such  a way  that  it 
functioned  as  an  inertial  reference  for  the 
purpose  of  these  measurements. 

The  probe  was  calibrated  at  each  frequency 
against  the  accelerometer  used  to  measure  the 
axial  motion.  The  calibration  procedure  was  to 
mount  the  probe  so  that  it  measured  the  vertical 
displacement  of  the  accelerometer  whilst 
simultaneously  the  accelerometer  output  and  the 
frequency  were  observed.  This  was  repeated  at 
several  amplitudes  to  check  system  linearity. 

The  displacement  probe  was  then  re-positioned 
to  measure  the  transverse  motion. 

These  observations  showed  that  trans- 
verse displacement  of  the  capsule  throughout  the 
frequency  range  of  98-520  Hz  was  less  than  one 
per  rent  of  the  axial  displacement. 


CONCLUDING  REMARKS 

The  performance  tests  of  the  calibrator 
described  in  this  paper  indicate  that  unwanted 
transverse  motion  can  be  constrained  by  using 
annular  spring  diaphragms.  The  use  of  such 
diaphragms  makes  it  possible  to  design  the 
system  for  operation  at  a resonance  frequency  in 
a useful  range  for  calibration  purposes,  and  the 
incorporation  of  means  of  varying  the  stiffness 
properties  of  the  set  of  diaphragms  makes  other 
resonance  frequencies  available. 

If  one  per  cent  total  harmonic  distortion 
and  one  per  cent  transverse  acceleration  are 
considered  to  be  acceptable,  the  range  of  the 
calibrator  in  its  present  form  is  determined  by 
a displacement  limitation  of  0.25  m peak-to- 
peak  in  the  range  98-520  Hz.  The  corresponding 
limit  of  peak  acceleration  is  47-1330  ms"1  (or 
5-136  g,  where  g * acceleration  due  to  gravity). 

A further  development  now  in  progress  will 
provide  for  a set  of  diaphragms  of  increased 
thickness  (t  » 3.2  mm)  and  improved  clamping. 
This  is  expected  to  raise  the  frequency  range  to 
about  1 kHz,  which  is  probably  the  practical 
limit  for  this  particular  design. 

A smaller  and  simpler  version  of  the 
calibrator  is  being  developed  for  use  at  test 
sites  to  provide  for  comparison  calibrations  at 
tW5  or  three  frequencies. 
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The  transverse  sensitivity  characteristics  of  the  test  accelerometer 
are  detrained  by  mounting  the  accelerometer  with  its  sensitive  axis 
vertical,  on  a horizontal  surface  which  is  caused  to  oscillate  with 
simple  harmonic  motion  of  constant  acceleration  amplitude  A in  any 
direction  6 in  the  horizontal  plane.  In  the  prototype  calibrator 
described,  the  horizontal  surface  is  the  free  end  of  a vertical 
cantilever  beam  that  is  driven  in  transverse  flexural  vibration  by 
two  small  electrodynamic  vibrators  arranged  and  operated  in  such  a 
way  that  by  varying  the  relative  amplitudes  of  the  driving  forces 
any  desired  direction  6 can  be  set.  '.V;  acceleration  vector  A(6) 
is  determined  by  measuring  its  two  orthogonal  horizontal  components 
with  a pair  of  monitoring  accelerometers  attached  to  the  beam  adjacent 
to  the  test  accelerometei . The  procedure  is  described  for  determining 
the  relationship  between  the  transverse  sensitivity  ratio  (TSR)  and  6, 
and  hence  evaluating  the  two  quantities  of  most  practical  significance:  | 
the  maxima  TSR,  and  the  orientation  of  the  accelerometer  for  minimum 
TSR.  Examples  are  given  of  typical  transverse  calibrations  at  100  Hz 
and  A * 200  ms"* , with  0 stepped  through  10  degree  intervals.  Sources 
of  error  are  discussed  and  uncertainties  estimated. 


NOMENCLATURE 

A « amplitude  of  transverse  (ms-1 ) 
acceleration 

f * frequency  of  transverse  (Hz) 
vibration 

6 ■ direction  of  vibration  (deg) 
in  xy  plane 

TSR  = transverse  sensitivity 
ratio 

Z « identification  of  test 
accelerometer 

Ez  * output  of  test  (V) 

accelerometer 

S.  * axial  sensitivity  of  (Vm*1  s’ ) 

2 test  accelerometer, 
along  axis  Oz 

Cz  « sensitive  axis  of 
test  accelerometer, 
normal  to  mounting 
surface  (see  Fig.l) 

0(TSRmin)  « value  of  0 for  which  (deg) 
TSR  is  minimum 


X,Y  » identifications  of 

transverse  monitoring 
accelerometers 

0x,0y  * co-ordinate  directions 
in  horizontal  plane 
(see  Fig.l) 

Ex,Ey  ■ outputs  of  X and  Y (V) 

* Y accelerometers 

respectively 

SY,SV  ■ axial  sensitivities  (Vm* 1 s’ ) 

* Y of  X and  Y 

accelerometers 

respectively 

Other  symbols  are  defined  in  context  as  required 
INTRODUCTION 

Accelerometers  are  in  widespread  use  for 
the  measurement  of  the  mechanical  vibration  of 
vehicles,  structures,  machinery  and  many  other 
test  objects.  In  the  most  common  type,  a piezo- 
electric element  having  a seismic  mass  attached 
to  one  end,  has  its  other  end  fixed  to  the  body 
of  the  accelerometer  which  is  attached  to  the 
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object  whose  vibration  is  to  be  measured.  The 
assembly  is  designed  to  have  an  axis  of  symmetry 
normal  to  its  base,  and  is  intended  to  measure 
vibratory  acceleration  applied  in  the  direction 
of  that  axis,  which  is  referred  to  as  the 
sensitive  axis  of  the  accelerometer. 

If  an  accelerometer  is  being  used  to 
measure-vibration  in  the  direction  of  its 
sensitive  axis  while  it  is  also  being  subjected 
to  acceleration  having  components  normal  to  that 
axis  (i.e.  transverse  acceleration)  it  is 
desirable  that  the  transverse  acceleration  should 
not  contribute  to  the  electrical  output  of  the 
accelerometer.  This  is  unfortunately  an  ideal 
which  is  not  usually  attained  in  practice. 

Minor  inaccuracies  involved  in  the  processes  of 
manufacturing  and  assembling  the  various 
components  of  the  accelerometer,  and  some  non- 
uniformity and  asymmetry  of  the  properties  of 
the  materials  used,  cause  the  transverse 
sensitivity  of  most  accelerometers  to  be  in  the 
range  from  less  than  one  to  about  ten  per 
cent  of  the  axial  sensitivity.  Neglect  of  this 
factor  in  vibration  measurement  may  result  in 
significant  error,  therefore  it  is  essential  to 
determine  the  transverse  sensitivity  character- 
istics of  accelerometers. 

The  specific  objects  of  transverse 
calibration  are  to  determine: 

(i)  the  maximum  transverse  sensitivity  ratio, 
denoted  by  TSRmax, 

(ii)  the  transverse  angular  direction  on  the 
accelerometer  that  corresponds  to  the 
minimum  transverse  response  denoted  F < 
e(TSRmin). 

In  a simplified  conceptual  model  of  an 
accelerometer  which  assumes  that  the  transverse 
response  is  attributable  to  a simple  lack  of 
coincidence  of  the  sensitive  axis  Oz,  as  defined, 
and  some  axis  Oz'  along  which  the  sensitivity  is 
a maximum  (see  Fig.l),  the  direction  of  TSRmax 
lies  in  a plane  containing  these  two  axes,  and 
e(TSRmin)  is  normal  to  this  plane. 

When  an  accelerometer  is  to  be  used  to 
measure  vibration  in  a specified  direction, 
TSRmax  is  a measure  of  the  greatest  output 
likely  to  arise  from  vibrations  in  other 
directions.  If  the  direction  of  a particular 
source  of  transverse  acceleration  is  known,  then 
knowledge  of  0(TSRmin)  will  allow  the  accelero- 
meter to  be  oriented  so  that  the  contribution 
from  transverse  components  is  greatly  reduced. 

In  the  simplified  model,  for  example,  aligning 
e(TSRmin)  within  ± 6 deg  of  the  direction  of 
the  transverse  acceleration  will  reduce  the 
effect  by  a factor  of  about  10,  and  a ± 2 deg 
alignment  will  reduce  it  by  a factor  of  about 
30,  with  respect  to  the  maximum  possible  if  no 
alignment  precautions  are  taken. 

This  paper  describes  r new  and  simple 
method  of  transverse  calibration.  Although  the 
paper  refers  to  piezoelectric  accelerometers, 
obviously  the  principle  is  applicable  to  the 


transverse  calibration  of  vibration  transducers 
in  general. 

PRESENT  METHODS  OF  TRANSVERSE  CALIBRATION 

The  transverse  sensitivity  ratio  TSR  is 
defined  as  E,  /Ea , where  Et  is  the  output  of  a 
pickup  when  oriented  with  its  sensitive  axis 
transverse  to  the  direction  of  the  input  motion, 
and  Ea  is  the  output  when  the  sensitive  axis  is 
aligned  in  the  direction  of  the  same  input 
motion  [ 1 ).  The  transverse  sensitivity 
calibration  is  performed  by  applying  a transverse 
sinusoidal  motion  at  a single  frequency  bsiow 
500  Hz  where  it  is  known  that  the  transverse 
excitation  is  at  least  100  times  that  in  the 
direction  of  the  sensitive  axis  [ 2 ).  Because 
the  transverse  sensitivity  is  known  to  vary 
considerably  with  the  direction  of  oscillation 
in  the  transverse  plane,  the  accelerometer  is 
rotated  about  its  sensitive  axis  through  360  deg 
in  increments  of  45  deg  or  less  to  determine  the 
maximum  transverse  response. 

Rockwel 1 and  Ramboz  ( 3 ) have  described 
transverse  calibration  methods  which  account  for 
the  vector  components  of  test  accelerometer 
output  due  to  motion  ir.  other  than  the  desired 
direction.  Ramboz  [ 4 ] described  a development 
in  which  the  full  polar  plot  of  transverse 
sensitivity  is  derived  from  observations  at  only 
four  values  of  0.  Finley  ( 5)  described  a 
comparison  method  using  a reference  accelero- 
meter whose  transverse  sensitivity  had  been 
carefully  pre-determined . The  comparison  was 
made  at  eight  angular  settings  about  the 
accelerometer  axis  using  an  octagonal  fixture  as 
shown  in  Fig. 2(a).  In  Ref.  [6]  a mechanically 
driven  apparatus  is  described  which  provides  for 
rotation  of  an  accelerometer  about  its  sensitive 
axis  while  it  is  being  vibrated  transversely. 

In  the  method  described  in  this  paper. 
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physical  resetting  of  the  acceleraeeter  about 
its  sensitive  axis  is  not  necessary;  instead 
the  accelerometer  regains  attached  to  a 
horizontal  surface  in  the  xy  plane  throughout 
as  shown  in  Fig. 2(b)  and  the  surface  itself  is 
caused  to  oscillate  parallel  to  any  direction 
9 in  the  plane. 


FIGURE  2 Basic  Difference  Between 

(a)  Existing  Method  and 

(b)  New  Method 


PRINCIPLE  OF  THE  NEW  TRANSVERSE  CALIBRATOR 

The  new  method  makes  use  of  the  fact  that 
two  rectilinear  coplanar  sinusoidal  motions  of 
the  same  frequency,  acting  in  phase  along 
mutually  perpendicular  directions,  have  as  their 
resultant  a rectilinear  motion  acting  in  a 
direction  which  is  determined  by  the  relative 
amplitudes  of  the  two  input  components  and  which 
can  be  varied  by  adjusting  the  relative  am- 
plitudes. Referring  again  to  Fig. 2(b)  the  input 
components  acting  along  Ox  and  Oy  are  adjusted 
to  produce  motion  along  a line  at  any  required 
angle  9 in  the  xy  plane. 

This  is  implemented  in  the  new  calibrator 
in  the  manner  shown  schematically  in  Fig. 3.  A 
vertical  cantilever  beam  of  circular  or  annular 
cross  section  is  clamped  near  its  lower  end. 

The  test  accelerometer  is  attached  to  the  upper 
end  of  the  beam  with  its  sensitive  axis  Oz  co- 
axial with  the  beam.  The  beam  is  driven  in 
transverse  vibration  by  two  small  electrodynamic 
vibrators  having  their  axes  mutually  perpendicu- 
lar and  in  a horizontal  plane  a short  distance 
above  the  clamp.  The  vibrators  are  operated  in 
phase  from  a single  oscillator  which  is 
adjusted  to  the  transverse  resonance  frequency 
of  the  system.  By  varying  the  relative 
amplitudes  of  the  two  driving  forces  the  test 
accelerometer  can  be  subjected  to  transverse 


FIGURE  3 Schematic  Diagram  of 
New  Transverse  Calibrator 


vibration  in  any  direction  in  the  plane  xy  of 
its  mounting  face. 

The  TSR  is  determined  by  measuring  the  output 
(Ez)  of  the  test  accelerometer  when  a transverse 
acceleration  of  known  amplitude  A is  applied  at 
successive  values  of  8.  The  acceleration  vector 
X(8)  is  determined  from  its  orthogonal  horizon- 
tal components,  which  are  measured  with  a pair 
of  monitoring  accelerometers  X and  Y attached  to 
flat  surfaces  ground  on  the  beam  near  the  test 
accelerometer.  The  outputs  Ex  and  Ey  of  these 
accelerometers  are  connected  to  an  oscilloscope 
so  that  an  analogue  of  the  transverse  motion  is 
displayed  as  a Lissajous  figure. 

DESCRIPTION  CF  NEW  TRANSVERSE  CALIBRATOR 

In  the  prototype  calibrator  illustrated  in 
Figs.  3 and  4 the  cantilever  beam  is  a 25  mm 
diameter  steel  bar  of  total  length  0.7  m.  The 
beam  is  clamped  near  its  lower  end  in  a housing 
which  also  supports  the  two  vibrators  and  is 
bolted  to  a vertical  face  of  a heavy  concrete 
block.  For  clarity  the  photograph  shows  the 
calibrator  removed  from  the  concrete  block. 

The  vibrators,  which  are  of  a commercially 
available  type  with  a peak  force  rating  8.5  N, 
have  their  axes  in  a horizontal  plane  0.1  m 
above  the  upper  face  of  the  clamp.  The  signals 
to  the  two  vibrators  are  obtained  through  a two- 
channel  amplifier  (see  Fig. 5)  from  a common 
oscillator  operated  at  the  frequency  of  reson- 
ance of  the  beam.  Each  channel  has  an  amplitude 
control  and  a t 90  deg  variable  phase-shift 
control.  The  latter  is  necessary  to  compensate 
for  the  effects  of  small  variations  in  the 
transverse  mechanical  stiffness  of  the  system  as 
the  direction  of  vibration  is  altered;  these 
variations  cause  differing  phase  shifts  in  the 
two  component  directions,  giving  rise  to 
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r ACCELEROMETER  X 


FIGURE  S Block 


FIGURE  4 Prototype  Calibrator 
Removed  from  Supporting  Block 


elliptic  rather  than  rectilinear  motion  in  the 
transverse  plane. 

The  test  accelerometer  Z is  mounted  on  the 
upper  end  surface  of  the  beam  with  its  sensitive 
axis  coaxial  with  the  beam.  Also  at  the  upper 
end  are  four  flat  ground  surfaces,  orthogonal 
to  each  other  and  to  the  flat  end  of  the  beam, 
on  two  of  which  are  attached  the  monitoring 
accelerometers  X and  Y.  These  measure  the 
transverse  acceleration  components  Acos0  and 
Asine  respectively.  The  effects  of  transverse 


Diagram  of  System 

excitation  of  the  X and  Y accelerometers  them- 
selves and  of  their  offset  6L  from  the  xy  plane 
(Fig. 3)  are  discussed  in  a later  section. 

The  outputs  Ex  and  Ey  of  the  two 
monitoring  accelerometers  are  applied  to  the  'X' 
and  'Y*  plates  respectively  of  an  osci lloscope, 
thus  A and  ecan  be  measured  by  using  the  circles 
and  lines  of  a polar  scale  fitted  over  the 
oscilloscope  screen  (Fig.S). 

With  the  test  and  monitoring  accelero- 
meters in  position,  the  resonance  frequency  can 
be  adjusted  to  near  the  desired  value  by  varying 
the  clamping  position  of  the  bar.  For  fine  ad- 
justment a small  ’tuning'  mass  in  the  form  of  a 
collar  can  be  clamped  at  any  position  (Li , 

Fig. 3)  between  the  driving  collar  and  the 
monitoring  accelerometers. 

With  the  system  as  illustrated  and  the 
resonance  frequency  adjusted  to  100  Hz  the 
free  length  of  the  beam,  (Li ),  is  approximately 
0.4  m. 

CALIBRATION  TECHNIQUES 

With  the  test  accelerometer  in  place,  its 
orientation  with  respect  to  Ox  (Fig. 2(b))  is 
noted,  making  use  of  a suitable  mark  or  feature 
on  the  accelerometer  and  using  a protractor 
which  fits  over  the  top  of  the  bar.  The  drive 
is  applied,  and  while  observing  the  oscilloscope 
screen  the  phase  and  amplitude  of  the  driving 
signals  are  adjusted  to  produce  on  the  screen  a 
straight  line  of  length  representing  the  desired 
acceleration  amplitude  A and  inclined  nominally 
at  the  direction  0 as  read  from  the  polar  scale. 

The  output  E£  of  the  test  accelerometer 
and  the  output  Ey  of  the  monitoring  accelero- 
meter X are  displayed  on  a two-beam  oscillo- 
scope (Fig. 5),  from  which  is  obtained  the  phase 
of  Ez  with  reference  to  the  positive  direction  of 
the  applied  transverse  acceleration.  When  0 
is  in  the  vicinity  of  90  deg  or  270  deg,  and 
output  Ey  is  consequently  very  small,  output 
Ey  may  be  substituted.  The  outputs  Ey.Ey  and 
E 7 are  also  read  from  digital  voltmeters.  This 
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is  repeated  for  say  10  deg  increments  of  6,  and 
finally  corresponding  values  of  TSR  and  6 are 
calculated  using  the  following  equations: 


where 


TSR  * Ez/SzA 


6 « arc  tan 


(Ey/Sy) 


(Ex/Sx) 

A - {(Ex/sx)J  ♦ (EY/Sy)1  I*1 


(1) 

(2) 

(3) 


and  S_,  Sy  and  Sz  are  the  axial  sensitivities 
of  the  three  accelerometers. 


The  convention  adopted  ror  Eq.(l)  is  that 
the  sign  of  TSR  is  « or  - when  Ez  is  in-phase  or 
anti-phase  respectively,  with  respect  to  the 
applied  transverse  acceleration  (see  Fig. 6(a), 
(b)).  In  applying  Eq.(2),  Ex  is  assigned  ♦ or 
- when  it  is  in-phase  or  anti-phase  respectively 
with  respect  to  Ey.  Two  quadrants  are  sufficient 
khen  using  this  new  technique. 

The  test  data  are  conveniently  reduced  by 
using  a staple  computer  program  which  can  be 
made  to  yield  tabulated  and/or  plotted  results. 

A more  rapid  check  may  be  made  by 
scanning  through  the  range  of  6,  while  adjusting 
amplitude  and  phase  to  maintain  a constant  value 
of  A,  until  a minimum  value  of  E7  is  found.  The 
outputs  Ex  and  Ey  are  read  and  the  angle 
6(TSRmin)  calculated.  TSRmax  may  then  be  found 
by  increasing  or  decreasing  6 by  90  deg. 

RANGE  AND  PERFORMANCE  OF  THE  PROTOTYPE 
CALIBRATOR 


The  prototype  was  designed  for  calibrat- 
ions at  100  Hz.  With  the  free  length  li  * 0.4  m 
and  a test  accelerometer  of  mass  40  grim,  clamp- 
ing the  adjustable  collar  of  mass  104  gram  at 
different  positions  U along  the  beat,  allowed 
the  resonance  frequency  to  be  varied  over  the 
range  95  Hz  to  103  Hz, 

The  total  range  of  resonance  frequency 
available  with  the  prototype  by  varying  It  in 
the  range  0.3  n to  0.6  m and  U over  its  full 
range  is  39  Hz  to  155  Hz. 

The  maximum  acceleration  amplitude 
attainable  at  100  Hz  with  the  prototype  as 
described  is  200  ms'1 

Two  examples  of  calibration  results  are 
presented.  In  Fig. 6(a)  the  full  line  shows  TSR 
versus  0 for  a relatively  inexpensive  accelero- 
meter having  a high  value  of  TSRmax  (13  per 
cent).  For  comparison  the  broken  line  shows  the 
result  obtained  for  an  accelerometer  of  reference 
standard  quality.  As  an  indication  of  the 
resolution  of  the  method.  Fig. 6(b)  is  a tracing 
of  the  computer  printout  for  the  latter. 

DISCUSSION  OF  UNCERTAINTIES 

Input  Transverse  Vibration 

The  motion  of  the  point  0 at  the  free  end 


FIGURE  6 Typical  Transverse  Calibration 
Results 


of  the  cantilever  actually  traces  out  a surface 
which  is  approximately  spherical,  rather  than  as 
idealized  in  the  xy  plane.  The  uncertainties 
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which  arise  from  the  consequent  non- recti linear- 
ity of  the  input  transverse  notion  are  discussed 
in  the  next  section. 

The  X and  Y accelerometers  which  are  used 
to  Measure  the  transverse  acceleration  are  then- 
selves  subjected  to  transverse  excitation. 
Uncertainty  froa  this  source  is  discussed  in 
detail  in  a later  section. 

The  evaluation  of  A,  Eq.fi),  is  influenced 
by  uncertainties  in  Ex,  Ey,  Sx  and  Sy.  To 
nininize  these  uncertainties  the  X and  Y 
acceleroneters  are  calibrated  k;  comparison  with 
a reference  acceleroneter,  Measuring  the 
electrical  output  ot  each  acceleroneter  with  the 
sane  digital  voltneter  that  is  used  with  it  on 
the  transverse  calibrator.  The  absolute 
calibration  of  the  voltneters  thereby  becoaes 
inconsequential. 

The  frequency  f is  easily  Measured  to 
within  0.1  Hz,  although  a less  precise  value 
would  be  acceptable  as  the  TSR  is  independent  of 
frequency  over  a wide  range  [2,31. 

Angle 

The  evaluatifm  of  the  direction  0 of  the 
input  transverse  vibration  is  influenced  by  the 
uncertainties  in  E»,  Ey,  Sx  and  Sy  (Eq.  (2)),  which 
have  been  discussed  above  . The  residual  un- 
certainty from  these  sources  is  mostly  attribut- 
able to  random  reading  errors. 

The  polar  setir  on  the  oscilloscope  is 
used  for  setting  ror.inal  values  of  6.  If  0 is 
read  directly  from  that  scale  the  additional 
error  does  not  exceed  2 deg. 

Some  further  uncertainty  is  incurred  in 
relating  the  test  value  of  0(TSRmin)  to  a 
reference  mark  on  the  body  of  the  test  accelero- 
neter. 


The  departure  of  the  X or  the  Y accelero- 
meter mounting  face  f.'om  squareness  with  the  Z 
accelerometer  mounting  face,  and  the  departure 
from  mutual  squareness  between  the  X and  Y 
mounting  faces  of  the  prototype  ca' tilever  beam 
are  all  less  than  two  minutes  of  arc.  The 
effects  of  these  small  discrepancies  can  be 
ignored. 

Output  of  Test  Accelerometer 

Uncertainty  in  the  measurement  of  the 
electrical  output  E^  resulting  from  the  trans- 
verse input  acceleration  A is  influenced  by 
signal-to-noise  ratio  and  by  harmonic  distortion, 
particularly  from  the  second  harmonic.  For 
improved  accuracy  a true-rms  voltmeter  is  used; 
this  is  explained  later  in  the  context  of  the 
discussion  of  non-rectilinearity  of  transverse 
acceleration. 

To  finu  the  electrical  out  p it  of  the  test 
accelerometer  resulting  from  vibration  of  the 
same  amplitude  A and  frequency  f applied  axially, 
the  accelerometer  is  calibrated  by  comparison 
with  the  same  reference  accelerometer  as  used  for 


the  determination  of  Sx  and  Sy,  measuring  the 
output  with  the  saMe  true-rms  meter  that  is  used 
in  the  transverse  calibration.  Only  relative 
values  of  Sy,  Sy  and  are  required  in  calculat- 
ing TSR  from  Eq.(l). 

Over All  Uncertainty 

For  an  accelerometer  with  axial  sensitiv- 
ity of  the  order  of  10  millivoi*s  per  g 
(where  g = 9.8  ms"1 ),  and  TSRmax  of  the  order  of 
one  per  cent,  the  uncertainty  in  TSR.-?ax  was 
estimated  to  be  i 2 per  cent  (of  TSRmcx) , and 
0(TSRminj  was  determined  with  an  uncertainty  f 
t 2 deg,  using  the  prototype  calibrator  sitn 
A * 200  ns"1 . 

Non-rectilinearity  of  Transverse 
Acceleration 

The  path  of  the  test  accelerometer  is 
assumed  to  be  a circular  arc,  as  shown  with 
greatly  exaggerated  curvature  in  Fig. 7. 


Z 


FIGURE  7 Non-rectilinearity  of  the 
Transverse  Motion 


Assuming  the  displacement (s)  along  the 
curved  line  to  be  represented  by  s = Ssinuit 
the  tangential  acceleration  (at)  is  given  by 

at  = - u>2  Ssinuit 

» Asinuit 

The  radial  acceleration  (a  ) acting 
along  the  sensitive  axis  of  the  test  accelero- 
meter is  given  by 

ar  * (uiScosuit)2  /R 

= u>2  S2  (cos2u>t  + 1)/2R 

where  R is  the  radius  of  curvature  of  the  path. 

The  ratio  of  radial  to  tangential 
oscillatory  acceleration  amplitude  is  therefore 
S/2R. 


If  the  deflected  form  of  the  beam  is 
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assumed  to  be  that  of  a uniformly  loaded 
cantilever  of  length  Lt  then,  at  maximimi  dis- 
placement S,  R can  be  shown  to  equal  3U/4. 

On  the  prototype,  the  value  of  R was  0.73  Li , 
the  difference  being  attributed  to  the  nun- 
uniform  loading  and  the  method  of  driving. 

As  an  example,  in  a transverse  calibration 
at  f - 100  Hi,  with  A - 200  ms*'  ,S»S»  10*4m. 
If  L,  is  0.4  m,  then  S/2R  * 0.0008.  That  is, 
the  oscillatory  component  of  ar  has  a peak  value 
which  is  6.08  per  cent  of  A,  and  is  of  frequency 
2f  * 200  Hi.  This  is  negligibly  small  in 
comparison  with  the  one  per  cent  limit  specified 
in  the  American  Standard  ( 2 J. 

This  effect  may  bj  measured  directly  if 
S(TStain)  is  first  . id;  with  this  direction 
of  transverse  excitation,  Ez  has  the  dominant 
component  at  frequency  2f,  as  well  as  random 
noise,  which  condition  is  easily  checked  on  the 
second  oscilloscope  screen  (Fig. 5).  If  the 
test  accelerometer  output  is  measured  through- 
out the  calibration  using  a true-rms  meter, 
then  all  E^  readings  can  be  corrected  for  these 
unwanted  components,  thus  improving  accuracy. 

The  axes  of  the  X and  Y accs.erometers 
are  displaced  a small  distance  4L  below  the 
mounting  surface  of  the  Z accelerometer  (Fig. 3). 
In  the  prototype  4L  * 8.7  bib,  thus  the  correct- 
ion factor  to  be  applied  to  A as  found  from 
Eq. (3)  is  1.030,  for  R - 0.302  m.  This  factor, 
the  ratio  of  the  transverse  vibration  amplitude 
at  the  Z mounting  surface  to  the  amplitude  at 
the  junction  of  the  Ox  axis  and  the  X mounting 
surface,  has  been  determined  from  direct  measure- 
ment  on  the  prototype,  using  an  image-shearing 
microscope.  The  uncertainty  in  the  determinat- 
ion does  not  exceed  0.5  per  cent. 

An  end  fixture  could  be  designed  to  make 
6L  « 0,  but  this  would  involve  increasing  the 
distance  of  the  X and  Y mounting  faces  from  the 
beam  axis,  which  is  undesirable. 

Curvature  of  the  path  also  produces  an  un- 
wanted tangential  component  due  to  tilting  of 
the  axis  of  the  accelerometer  in  the  earth's 
gravitational  field,  but  this'  effect  is 
negligible:  in  the  example  t - S/R  *0.0016 

(Fig. 7),  hence  the  tangential  component  from 
this  source  is  about  0.016  ms"’ , i.e.  0.008  per 
cent  of  at.  Likewise  the  change  (1  - cos$)  in 
the  radial  acceleration  is  negligible. 

Transverse  Excitation  of  the  X and  Y 
Accelerometers 

With  the  test  accelerometer  subjected  to  a 
transverse  simple  harmonic  motion  represented  by 
the  vector  X(8)  in  the  xy  plane  in  Fig. 3,  ideal 
X and  Y accelerometers  would  have  output 
voltages  proportional  to  Acos8  and  Asin8  respec- 
tively (see  Fig. 8).  Because  the  actual  X and  Y 
accelerometers  do  have  some  transverse 
sensitivity  the  outputs  are  respectively: 

Ex  = A Sjj(cos8  ♦ TSRySin8  isinut 


Ey  « A Sy(sin6  * TSRyCos6)sinut 

where  TSRX,  TSRy,  the  transverse  sensitivities 
of  the  X ard  Y accelerometers  respectively  in 
the  xy  plnne,  may  have  either  * or  - sign 
depe.'dii^  on  accelerometer  orientation. 


FIGURE  8 Transverse  Excitation  of 
Monitoring  Accelerometers 


If  the  sign  and  magnitude  of  TSR^  and 
TSRy  are  known,  a correction  may  be  applied 
when  calculating  A and  6 . For  the  case  where 
Sx  - Sy  and  | TSRy | « |TSRy|  = 0.01,  the 
correction  to  A is  $ 1.0  per  cent,  and  to 
8 < 0.6  deg. 

Alternatively  the  X and  Y accelerometers 
may  be  selected  for  low  TSPmax  (<  2 per  cent) 
and  oriented  so  that  their  8(TSRmin)  directions 
lie  as  nearly  as  practicable  in  the  xy  plane. 

If  they  are  ai'gned  in  this  way  to  within,  say, 

± 6 deg,  it  may  be  shown  that  the  error  due  to 
transverse  sensitivity  in  calculating  A does 
not  exceed  0.2  per  cent,  and  the  error  in 
calculating  8 does  not  exceed  0.1  deg.  This 
assumes  that  for  each  accelerometer  TSRmin  * 0. 
The  signs  of  TSRy  and  TSRy  are  indeterminate  in 
this  case,  hence  these  remaining  systematic 
errors  must  be  treated  as  further  sources  of 
uncertainty. 

The  effect  of  radial  acceleration  ar  on 
the  X and  Y accelerometers  is  negligible. 

CONCLUDING  REMARKS 

A new  method  of  transverse  calibration 
has  been  presented  which  allows  the  transverse 
response  properties  of  an  accelerometer  to  be 
evaluated  without  the  resetting  of  the  accelero- 
meter or  its  mounting  fixture  that  is  necessary 
in  methods  normally  used.  Experience  to  date 
in  using  the  prototype  to  calibrate  various 
piezoelectric  accelerometers  indicates  the  new 
method  to  be  accurate,  sensitive  and  easily 
applied. 

When  transverse  calibration  is  done 
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directly  on  an  electrodynamic  vibrator, 
unwanted  transverse  notion  of  the  vibrator  drive 
■ay  introduce  soae  difficulty:  with  the  new 
aethod  there  is  also  an  unwanted  component  of 
notion  but  this  is  of  small  aagnitude  and  is 
identifiable  and  Measurable  as  second  harmonic 
of  the  transverse  calibration  frequency, 

A Modification  of  the  calibrator  that 
would  appear  to  be  desirable  is  to  Mount  the 
Monitoring  accelerometers  on  a rotatable  fixture 
so  that  the  transverse  acceleration  in  any 
direction  is  Measured  with  an  accelerometer 
aligned  in  that  direction  and  hence  not  subjected 
to  transverse  excitation.  After  Making  tests 
with  such  an  arrangement  it  was  decided  to  adopt 
the  simpler  design  and  procedure  described  in 
the  paper. 

Further  developments  in  progress  will 
provide  improved  clamping  of  the  beam  and  will 
extend  the  frequency  range.  Obviously,  higher 
resonance . frequencies  can  be  attained  by 
changing  the  length,  material  and  cross  section 
of  the  beam.  Another  and  more  interesting 
possibility  is  to  operate  a cantilever  beam  at 
resonance  in  its  higher  modes.  Preliminary 
observations  with  the  prototype  system  driven 
at  its  resonance  frequency  (600  Hz)  in  its 
second  flexural  mode  have  confirmed  that  the 
range  can  be  extended  in  this  manner. 
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